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Abstract 

The  performance  of  biphase  direct-sequence  spread-spectrum  multiple- 
access  (SSMA)  communication  for  a  general  class  of  fading  channels  is 
investigated.  The  channels  considered  are  those  for  which  the  channel 
output  consists  of  a  strong  stable  specular  signal  plus  a  faded  version  of 
this  signal.  Such  channels  are  the  result  of  a  transmission  medium  which 
gives  rise  to  a  major  stable  communication  path  and  a  number  of  additional 
weaker  communication  paths.  The  fading  channel  is  modeled  as  a  general 
wide-sense-stationary  uncorrelated-scattering  (WSSUS)  channel  --  a  model 
which  is  general  enough  to  exhibit  both  time  and  frequency  selectivity 
and  to  impose  no  restrictions  on  the  fading  rate.  A  discussion  of  the 
important  parameters  of  the  WSSUS  channel  is  given  and  two  important 
classes  of  WSSUS  channels  are  developed  from  the  general  fading  channel 
model:  time-selective  fading  channels  and  frequency-selective  fading 
channels.  In  analyzing  the  performance  of  direct-sequence  SSMA  communica¬ 
tions  via  fading  channels  two  measures  of  system  performance  that  are 
considered  are  average  signal-to-noise  ratio  at  the  receiver  output  and 
the  average  probability  of  error. 


For  the  general  WSSUS  model,  results  are  obtained  for  the  average 
signal-to-noise  ratio  at  the  receiver  output  in  terms  of  the  spread- 
spectrum  signature  sequences  and  the  covariance  function  for  the  fading 
process.  The  results  are  then  specialized  to  both  time -selective  and 
frequency-selective  fading  channels.  For  these  two  classes  of  channels, 
expressions  are  obtained  for  the  correlation  receiver  output  signal-to- 
noise  ratio  in  terms  of  the  aperiodic  autocorrelation  functions  of  the 
signature  sequences,  the  covariance  function  of  the  fading  process,  and 
the  additive  white  Gaussian  noise  spectral  density.  Numerical  evaluations 
are  presented  for  specific  examples  of  each  of  these  two  types  of  channels. 
The  effects  of  fading  on  single-user  direct-sequence  spread -spectrum 
and  phase-shift-keyed  systems  is  discussed  and  it  is  shown  that  the 
spread-spectrum  system  yields  markedly  improved  performance  over  the 
phase-shift-keyed  system  when  frequency-selective  fading  channels  are 
considered.  Analytical  expressions  for  the  average  signal-to-noise  ratio 
are  derived  for  a  SSMA  system  with  random  signature  sequences,  and  the  use 
of  these  expressions  in  preliminary  system  design  is  discussed. 

An  expression  is  derived  for  the  probability  of  error  at  the  output 
of  a  correlation  receiver  in  terms  of  expectations  over  the 
faded  signal  and  the  multiple-access  interference.  Due  to  analytical 
difficulties  in  evaluating  the  expectations,  however,  a  moment  space 
bounding  technique  is  used  to  derive  bounds  on  the  probability  of  error. 

For  a  single  user  direct-sequence  spread-spectrum  system,  exact  expressions 
for  the  probability  of  error  are  given  for  both  time -selective  and 
frequency-selective  fading  channels.  These  expressions  are  compared 
with  the  results  obtained  from  Nth  moment  space  bounds. 
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CHAPTER  1 
INTRODUCTION 

1.1.  Spread-Spectrum  Multiple-Access  Communication  Systems 

The  advent  of  the  synchronous  communication  satellite  in  1963 
[Spilker,  1977]  has  made  possible  the  realization  of  practical  multiple- 
access  communication  systems.  This  realization  of  multiple-access 
communication  systems  has  been  achieved  principally  because  of  the  wide 
coverage  area,  small  duration  of  outage  times,  and  simplicity  of  antenna 
systems  afforded  by  the  use  of  synchronous  orbit  satellites.  To  date, 
three  major  types  of  multiple-access  schemes  are  either  being  utilized  by, 
or  are  being  proposed  for  use  with,  synchronous  satellites:  frequency- 
division  multiple-access,  time-division  multiple-access,  and  code-division 
multiple- access  [Pritchard,  1977].  Frequency-division  multiple-access 
(FDMA)  achieves  its  multiple-access  capability  through  the  use  of  a 
separate  carrier  frequency  for  each  user.  Its  principal  advantages  are 
that  it  is  compatible  with  existing  analog  trunk  systems  and  that  the 
satellite  design  for  a  FDMA  system  is  very  simple.  Time-division  multiple- 
access  (TDMA)  achieves  its  multiple-access  capability  through  the  use  of 
a  dedicated  time-slot  for  each  user.  The  principal  advantages  of  such 
a  system  are  that  it  fits  in  naturally  with  digital  communication  systems 
and  that  it  exhibits  increased  capacity  over  FDMA.  The  third  multiple- 
access  scheme,  code-division  multiple-access  (CDMA),  achieves  its 
multiple-access  capability  primarily  through  coding.  Unlike  FDMA  and 
TDMA,  however,  CDMA  requires  no  precise  time  or  frequency  coordination 


among  the  various  users.  CDMA  techniques  have  been  considered  for  use  in 
a  wide  variety  of  applications  including  the  NASA  tracking  and  data-relay 
system  (TDRS)  ([Stampfl  and  Jones,  1970],  [Chen  and  Burnett,  1977]), 
systems  to  provide  communication  to  aircraft  [Lebow,  et.  al.,  1971],  air 
traffic  control  systems  [Stiglitz,  1973],  systems  to  control  remotely- 
piloted  vehicles  [Malm  and  Schreder,  1973],  and  numerous  satellite 
communication  systems  (e.g.,  [Drouilhet  and  Bernstein,  1969],  [Kochevar, 
1977],  and  [Pritchard,  1978]).  In  fact,  a  recent  survey  article  of  29 
satellite  communication  systems  lists  five  proposed  military  satellite 
communication  systems  that  will  use  CDMA  [Pritchard,  1978]. 

The  most  common  form  of  CDMA  is  spread-spectrum  multiple-access  (SSMA) 
which  is  characterized  by  the  use  of  a  unique  code  sequence  assigned  to  each 
user  and  modulated  onto  the  carrier  along  with  the  digital  data.  By  a 
spread-spectrum  system,  we  mean  any  system  by  which  a  data  signal  is 
modulated  onto  a  wideband  carrier  so  that  the  resultant  transmitted  signal 
has  a  bandwidth  which  is  much  larger  than  the  data  signal  bandwidth 
[Scholtz,  1977].  The  reason  for  using  SSMA  is  that,  in  addition  to  its 
multiple-access  capabilities,  other  desirable  qualities  are  simultaneously 
provided  including  increased  immunity  from  interference  and  jamming,  low 
detectability,  compatability  with  other  (non- spread- spectrum)  systems 
operating  in  the  same  frequency  band,  and,  with  the  proper  choice  of 
modulation  method,  increased  immunity  from  the  effects  of  fading  and 
multipath  distortion  [Cahn,  1973] .  The  most  commonly  used  forms  of  SSMA 
are:  direct  sequence  SSMA  (DS/SSMA),  in  which  a  high  rate  code  is  used 
to  phase  modulate,  together  with  the  data  signal,  the  carrier-signal; 


frequency-hopped  SSMA  (FH/SSMA) ,  in  which  a  code  sequence  is  used  to 
control  a  frequency- synthesized  carrier  signal  onto  which  data  is 
modulated,  and  hybrid  SSMA,  which  can  be  a  combination  of  DS/SSMA  and 
FH/SSMA.  By  way  of  example,  DS/SSMA  will  be  used  in  the  NASA  TDRS  system 
[Chen  and  Burnett,  1977],  FH/SSMA  has  been  used  in  the  TATS  modulation 
system  for  the  Lincoln  Experimental  Satellites  [Drouilhet  and  Bernstein, 
1969],  and  hybrid  SSMA  will  be  used  in  the  proposed  Joint  Tactical 
Information  Distribution  System  (JTIDS),  a  military  communication, 
navigation,  and  identification  system  [Smith,  1978]. 

For  the  remainder  of  this  thesis  we  shall  be  concerned  only  with 
DS/SSMA,  also  known  in  the  literature  as  phase-coded  SSMA  [Pursley,  1974] 
or  pseudonoise  SSMA  [Anderson  and  Wintz,  1969].  In  particular,  we  shall 
be  concerned  with  analysis  of  the  communication  performance  of  an 
asynchronous  DS/SSMA  communication  system  operating  over  a  fading  channel. 
Such  a  system  has  been  previously  analyzed  elsewhere  for  communication 
over  additive  white  Gaussian  noise  (AWGN)  channels  (e.g.,  [Pursley,  1974], 
[Pursley,  1977],  [Yao,  1977]).  In  the  next  section  we  shall  present  a 
review  of  some  of  the  results  obtained  on  the  performance  of  DS/SSMA 
communications  via  AWGN  channels. 

1.2.  DS/SSMA  Communications  via  AWGN  Channels 

In  this  section  we  will  initially  review  the  results  obtained  by 
Pursley  (1974,  1979)  for  the  communications  performance  of  an  asynchronous 
DS/SSMA  communication  system  operating  over  an  additive  white  Gaussian 
noise  (AWGN)  channel.  The  purpose  of  this  review  is  threefold:  to  review 


the  principles  of  operation  of  a  DS/SSMA  communication  system,  to 
introduce  notation  that  will  be  useful  in  the  analysis  of  DS/SSMA 
communications  via  fading  channels ,  and  to  provide  a  basis  of  comparison 
of  performance  for  DS/SSMA  communications  via  fading  channels.  Later 
in  this  section  we  shall  present  a  brief  review  of  other  published 
literature  on  the  performance  of  DS/SSMA  systems  over  AWGN  channels. 

We  shall  consider  the  DS/SSMA  system  model  shown  in  Fig.  1  for  K. 
users.  The  k-th  users'  data  signal  b^(t)  is  a  sequence  of  unit  amplitude, 
positive  and  negative,  rectangular  pulses  of  duration  T,  given  by 

09 

bk(t)  =  s  \  *Mt_jiT)’  (I*1) 

2=-®  ’ 

where  b^  ^  €  {+l,-l}  denotes  the  k-th  user's  information  sequence  and 
PT(t)  =  1  for  0  <  t  <  T  and  pT(t)  =0  otherwise.  Each  user  is  assigned 
a  code  waveform  a^(t)  which  consists  of  a  periodic  sequence  of  unit 
amplitude,  positive  and  negative,  rectangular  pulses  of  duration  T  .  The 
code  waveform  for  the  k-th  user  may  therefore  be  written  as 


*k(t) 


-  Z  af°p  (MTJ. 

j=-»  J  c 


(1.2) 


(k) 

where  (a^  )  is  the  discrete  periodic  signature  sequence  assigned  to  the 

k-th  user.  We  assume  that  each  signature  sequence  has  period  N  =  T/T 


so  that  there  is  one  code  period  per  data  symbol. 


Figure  1.  DS/SSMA  system  model  with  AWGN  channel. 
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The  data  signal  bk(t)  m°dulated  onto  the  phase-coded  carrier 
c^(t) ,  which  is  given  by 

cfe(t)  =  Jrs  ak(t)cos(ojct  +  9k)  (1.3) 

so  that  the  transmitted  signal  for  the  k-th  user  is 

sk(t)  =  »/2P  ak(t)bk(t)cos(uoct  +  0^  .  (1.4) 

Here  P  represents  the  common  signal  power,  represents  the  common  center 
frequency,  and  ©k  represents  the  phase  of  the  k-th  carrier. 

For  asynchronous  systems,  the  received  signal  r(t)  at  the  input  to 
the  receiver  in  Figure  1  is  given  by 


K  _ 

r(t)  =  n(t)  +  Z  JT? 

k=l 


afc  ( t -T  k)  bk  ( t -T  k}  c  °s  (u)c  t-KPfc) 


(1.5) 


where  cpk  =  ®^"(JucTjc»  accounts  for  the  nominal  propagation  time  for  the 
k-th  signal,  and  n(t)  is  additive  white  Gaussian  noise  with  a  two-sided 
spectral  density  NQ/2. 

If  the  received  signal  r(t)  is  the  input  to  a  correlation  receiver 
matched  to  s^(t),  the  output  at  the  sample  moment  t  *  T  is  given  by 
[Pursley,  1977]  as 


X 

Z±  =  J  r(t)ai(t)cos  u^tdt 

“  v'WI  [b1>0x  +  ji[»k>.l»k>1<V(\,Afi<Vl'o' 

k^i 

T 

+  J  n(t)a, (t)cos  cju  tdt. 

0  c 


(1.6) 
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Since  we  are  concerned  with  relative  phase  shifts  modulo  2rr  and 
relative  time  delays  modulo  T,  we  have  assumed  in  (1.6),  without  loss  of 
generality,  that  9^  =  0  and  t  =  0.  Furthermore,  as  in  [Pursley,  1977], 
in  writing  (1.6),  we  have  neglected  the  double  frequency  components.  The 
two  terms  R^  ^  and  ^  appearing  in  (1.6)  are  the  continuous -time  partial 
cross-correlation  functions  of  the  code  waveforms,  defined  by 


Vi<T)  *  ■I',,  Vt-T)1i(c,dt 


(1.7) 


\  i(T)  =  J*  \(t-T)a.(t)dt 

’  T 


(1.8) 


for  0  <  t  <  T.  For  convenience,  we  shall  denote  the  continuous- time 
partial  autocorrelation  functions  R^  ^  and  R^  ^  by  R^  and  respectively. 
For  values  of  r  in  the  range  0  <  XT£  —  T  S  (i,+l)Tc  £s  T,  the  two  cross¬ 
correlation  functions  R^  ^  and  R^  .  can  be  written  as 


\,i(T)  =  Ck,i(X"N)Tc  +  tCk,i(A+1‘N)"Ck,iU“N)1(T‘jiTC)  (1*9> 


&k,i(T)  =  Ck,i<A)Tc  +  [Ck,ia+1)  ‘  Ck,i('4)1(T"^Tc)  (1'10) 


where  C,  .  is  the  discrete  aperiodic  cross-correlation  function  for  the 

tC  ,  3. 

sequences  (a^k^)  and  (a^^)  defined  by 


Ck,i<*> 


N’s“'ea(k)a(i) 

A  i 

N-l+£.  n  .  , . . 
Z  a(k)a(l) 
j=0  2  1  2 


0  <  L  <  N-l 


1-N  <  l  <  0 


(1.11) 


0 


Ul  2  N. 
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The  discrete  aperiodic  autocorrelation  function  C,  ,  will  be  denoted  by  C,  . 

Up  to  this  point,  we  have  not  indicated  any  means  of  measuring  the 
DS/SSMA  system  performance.  Three  performance  measures  that  have  been 
considered  in  the  past  are  average  signal-to-noise  ratio  [Pursley,  1974], 
[Pursley,  1977],  worst  case  performance  [Pursley,  1977],  and  average 
probability  of  error  [Yao,  1977].  For  DS/SSMA  systems  with  AWGN  channels, 
the  work  of  Yao  (1977),  coupled  with  the  results  of  Pursley  and  Sarwate  (1977b) 
on  the  efficient  computation  of  signal-to-noise  ratio  establishes  that  for 
systems  of  interest,  the  signal-to-noise  ratio  is  an  accurate  measure  of  per¬ 
formance  which  is  relatively  easy  to  compute.  Using  the  average  signal-to- 
noise  ratio  (SNR)  at  the  output  of  the  i-th  correlation  receiver  as  the 
measure  of  system  performance,  as  in  [Pursley,  1977],  the  phase  angles, 
time  delays,  and  data  symbols  for  the  k-th  signal  (k^i)  are  modeled  as 
mutually  independent  random  variables  which  are  uniformly  distributed  on 
[0,2tt],  [0,T],  and  {+l,-l},  respectively.  SNR.^  is  then  defined  by 

SNR.  A  Etzjb.  >0  =  +l}/(Var{Z.|bij0  =  +l})*,  (1.12) 


where,  without  loss  of  generality. 


we  have  assumed  b.  _  =  +1. 

1,0 


we  find  that 


From  (1.6) 


stzjbi  o  =  +1]  =>  Vp72  T 


(1.13) 


and 


VarlZilbl,0  -  +l5 


(T)  +  8i,i 


(T)]dT  +  kNQT 


K  N-l 

2  £ 


k-1  0 

k^i 


(X+1)T 

I 

IT 


'^,1<T)  +  ^,i<T>idT  +  w0T- 


(1.14) 


Using  (1.9)  and  (1.10)  in  (1.14),  (1.14)  reduces  to 


Var{Z£ 


+1} 


woT 


(1.15 


where 


and 


rk,i  *  2^k,i(0>  +  *k,i<l> 


(1.16 


N-l 


,(n)  =  2  Cki(1)Cki(l+n)l 

K>-L  1_n 


(1.17 


which  is  a  function  of  the  discrete  aperiodic  cross-correlation  functions 
^ .  In  [Purs ley  and  Sarwate,  1977b],  it  was  shown  that  ^(*0  can  be 
also  expressed  by 


N-l 

.(n)  =  L  C.  (i)  C  (X+n)  (1. 

f=l-N  K  1 


which  is  a  function  of  the  discrete  aperiodic  autocorrelation  functions  of 
the  code  signature  sequences.  From  (1.12),  (1.13),  and  (1.15),  the 
average  signal-to-noise  ratio  at  the  output  of  the  i-th  correlation 
receiver  is  found  to  be 

SNRi  * 

where  3  =  PT  is  the  signal  energy  per  bit. 


3  -1 
(6NJ) 


Nr 


,S.rk,i  +  23 
k=l 


(1.19 


k4i 


In  addition  to  the  results  obtained  by  Pursley  (1974,  1977)  indicated 


10 


above,  many  ocher  results  have  been  obtained  on  the  performance  of 
DS/SSMA  communications  via  AWGN  channels.  In  1977,  Yao  (1977)  presented 
several  approaches  for  evaluating  bounds  on  the  average  probability  of  error 
(P  )  of  a  DS/SSMA  system  using  an  isomorphism  theorem  from  the  theory  of 
moment  spaces.  Under  a  certain  set  of  conditions  (K  large,  N  large,  and 
N  »  K) ,  Yao  verified  that  Pg  could  be  approximated  quite  well  by  the  expres¬ 
sion  Pg  =  l-f(SNR^),  where  $  is  the  standard  Gaussian  cumulative  distribution 
function  and  SNR^  is  given  by  (1.19).  Further  discussion  of  Yao's  work 
may  be  found  in  Chapter  4  of  this  thesis.  Also  in  1977,  Roefs  and  Purs  ley 
(1977)  evaluated  the  performance  of  a  DS/SSMA  system  employing  binary 
random  sequences  as  signature  sequences .  They  found  that  the  expression 
obtained  for  the  average  signal- to-noise  ratio  using  random  binary  signa¬ 
ture  sequences, 


SNR 


K-l 

3N 


> 


) 


(1.20) 


was  a  very  accurate  approximation  to  SNR^  for  typical  system  values  of 
£?/Nq,  N,  and  K  and  thus  is  useful  in  the  preliminary  design  of  a  DS/SSMA 
system.  In  1978,  the  results  of  (Pursley,  1977]  were  extended  to  the 
quadriphase  DS/SSMA  system  ease  ([Pursley  and  Garber,  1978],  [Garber, 
1978a]).  Generally  speaking,  the  results  obtained  for  the  quadriphase 
DS/SSMA  system  case  are  very  similar  to  the  biphase  DS/SSMA  system  case 
with  the  exception  that  a  broader  class  of  signature  sequences  can  be 
employed  by  the  quadriphase  system,  offering  promise  of  improved  perfor¬ 
mance  over  AWGN  channels.  Other  results  that  have  been  obtained  on  the 
performance  of  DS/SSMA  communication  systems  over  AWGN  channels  deal 
primarily  with  properties  of  the  signature  sequence  correlation  functions 
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appearing  in  (1.19)  (e.g.,  see  [Roefs,  1977],  [Pursley  and  Sarwate, 
1977a],  [Pursley  and  Roefs,  1979],  [Sarwate  and  Pursley,  1979]).  A 
tutorial  treatment  of  many  of  these  properties  as  well  as  an  extensive 
bibliography  of  other  papers  treating  correlation  properties  of  signature 
sequences  may  be  found  in  [Sarwate  and  Pursley,  1979]. 

1.3  Direct  Sequence  Spread  Spectrum  Communications  via  Fading  Channels 

In  Section  1.1  we  noted  that  one  of  the  reasons  for  considering 
SSMA  as  a  multiple-access  technique  was  that,  with  the  proper  choice  of 
modulation  method,  SSMA  could  be  designed  to  reduce  the  effects  of  fading. 
In  fact,  one  of  the  earliest  applications  of  spread  spectrum  techniques 
was  in  communicating  over  multipath  channels  by  employing  a  pseudonoise 
modulated  frequency-shift  keyed  (FSK)  signal  [Price  and  Green,  1958].  In 
this  section  we  shall  review  some  of  the  literature  on  the  performance  of 
direct  sequence  spread-spectrum  (DS/SS)  communications  via  fading 
channe Is . 

One  of  the  earliest  studies  on  the  performance  of  DS/SSMA  over 
multipath  channels  was  one  conducted  by  Massey  and  Uhran  (1969)  for  NASA 
in  1969.  This  study  was  concerned  with  the  design  of  good  sets  of 
signature  sequences  which  would  permit  reliable  communication  over 
multipath  channels.  Essentially  this  study  culminated  in  the  realization 
that  the  odd  cross-correlation  function  of  the  signature  sequences, 
defined  in  [Pursley,  1977],  is  just  as  important  as  the  periodic  cross¬ 
correlation  function  for  communication  over  multipath  channels  or  for 
asynchronous  DS/SSMA  communications.  Some  of  these  results  have  also 
been  reported  in  [Massey  and  Uhran,  1975].  Kadar  and  Schreiber 


(1971)  considered  the  effects  of  a  single  fading  path  together  with  a 
single  nonfading  path  on  the  performance  of  a  DS/SS  system  employing  a 
post-detection  correlation  receiver.  Their  primary  result  was  an  expres¬ 
sion  for  the  signal- to-interference  ratio,  neglecting  thermal  noise,  at 
the  output  of  the  post-detection  correlation  receiver.  Cahn  (1973)  presented 
"heuristic”  analyses  of  the  effects  of  direct-plus-specular  multipath  and 
direct-plus-diffuse  multipath  on  the  performance  of  a  single-user  DS/SS 
communication  system.  Cahn  assumed  that  the  fading  components  of  the 
received  signal  were  Gaussian  and  hence  could  be  treated  as  additional 
noise  in  deriving  an  expression  for  the  probability  of  error  at  the  output 
of  a  correlation  receiver.  Jacobs  (1974),  in  a  tutorial  paper  considered 
the  effects  of  a  direct-plus -specular  multipath  channel  model  on  a  DS/SS 
system.  Finally,  Chang  (1979),  considered  the  effects  of  fading  on  the 
performance  of  a  digital  matched  filter  in  a  communication  system  employing 
amplitude  shift  keying  (ASK)  together  with  a  maximal -length  shift  register 
signature  sequence. 

Within  the  past  several  years,  there  have  appeared  several  studies 
on  the  effects  of  fading  on  DS/SSMA  cotmminication  systems  ([Orr,  1977], 
[Welch,  1978],  [Gardner  and  Orr,  1979]).  All  of  these  studies  have 
concentrated  on  determining  the  system  signal-to-noise  ratio  and/or 
probability  of  error  for  multiple-user  DS/SS  systems  communicating  via 
a  class  of  channels  known  as  nondispersive,  slow-fading  channels  (see 
Chapter  2) . 

The  literature  survey  given  above  is  intended  to  indicate  the 
amount  and  direction  of  research  that  has  already  been  performed  on 
DS/SS  communications  via  fading  channels.  It  is  also  intended  to 
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indicate  some  of  the  deficiencies  in  the  analysis  of  DS/SSMA  communications 
over  fading  channels.  In  particular,  in  all  of  the  analyses  presented  to 
date,  little,  if  any  attention  has  been  given  to  detailed  consideration  of 
the  fading  channel  model  being  used.  While  the  assumptions  of  direct- 
plus-specular  fading  or  nondispersive  fading  may  be  valid  for  some  physical 
channels,  the  channels  over  which  spread-spectrum  systems  are  employed 
often  exhibit  radically  different  channel  characteristics.  For  example, 
because  the  bandwidth  occupied  by  a  DS/SSMA  system  is  typically  quite 
large,  the  channel  used  may  exhibit  frequency  selectivity  due  to  the 
presence  of  a  distribution  of  channel  paths  over  which  the  signal  propa¬ 
gates.  Clearly,  such  a  characteristic  is  not  incorporated  in  either  the 
direct-plus-specular  channel  model  or  the  nondispersive  fading  channel 
models.  Furthermore,  the  published  results  noted  above  have  either  made 
certain  very  specialized  assumptions  in  deriving  the  results  above  or 
have  not  used  correlation  receivers  in  the  DS/SSMA  system  model.  In 
this  thesis  we  will  analyze  the  performance  of  the  DS/SSMA  system  model 
shown  in  Figure  1  for  a  broad  class  of  fading  channels.  In  the  next 
section  an  outline  of  the  thesis  research  is  given. 

1.4  Outline  of  the  Thesis 


Chapter  2  describes  a  comprehensive  model  of  a  fading  channel 
in  terms  of  the  statistical  properties  of  the  channel  with  the  viewpoint 
of  characterizing  typical  fading  channels  over  which  DS/SSMA  systems 
might  be  used.  This  characterization  includes  development  of  the  impor¬ 
tant  wide-sense-stationary  uncorrelated-scattering  (WSSUS)  channel  model 
and  three  subclasses  of  the  WSSUS  channel  model:  the  time-selective 


fading  channel  model,  the  frequency-selective  fading  channel  model,  and 
the  nondispersive  fading  channel  model. 

Chapter  3  analyzes  the  performance  of  DS/SSMA  communications  via 
fading  channels  using  the  average  signai-to-noise  ratio  at  the  output  of 
a  correlation  receiver  as  the  system  performance  measure.  In  this  analysis, 
the  DS/SSMA  system  model  of  Chapter  1  is  used  together  with  the  WSSUS  fading 
channel  model  of  Chapter  2.  After  analysis  of  system  performance  over 
doubly-spread  channels  is  considered,  performance  of  the  DS/SSMA  system 
over  time- selective  and  frequency-selective  fading  channels  is  analyzed. 

A  comparison  of  the  performance  of  single-user  phase-shift  keyed  (PSK)  and 
DS/SS  communication  systems  over  fading  channels  is  given.  An  analysis  of 
the  performance  of  a  DS/SSMA  system  using  random  binary  sequences  as  the 
signature  sequences  is  also  given. 

Chapter  4  analyzes  the  performance  of  DS/SSMA  communications  via 
fading  channels  using  probability  of  error  at  the  output  of  a  correlation 
receiver  as  the  system  performance  measure.  Bounds  on  the  probability  of 
error  are  obtained  using  an  isomorphism  theorem  from  the  theory  of  moment 
spaces.  The  "Gaussian  approximation"  for  the  probability  of  error, 

P  »  1  -  ^(SNR^),  for  DS/SSMA  communications  via  fading  channels  is  dis¬ 
cussed. 

In  Chapter  5,  a  summary  of  the  major  results  of  this  thesis  is  given. 
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CHAPTER  2 

INTRODUCTION  TO  FADING  CHANNELS 
2.1.  A  Description  of  Fading 

In  many  analyses  of  the  performance  of  various  types  of  communication 
systems,  the  communication  channel  is  frequently  modeled  as  a  linear  time- 
invariant  system  whose  transfer  function  consists  of  a  frequency-independent 
magnitude  less  than  unity  proportional  to  the  propagation  loss  and  a 
delay  term  proportional  to  the  propagation  delay  between  the  channel 
modulator  and  the  channel  demodulator.  In  addition,  the  channel  is  usually 
considered  to  be  corrupted  by  additive  stationary  Gaussian  noise.  While 
this  simple  additive  white  Gaussian  noise  model  is  quite  accurate  for 
channels  such  as  deep  space  communication  channels,  it  is  often  an  overly- 
simplified  model  for  high-frequency  (HF)  long  distance  communications 
achieved  via  the  ionosphere  and  for  microwave  communications  beyond  the 
horizon  achieved  through  the  use  of  tropospheric  scatter.  In  the  latter 
two  channels,  the  received  signal  has  been  experimentally  shown  to  undergo 
a  process  known  as  fading.  Fading  Is  a  term  generally  used  to  describe 
any  linear  channel  whose  performance  is  other  than  that  of  the  ideal 
channel  described  above.  A  fading  channel  may  exhibit  such  properties  as 
selective  frequency  response,  intersymbol  interference  in  digital 
communications,  spreading  of  signals  in  the  frequency  domain,  a  time- 
varying  amplitude  response,  or  any  combination  of  these  attributes.  The 
description  of  a  comprehensive  model  for  a  fading  channel  is  given  in  this 
chapter.  The  remainder  of  this  section  will  describe  fading  from  a 
phenomenological  point  of  view  and  give  special  examples  of  fading 


channe Is . 
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In  addition  to  the  two  examples  of  fading  channels  given  above,  the 
HF  ionospheric  channel  and  the  microwave  tropospheric  scatter  channel, 
there  exist  several  other  types  of  channels  which  exhibit  fading.  One 
example  is  a  very  high  frequency  (VHF)  communication  link  between  an 
aircraft  and  a  synchronous  satellite  relay  [Bond  and  Meyer,  1966].  Such 
a  link  has  been  analyzed  and  has  been  shown  to  exhibit  fading  due  to  the 
presence  of  a  secondary  propagation  path  between  the  satellite,  the 
earth's  surface,  and  the  airplane,  in  addition  to  the  primary  propagation 
path  between  the  airplane  and  the  satellite.  Another  example  of  a 
fading  channel,  this  time  an  artificially  created  one,  is  the  communica¬ 
tions  channel  temporarily  created  by  an  experiment  known  as  the  West 
Ford  Project  [Lebow,  et.  al.,  1964].  In  this  project,  20  Kg  of  2  cm 
long  copper  dipoles  were  injected  into  an  orbit  about  the  earth  and 
transhorizon  communications  were  conducted  at  8GHZ  using  this  orbiting 
dipole  belt  as  a  scattering  mechanism.  Fading  was  predicted  prior  to 
the  experiment  and  confirmed  experimentally  during  the  course  of  the 
experiment.  A  third  example  of  a  fading  channel  is  that  of  line-of- 
sight  microwave  communication  links.  Although  such  links  are  designed 
to  provide  reliable  point-to-point  communications,  occasionally  these 
links  undergo  severe  fading  due  to  the  formation  of  tropospheric  inversion 
layers  permitting  multiple  transmission  paths  between  the  transmitting 
and  receiving  antennas  [Jakes,  1978].  As  a  result  of  these  multiple 
paths,  the  received  signal  is  a  vector  sum  of  several  delayed  versions 
of  the  transmitted  signal,  and  consequently,  depending  on  the  phase 
relationship  of  the  received  multiple  signals,  the  channel  exhibits 


fading.  Recently,  Rummler  (1978,  1979)  has  shown  that  the  fading  effects 
exhibited  by  line-of-sight  microwave  channels  could  be  modeled  well  by  a 
three-path  channel  model.  A  final  example  of  another  channel  exhibiting 
fading,  and  one  of  much  current  interest,  is  that  of  communication  at* 
millimeter  to  optical  wavelengths  in  line-of-sight  paths  through  the 
non-ionized  atmosphere  [Strohbehn,  1968].  In  this  channel,  fading  is 
present  due  to  random  fluctuations  in  the  dielectric  constant  of  the 
atmosphere.  These  four  examples  serve  to  demonstrate  that  fading  is  not 
limited  to  just  the  two  "classical"  fading  channels,  the  HF  ionospheric 
and  microwave  tropospheric  scatter  channels. 

Fading  encountered  over  an  HF  ionospheric  channel,  for  example,  has 
been  experimentally  verified  to  be  of  two  types  —  short  duration  rapid 
fading  over  time  spans  less  than  a  second  and  long  duration  slow  fading 
over  time  spans  from  one  second  to  an  hour  or  longer.  The  statistics 
of  the  two  fading  processes  are  different;  hence  these  two  types  of  fading 
must  be  accounted  for  in  the  channel  model.  It  will  be  shown  in  Section  2.3 
that  the  two  types  of  fading  lead  to  a  "quasi-wide- sense- stationary" 
fading  channel  model  for  practical  radio  channels. 

The  origin  of  the  fading  mechanism  for  most  of  the  fading  channels 
mentioned  above  may  be  traced  to  the  scattering  of  an  electromagnetic 
wave  off  a  random  medium.  To  see  how  this  leads  to  fading,  consider  the 
following:  let  a  single  continuous  sine  wave  be  allowed  to  be  scattered 
by  a  random  medium.  The  scattered  components  may  be  resolved  into 
in-phase  and  quadrature  components.  The  instantaneous  values  of  the  two 
types  of  components  may  be  shown  to  be  uncorrelated.  Using  the  central 
limit  theorem,  as  the  number  of  in-phase  and  quadrature  components  becomes 
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large,  the  sum  of  the  in-phase  components  approaches  a  Gaussian  random 
process.  Similarly,  the  quadrature  components  add  to  form  an  identically 
distributed  Gaussian  random  process.  Hence,  the  in-phase  and  quadrature 
random  processes  collectively  form  a  zero-mean  complex  Gaussian  random 
process  [Van  Trees,  Sec.  A. 3.1,  1971].  If  the  random  medium  is  a  single 
surface  and  is  time-invariant,  the  received  signal,  after  scattering,  can 
be  shown  to  have  a  Rayleigh-distributed  amplitude  and  uniformly  distributed 
phase,  i.e.,  the  signal  is  undergoing  fading.  Further  discussion  of  the 
modelling  of  scattering  of  sinusoids  off  of  random  media  as  a  complex 
Gaussian  random  process  may  be  found  in  [Stein,  1966];  in  the  following 
we  will  accept  the  fact  that  scattering  results  in  a  zero-mean  complex 
Gaussian  random  process,  provided  that  a  "sufficient  number"  of  random 
scatterers  exists  for  the  particular  geometry  under  consideration. 

2.2.  Statistical  Model  of  Fading 

In  this  section,  we  develop  the  most  general  model  of  the  fading 
channel  with  the  idea  in  mind  that  the  model  should  be  applicable  to  the 
analysis  of  DS/SSMA  communications  over  fading  channels.  In  subsequent 
sections,  the  model  developed  here  will  be  simplified  by  making  appropriate 
assumptions . 

To  begin  with,  since  the  signals  under  consideration  are  bandpass  signals 
(i.e.,  narrow-band  signals  centered  at  some  frequency  ua^) ,  complex 
envelope  notation  is  the  most  convenient  method  of  describing  these 
signals.  To  clarify  the  notation  to  be  used,  several  standard  properties 
([Van  Trees,  1971],  [Stein,  1966],  [Stein  and  Jones,  1967])  of  bandpass 
signals  are  given  below  in  complex  envelope  notation. 


The  transmitted  signal  may  be  represented  in  complex  envelope 


notation  as 


s0(t)  =  Re[u0(t)exp(j2Trf0t)  ] ,  (2.1 

where  u^(t)  is  a  lowpass  signal  having  a  Fourier  transform  U^Cf), 

00 

uQ(t)  =  J  UQ(f)exp(j2  n ft)df .  (2.2 
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By  narrowband  signals,  we  mean  that  if  we  define  the  normalized  bandwidth 


of  uQ(t)  as 


J>0(f)|2df 


Bn  - - " - Y 

|u0(0)| 


B  v<  f 
n  0 


Given  a  linear,  time-invariant  system  with  an  impulse  response  h(t) 
and  a  transfer  function  H(f),  where 
00 

h(t)  -  J*  H(f)e4'j2  TTftdf  =  Re(2h  (t)exp(j2TTf  t)}  (2.5) 


and  H(f)  is  a  bandpass  function  around  fg,  the  output  y(t)  due  to  an 


input  Sg(t)  given  by  (2.1)  is 


where 


y(t)  »  Re 
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Equations  (2.6)  and  (2.7)  illustrate  the  advantage  of  using  complex 
envelope  notation  for  bandpass  signals  and  systems;  to  evaluate  the  output 
of  a  bandpass  system  due  to  a  bandpass  signal  input,  we  simply  convolve 
the  impulse  response  envelope  with  the  input  signal  envelope,  multiply  by 
exp(j2irfgt)  and  take  the  real  part  of  the  resulting  product. 


For  a  linear  time-varying  bandpass  system  with  impulse  response  h(t,t) 


where 


h(t,T)  =  Re[2h£(t,T)exp(j2  tt(  t-T))], 


(2.8) 


the  expression  corresponding  to  (2.7)  is 


yg(t)  =  J  he(t,u)uQ(u)du. 


(2.9) 


In  (2.8),  h(t,T)  denotes  the  output  at  time  t  due  to  an  impulse  at  time  t. 
Finally,  we  will  need  the  following  algebraic  identity:  Given  two 


complex  numbers  X  and  Y,  then 


Re[X]Re[Y]  =  ^  Re[XY]  +  %  Re[XY  ] 


(2.10) 


Other  properties  of  complex  envelope  notation  will  be  developed  as  needed. 

We  now  proceed  to  the  development  of  the  fading  channel  model. 
Consider  the  propagation  model  shown  in  Figure  2.  In  the  figure,  a  signal 
Sg(t)  given  by  (2.1)  is  transmitted  and  scattered  by  a  moving  random 
medium  which  is  assumed  to  be  able  to  be  modeled  as  a  layered  scatterer, 
where  each  layer  has  an  incremental  thickness  dr.  Associated  with  each 
layer  is  a  propagation  delay  t  which  is  in  addition  to  the  nominal 
propagation  delay  t^  between  the  transmitter  and  the  receiver.  The 
received  signal  due  to  scattering  from  the  layer  whose  incremental 


propagation  delay  is  t  is  given  by 


s(t,T)  =  Re[P  (t  ,  t-t0)u0(t-tQ- t)  exp(j2  tt  fQ(t-t0- t  ))  ]dT  (2.11) 

where  |3  (r ,  t)  is  the  time-varying  transmission  coefficient  for  a  wave 
scattering  of  a  layer  whose  incremental  propagation  delay  is  t.  Several 
comments  about  the  form  of  (2.11)  are  in  order.  First,  (2.11)  assumes  that 
the  thickness  dr  and  its  corresponding  propagation  delay  dr  may  be  chosen 
small  enough  such  that  j3(x,t)  is  constant  in  its  first  argument  over  the 
delay  dr.  This  allows  the  effect  of  scattering  off  a  single  layer  to  be 
modeled  as  a  simple  time-varying  multiplicative  factor  as  opposed  to  some 
form  of  superposition  of  responses.  Second,  it  is  important  to  note  that 
the  received  signal  s(t,T)  is  a  function  of  two  arguments:  the  first 
denotes  time  and  the  second  denotes  the  incremental  propagation  delay 
due  to  the  layer  that  the  transmitted  signal  is  being  scattered  from. 

Hence  the  received  signal  is  both  time-varying  due  to  the  fact  that  the 
scattering  volume  was  assumed  to  be  moving,  and  layer  dependent  since  the 
scattering  medium  was  assumed  to  be  able  to  be  modelled  as  consisting  of 
differential  layers.  Third,  the  time  origins  of  the  various  terms  of 
(2.11)  should  be  noted.  For  a  received  signal  at  time  t  due  to  scattering 
from  a  layer  having  a  total  propagational  delay  of  t^  +  t  associated  with 
it,  the  signal  would  have  to  be  transmitted  at  time  t-t^-x  .  Thus  the 
arguments  of  Uq(’)  and  exp (j2  tt  f  ( •) )  are  the  transmitted  time  origin. 

The  second  argument  of  j3(T,t-tQ)  refers  to  the  nominal  time  at  which 
the  received  signal  was  transmitted.  Such  a  time  origin  is  by  no  means 
unique;  other  authors  have  made  the  time  origin  of  the  scattering  process 


23 


the  actual  time  origin  of  the  transmitted  signal  t-tg-T  [Kailath,  1961] 
or  the  time  at  which  the  transmitted  wave  was  scattered  off  a  reasonably 
well-defined  medium  [Van  Trees,  Ch.  13,  1971].  All  such  time  origins 
are  essentially  equivalent,  although  one  choice  of  time  origin  may  be 
more  appealing  to  the  intuition  for  a  given  type  of  scattering  medium. 
Finally,  the  statistical  nature  of  0(T,t)  must  be  determined.  From  our 
discussion  in  Section  2.1,  we  assume  0(T,t)  to  be  a  zero-mean  complex 
Gaussian  random  process.  Inherent  in  this  assumption  is  the  fact  that 
j3(T,t)  has  a  Rayleigh- fading  envelope  and  a  uniformly  distributed  phase, 
both  of  which  are  time  varying  due  to  the  presence  of  a  moving  random 
medium.  Hence,  the  term  exp(j2TrfgT)  may  be  absorbed  into  0(T,t)  in  (2.11). 

The  total  received  signal  is  a  superposition  of  the  responses  due 
to  all  the  scattered  layers : 


OB  OB  -I 

s(t)  *  J  s(t,T)dT  =  Re  J  p  (T,t-t0)u0(t-tQ-T)exp[j2  TT  f0(t-tQ)  ]dr  .  (2.12) 
-00  -00  J 

Equation  (2.12)  may  also  be  written  in  the  form 


s(t)  »  Re[u(t-tg)exp  j  (2  tt  fQ (t- tQ) )  ] 


(2.13) 


where 


u(t)  =  J  0  (T,  t)uQ(t-T)dT  =  J  0  (t-T  ,t)uQ(T)dT 


(2.14) 


Comparing  (2.13)  and  (2.14)  with  (2.6)  and  (2.9),  it  is  readily 
seen  that  0(T,t)  represents  a  time-varying  equivalent  low-pass  impulse 
response  for  the  general  fading  channel.  Applying  an  impulse  at  time  a 
to  (2.14),  we  obtain 
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Uj(t)  =J  P(T,t)  6(t-a-T)dT  =  P(t-Of,t).  (2.15) 

-CD 

That  is,  3(t-Q(,t)  is  the  response  at  time  t  to  an  impulse  applied  at 
time  Qf . 

Equation  (2.12)  represents  the  most  general  model  of  a  fading  channel 
that  we  will  discuss  here.  Before  proceeding  further,  experimental 
evidence  justifying  the  use  of  this  model  as  well  as  suggesting  other 
models  will  be  examined  briefly.  To  begin  with,  if  the  random  medium  is 
slowly  varying  with  time  and  can  be  modeled  as  consisting  of  a  single 


layer,  (2.12)  reduces  to 


s(t)  =  Re[puQ(t-t0)exp  j 2  tt  f Q (t- tQ)  ] 


(2.16 


where  P  is  a  zero-mean  complex  Gaussian  random  variable.  The  assumptions 
necessary  for  the  transition  from  (2.12)  to  (2.16)  will  be  examined  in 
greater  detail  in  Section  2.7.  Letting  the  transmitted  signal  be  a 
continuous  sine  wave  at  a  frequency  f^,  (2.16)  becomes 


s(t)  =  Re [P exp  j2  tt  fQ(t-t0)  ] 


(2.17 


It  can  be  shown  that  the  amplitude  of  a  zero-mean  complex  Gaussian  random 
variable  is  Rayleigh-distributed  and  the  phase  uniformly  distributed. 
Hence  (2.17)  predicts  that  the  amplitude  of  the  received  sine  wave  is 
Rayleigh  distributed.  Experimental  evidence  (e.g.  [McNicol,  1949], 
[Grisdale,  et.al.,  1957])  for  the  HF  ionospheric  channel  tends  to  support 
this  prediction;  at  the  same  time,  experimental  data  [McNicol,  1949]  for 
many  fading  channels  exhibit  a  Rician  amplitude  distribution.  Rician 
fading  is  usually  considered  to  be  due  to  a  specular  component  in  the 


received  path,  i.e.,  the  presence  of  a  fixed  non-random  scatterer  in  the 
transmitter-receiver  propagation  path.  In  this  case  (2.12)  may  be 
modified  to  include  a  deterministic  component  in  the  received  signal: 

CO 

s (t)  =  Re[ (Au(t-tQ)  +  J  P(T,t-t0)uQ(t-t0-T))dT  exp(j2 n f0(t-tQ) ]  (2.18) 

—CO 

where  A  is  the  transmission  coefficient  associated  with  the  specular 
path.  The  presence  of  a  specular  component  is  necessary  for  coherent 
communications  including  DS/SSMA  communication  systems;  accordingly 
(2.18)  will  be  used  in  subsequent  work. 

The  above  comments  are  not  meant  to  suggest  that  all  fading  is 
Rician  or  Rayleigh  in  nature;  on  the  contrary,  experimental  data  for 
the  HF  ionospheric  channel  and  the  mm-to-optical  wave  line-of-sight 
channel  indicate  that  the  amplitude  of  the  received  signal  is  lognormally 
distributed  in  many  instances.  Stein  (1966)  suggests  that  deviations  from 
a  Rayleigh  amplitude  distribution  for  the  HF  ionospheric  channel  may  be 
due  to  an  insufficient  number  of  scatterers  in  a  scattering  "layer" 
to  support  the  complex  Gaussian  random  process  hypothesis  stated  in 
Section  2.1.  For  the  line-of-sight  mm-to-optical  wave  channel, 

Strohbehn  (1968)  has  shown  that  the  channel  may  be  modeled  as  consisting 
of  a  very  large  number  of  randomly  moving  slabs.  The  transmitted  signal 
is  then  modulated  by  the  product  of  the  random  transmission  coefficients 
of  each  of  the  slabs.  Assuming  the  transmission  coefficients  to  be 
independent  of  each  other,  the  logarithm  of  the  overall  channel  amplitude 
transmission  function  is  the  sum  of  the  logarithms  of  the  slab  transmission 
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functions  and,  using  a  central  limit  theorem  argument,  is  thus  normally 
distributed.  The  channel  transmission  function  is  then  seen  to  be 
lognormally  distributed  in  amplitude.  Despite  the  fact  that  some  channels 
may  be  modelled  as  exhibiting  a  lognormal  amplitude  distribution,  for  the 
remainder  of  this  study,  we  assume  that  the  fading  channels  of  concern 
exhibit  predominately  Rician  fading. 

In  determining  the  performance  of  communication  systems  over  fading 
channels,  frequent  use  is  made  of  the  second  order  statistical  properties 
of  the  channel  model.  In  particular,  assuming  |3 (t , t)  to  be  a  zero-mean 
complex  Gaussian  process,  given  Ug(t),  u(t)  given  by  (2.14)  is  also  a 
zero-mean  complex  Gaussian  process  which  is  completely  characterized 
statistically  by  its  covariance  function.  What  is  not  commonly  realized 
however,  is  that  for  a  bandpass  process  represented  by  complex  envelope 
notation,  two  covariance  functions  are  needed  to  completely  characterize 
the  second  order  properties  of  the  process.  To  see  this,  the  covariance 
of  s(t)  given  by  (2.12)  is  defined  in  the  usual  manner  as 

Rs(ti,t2)  »  EU(t1)s(t2)]  .  (2.19) 

Using  (2.12)  and  (2.10)  in  (2.19),  we  obtain 

['  OB  ® 

J*  I  ^E^(Ti»ti"to^  ^T2,t2"t0^U0^tl't0”Tl^uO^t2’C0~T2^ 

-00  .00 

- 

*  exp  [  j  rrf0(t1+t2-2tQ)]dT1dT2 

t 

r  00  CD  ^ 

+  Re-J*  J  ^^(Ti»trto)P*(T2't2'tO)}uO(trVTl)uO(t2'tO-T2) 

^-09  _S0 

•exp[  j2  TT  fQ(t1-t2)]dT1dT21  • 


(2.20) 
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Note  that  to  evaluate  (2.20),  two  covariance  functions  are  needed, 

^l’VW  =  %E[P(T1,t1)P(T2,t2)]  (2.21) 

and 

A(Tl*T2;tl*t2)  =  (T1,t1)P*(T2,t2)}.  (2.22) 

A(*)  given  by  (2.22)  is  defined  to  be  the  space-time  cross -covariance 
function  (or  simply  the  covariance  function)  of  the  fading  process 
[Stein,  1966].  In  most  applications,  the  narrow-band  process  P(T,t)  is 
so  constituted  that  ([Bello,  1963],  [Van  Trees,  1971]) 

A(Tl’T2;Cl,t2)  ^  °*  (2.23) 

Examples  of  processes  which  do  not  satisfy  (2.23)  are  given  in 
([Bello,  1961],  [Brown  and  Crane,  1969]).  It  is  easily  seen  that  (2.23) 
is  a  necessary  condition  for  stationary  bandpass  processes.  For  by  representing 
the  equivalent  bandpass  impulse  response  in  complex  envelope  notation  as 

B  (t  ,  t)  =  ReCp(T,t)exp(j2nfQt)}  (2.24) 

as  in  (2.8),  the  autocorrelation  function  of  B(T,t)  is  given  by 

fVTl,tl;T2,t:2^  =  Re[%ECP(T1>t1)P(T2»t2^exp^2TTf0^,:l  +  t2^'* 

+  Re[]jE{p  (T1,t1)P*(r2,t2)]exp(j2  n  fQ(t;L  -  t2))] 

=  Re[A(T1,72;t1,t2)exp(j2  rr  fQ  (tx  +  t2))] 

+  Re[A(T  .T.jt,  ,t  )exp(j2  tt  f  (t.  -  t_))].  (2.25) 
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In  order  for  R  (•)  to  be  a  function  only  of  the  time  difference  t1  -  t„, 

A(*)  must  vanish  in  (2.25).  A  different  proof  of  this  result  utilizing 
the  spectra  of  the  bandpass  processes  may  be  found  in  [Van  Trees,  Sec.  A. 3.1, 
1971].  Grettenberg  (1965)  has  proven  that  a  necessary  and  sufficient 
condition  for  (2.23)  to  hold  is  that  P(T,t)  and  eJ  (3  (t  ,  t)  be  identically 
distributed  for  all  real  9.  A  more  direct  proof  of  this  latter  result  may 
be  found  in  [Miller,  Th.  II. 4. 4,  1974]. 

In  the  sequel,  it  will  be  assumed  that  (2.23)  holds  for  all  complex 
envelope  processes  that  we  shall  encounter.  There  are  three  reasons  for 
this:  First,  in  virtually  all  the  literature  on  fading  channels  (e.g., 

[Bello,  1963],  [Kennedy,  1969],  [Stein,  1966],  [Van  Trees,  1971])  this 
assumption  is  made  either  explicitly  or  by  modelling  the  equivalent 
bandpass  impulse  response  as  a  stationary  process.  Second,  if  A(*)  is 
non-zero,  much  of  the  appeal  of  the  representation  of  bandpass  processes 
by  complex  envelope  notation  is  lost  since  two  covariance  processes  must 
now  be  specified  whereas  the  real  bandpass  process  requires  only  one. 

Finally,  and  probably  the  most  important  reason  from  a  practical  standpoint 
is  the  fact  that  even  if  A(*)  were  non-zero,  it  is  doubtful  that  a 
meaningful  measurement  of  A(’)  could  be  made  over  an  actual  fading 
channel,  considering  the  difficulty  in  measuring  A(»)  [Bello,  1971]. 

One  other  covariance  function  associated  with  the  general  fading 
channel  model  given  by  (2.12)  is  the  frequency-time  cross-covariance 
[Stein,  1966]  defined  by 

RF(f1,f2;t1,t2)  =  ^E[H(f1;t1)H*(f2;t2)}  (2.26) 

where 

00 

H(f;t)  *  J  (3  (t  ;  t)exp(-j2  rr  f  t  )dT 


(2.27) 


is  the  range  Fourier  transform  of  the  equivalent  low-pass  impulse  response 
P (t , t)  for  the  channel.  Using  (2.27)  in  (2.26)  we  obtain 

00  00 

Rp(f1,f2;t1,t2)  =  J  J  A(t  1,T2;t1>t2)exp(j2  tt  (f1T1-f2r2))dT1dT2  (2.28 

—00  —00 

which  is  recognized  to  be  the  double  Fourier  transform  of  the  space- time 
cross-covariance  function. 

In  practice  the  very  general  fading  channel  model  developed  above 
is  difficult  to  use  in  the  performance  analysis  of  communication 
systems  due  to  the  mathematical  complexities  involved.  Furthermore, 
as  has  been  shown  by  Bello  (1963),  a  simpler  model  is  warranted  for 
most  radio  channels.  In  the  next  section,  we  will  develop  such  a 
model . 

2.3.  The  WSSUS  Fading  Channel 

We  now  develop  the  wide-sense  stationary  uncorrelated  scattering 
(WSSUS)  fading  channel  model  from  the  general  fading  channel  model  (2.12) 
presented  in  Section  (2.2).  In  this  development,  we  will  also  define  the 
wide-sense  stationary  (WSS)  channel  and  the  uncorrelated  scattering  (US) 
channel. 

In  Section  2 . 1  we  noted  that  some  channels  exhibit  two  types  of 
fading:  long-term  fading  and  short-term  fading.  Short  term  fading  over 
these  channels  is  often  such  that  the  short  term  fading  statistics  are 
approximately  stationary  over  time.  Hence,  it  is  convenient  to  define 
a  subclass  of  the  general  fading  channel  model  known  as  wide-sense 
stationary  (WSS)  channels.  Bello  (1963)  defines  the  WSS  channel  as  a 
chattel  whose  correlation  functions  Rp(*)  and  A(*)  are  invariant  under 
a  translation  in  time.  The  space-time  cross-covariance  function  given 
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by  (2.22)  satisfies 

A(TrT2;trt2)  =  A(T1»T2itrt2)  (2-29) 

while  the  frequency- time  cross-covariance  function  given  by  (2.26)  satisfies 
RF(f1,f2;t1,t2)  =  RF(f1,f2;t1-t2)  (2.30) 

for  the  WSS  channel. 

While  developing  the  general  fading  channel  model  in  Section  2.2, 
we  assumed  that  the  scattering  medium  could  be  modeled  as  consisting 
of  differential  layers.  A  reasonable  assumption  for  many  channels  is  to 
assume  that  the  complex  Gaussian  process  f3(T,t)  is  independent  of  P(cr,t) 
for  t  ^  a.  Note  that  this  is  equivalent  to  assuming  that  the  effect  of 
scatterers  in  one  differential  layer  is  independent  of  the  effect  of  the 
scatterers  in  all  other  differential  layers.  The  space-time  cross¬ 
covariance  function  for  such  a  channel  is  given  by 

A(T1,T2;t1,t2)  =  A(Ti5ti’t2)6(Ti  '  T2)*  (2.31) 

Channels  whose  space- ^  me  cross -covariance  functions  satisfy  (2.31) 
are  known  as  uncorrelated  scattering  (US)  channels  [Bello,  1963].  We 
remark  in  passing  that  the  US  channel  is  the  wide -sense  dual  of  the  WSS 
channel  using  the  definitions  of  duality  given  by  Bello  (1964). 

Channels  which  exhibit  both  WSS  channel  characteristics  and  US 
channel  characteristics  are  known  as  wide-sense  stationary  >  correlated 
scattering  (WSSUS)  channels  [Bello,  1963].  Using  (2.29)  and  31),  the 
space-time  cross -covariance  function  for  a  WSSUS  channel  is  given  by 
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MTl’T2;tl,t2)  =  p(T1,t1-t2)6(T1-T2),  (2.32) 

where  for  convenience  we  have  let  p(T^,t^-t2)  =  A(T  ^  »T2  >  c]_»  c2)  f°r  the 
special  case  of  a  WSSUS  channel.  Note  that  for  WSSUS  channels,  two 
assumptions  on  A(«)  are  being  made:  (1)  the  scattering  processes  due  to 
different  layers  are  statistically  uncorrelated  and  (2)  the  scattering 
processes  in  each  layer  are  wide-sense  stationary.  Since  P(T,t)  is  a  zero- 
mean  complex  Gaussian  process  which  satisfies  (2.23),  wide-sense  stationarity 
of  3 (T , t)  implies  strict-sense  stationarity  of  3(T,t)  [Miller,  1974]. 

Although  the  above  channel  models  simplify  the  determination  of  the 
performance  of  conmuni cation  systems  considerably,  they  would  be  of 
little  value  unless  these  models  correspond  to  actual  fading  channels. 

Fortunately,  most  radio  channels  appear  to  exhibit  WSSUS  channel  properties 
([Bello,  1963],  [Stein,  1966]).  As  noted  in  Section  2.1,  radio  channel 
fading  is  often  characterized  by  the  superposition  of  short-term  fading 
on  long-term  fading.  The  short-term  fading  is  usually  found  to  exhibit 
stationary  statistics  while  the  long-term  fading  is  often  highly 
non- stationary,  depending  upon  the  time  interval  of  interest.  Bello  (1963) 
has  introduced  the  term  quasi-wide-sense  stationary  uncorrelated  scattering 
(QWSSUS)  to  describe  such  a  channel.  As  might  be  expected,  a  QWSSUS 
channel  has  WSSUS  channel  characteristics  over  time  intervals  on  the 
order  of  the  duration  of  short-term  fading.  Over  longer  time  intervals, 
the  channel  correlation  functions  no  longer  exhibit  stationarity. 

However,  Bello  (1963)  has  noted  that  the  performance  of  communication 
systems  over  QWSSUS  channels  may  be  evaluated  by  computing  system 


performance  by  modeling  the  channel  as  a  WSSUS  channel,  and  then  averaging 
the  short-time  system  performance  over  the  long-term  fading  statistics 
of  the  channel.  In  light  of  this,  for  the  remainder  of  this  chapter, 
we  will  consider  the  WSSUS  channel  model  only.  An  additional  reason 
for  considering  this  model  is  that  the  WSSUS  channel  is  the  simplest 
nondegenerate  channel  which  exhibits  both  time-  and  frequency-selective 
behavior . 

2.4.  Doubly-Spread  Channels 

In  the  previous  section,  we  developed  the  WSSUS  fading  channel 
model  from  the  general  fading  channel  model.  In  this  section,  we 
discuss  the  most  general  class  of  WSSUS  channels,  known  either  as 
doubly-spread  channels  [Van  Trees,  1971]  or  doubly-dispersive  channels 
[Kennedy,  1969].  Our  goal  is  to  characterize  the  various  parameters 
of  doubly-spread  channels  and  to  introduce  notation  often  used  in 
conjunction  with  discussions  of  fading  channels. 

A  word  about  terminology  is  first  in  order.  Doubly-spread  channels 
are  so  called  because  they  spread  the  time  and  frequency  waveforms  of 
a  signal  transmitted  through  the  channel.  Demonstration  of  this 
spreading  in  both  domains  must  wait  until  discussion  of  the  singly-spread 
degenerate  channels  which  exhibit  spreading  in  only  one  domain. 

For  the  doubly-spread  channel,  the  space-time  cross-covariance 

function  of  the  channel  is  given  by  (2.32).  The  scattering  function 

S  (T,f)  of  the  channel  is  defined  to  be 
DK 

00 

S  (T,f)  e"j2TTftp(T,t)dt,  (2.33) 

-00 

where  p(f,t)  is  defined  implicitly  by  (2.32).  Note  that  the  scattering 
function  is  the  temporal  Fourier  transform  of  p  (t , t)  . 


The  scattering  function  provides  a  means  of  characterizing  the 
Doppler  spread  of  the  channel,  as  will  be  seen  shortly.  It  should  be 
noted  that  most  of  the  results  to  this  point  and  for  the  rest  of  this 
chapter  could  just  as  well  have  been  developed  in  terms  of  the  channel 
scattering  function  instead  of  the  channel  covariance  function;  see 
Kennedy  (1969)  for  such  an  approach.  Because  the  channel  scattering 
function  and  the  channel  covariance  function  are  Fourier  transform  pairs 
such  parallelism  is  to  be  expected. 

If  the  scattering  function  is  concentrated  in  one  region  of  the 
t , f  plane,  it  may  be  characterized  grossly  in  terms  of  its  moments.  We 
make  the  following  definitions  [Van  Trees,  1971]: 

The  mean  delay  is  defined  to  be 

00 

■rii'/ T  sD8(T’f)dTdf-  <2- 

2c, 

b 

The  mean-square  delay  spread  is  defined  to  be 

A  X  p  f1  2  2 

L  =  “T  JJ  T  SDR(x,f)dTdf  -  .  (2. 

b 

The  mean  Doppler  shift  is  defined  to  be 

00 

"D  -  72  SDR(T-f)dTd£-  (2- 

2ab  - 

The  mean-square  Doppler  spread  is  defined  to  be 

B  "  -^JT  f2SDR(T,f)drdf  -  m£  .  (2. 
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In  (2 . 34) - (2 . 37) , 

00 

2ab  =  JT  SDR(T,f)dTdf  '  (2-38) 

■  CO 

In  addition,  we  will  define  the  duration  T  and  the  bandwidth  W  of  a 
narrowband  bandpass  signal  given  by  (2.1)  to  be 

03 

T  =  ^J  t2|u  (t)|2dt  (2.39) 

t  -00 

and 

00 

f2|u0(t)|2df  (2.40) 

t 

where 

Ec  3  J  lu0(c)|2dt  (2.41) 

—CD 

and  Uq (t)  is  defined  by  (2.2). 

This  definition  of  the  mean-square  properties  L,  B,  T,  and  W  is 
not  unique;  Kennedy  (1969)  defines  these  properties  using  Lerner's  (1959) 
definition  of  the  duration  of  a  signal.  The  definitions  given  above  are 
given  only  for  completeness.  For  strictly  time-  or  band-limited  scattering 
functions  and/or  bandpass  signals,  simpler  definitions  will  be  used  for 
L ,  B ,  T ,  and  W . 

An  underspread  channel  is  defined  to  be  one  for  which  [Van  Trees, 

1971] 

BL  <  1;  (2.42) 

similarly,  an  overspread  channel  is  defined  to  be  one  for  which 
[Van  Trees,  1971] 


BL  >  1  . 


(2.43) 
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In  reviewing  the  literature  on  doubly-spread  channels,  the  terms 
correlation  (or  coherence)  time  and  correlation  (or  coherence)  bandwidth 
are  often  encountered  in  characterizing  the  channel  [Kennedy,  1969].  The 
correlation  time  of  a  fading  channel  is  defined  to  be  the  time  separa¬ 
tion  t  beyond  which  samples  of  the  received  complex  envelope  are 
independent.  Since  the  channel  scattering  process  is  assumed  to  be 
modeled  as  a  zero-mean  complex  Gaussian  process,  given  the  transmitted 
signal,  the  received  signal  envelope  is  also  a  zero-mean  complex 
Gaussian  process.  Thus  independence  of  time  samples  is  implied  if  the 
correlation  function  of  the  envelope  is  zero.  Using  (2.32)  and  (2.14), 
the  correlation  between  time  samples  of  the  received  signal  envelope  is 

00 

RTE(tl,t2)  “  *  E[u(t1)u*(t  )3  »  J  PI(T,t1-t2)u0(t1-T)u*(t2-T)dT  .  (2.44) 

.00 

By  convention,  we  choose  the  correlation  time  for  the  channel  to  be  the 
smallest  time  separation  =  t^-t2  for  which 

RTE(t1,t2)  =  0.  (2.45) 

In  a  similar  fashion,  the  correlation  bandwidth  of  a  fading  channel 
is  defined  to  be  the  frequency  separation  beyond  which  samples  of  the 
Fourier  transform  of  the  received  complex  envelope  are  independent. 

From  the  comments  above,  the  received  signal  envelope  is  a  complex 
Gaussian  process  and  it  therefore  follows  that  the  Fourier  transform  of 
the  received  envelope  is  also  a  complex  Gaussian  process.  Thus 
independence  of  frequency  samples  is  implied  if  the  correlation  function 
of  the  Fourier  transform  of  the  received  envelope  is  zero.  From  (2.14), 
the  Fourier  transform  of  the  received  envelope  is 
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u(f)  =  J’j*  P  (T  ,  t)uQ(t-T)exp(-j2  TTft)dTdt.  (2.46) 

-00 

Using  (2.32),  the  correlation  between  frequency  samples  of  the  Fourier 
transform  of  the  received  signal  envelope  is  given  by 

RFE(fl’f2}  "  *  E[U(f1)U*(f2)} 

00 

=  JTJP  ^  » t1-t2)u0(t1-T)u*(t2-T)exp[j2TT(f2t2-f1t1)  JdTdt^^ 

—  00 

(2.47) 

By  convention,  we  choose  the  correlation  bandwidth  for  the  channel  to  be 
the  smallest  frequency  separation  =  f^-f  for  which 


^FE  1’  ^2^  *  °  ’ 


(2.48) 


In  Sections  2.5  and  2.6  simpler  expressions  are  derived  for  the 
correlation  functions  of  the  received  envelope  for  time  separations  and 
frequency  separations,  respectively. 

Finally,  we  note  that  in  the  above,  the  doubly-spread  channel 
was  characterized  in  terms  of  its  channel  covariance  function  and  the 
temporal  Fourier  transform  of  this  quantity,  the  channel  scattering 
function.  Alternatively,  we  could  have  characterized  the  channel  by  the 
spatial  Fourier  transforms  of  these  two  quantities: 


and 


RDR(v,t)  *  J  P  (T  ,t)exp(-j2  TTVT)dT 
—00 

00 

PDR^V,f)  “  i*  SDR^T  >f)exP<-J2  "vTjdT 

—00 


(2.49) 


50) 


(2. 
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The  quantity  R^(f,v)  defined  in  (2.49)  is  known  as  the  two-frequency 
correlation  function  [Van  Trees,  1971];  the  quantity  P  (f,v)  in  (2.50) 
is  defined  to  be  the  Doppler  cross-power  spectral  density  [Bello,  1963]. 

The  usefulness  of  having  the  four  quantities  p  (t , t) ,  SD^(T,f),  (v,t), 

and  P  (v,f)  comes  in  characterizing  doubly-spread  channels  which  have 
correlation  functions  which  are  concentrated  in  one  or  more  of  the 
variables  time  t,  Doppler  spread  f,  delay  spread  t,  or  delay  frequency  v. 

We  shall  find  in  Section  2.6  that  the  two-frequency  correlation  function 
is  useful  in  parameterizing  a  subclass  of  doubly-spread  channels. 

In  this  section  we  have  defined  and  characterized  the  most  general 
WSSUS  fading  channel,  the  doubly-spread  channel.  Often  fading  radio 
channels  exhibit  spreading  predominately  in  either  the  time  or  frequency 
domains  only.  In  still  other  cases,  the  fading  effects  are  such  that 
they  may  be  modeled  by  a  random  variab le  instead  of  a  random  process.  In 
such  cases,  it  is  convenient  to  define  subclasses  of  doubly-spread  channels 
having  specific  characteristics.  In  the  next  three  sections  we  will 
develop  the  models  for  three  subclasses  of  doubly-spread  fading  channels, 
also  known  as  degenerate  channels  [Van  Trees,  1971]. 

2.5.  Time -Selective  Fading  Channels 

In  Section  2.2,  we  developed  the  general  model  of  a  fading  channel  by 
assuming  that  the  scattering  medium  could  be  modeled  as  a  randomly  moving , 
layered  volume  of  scatterers;  each  layer  of  which  could  be  modeled  as  a 
complex  Gaussian  process.  In  this  section  we  will  assume  that  the 
scattering  medium  can  be  modeled  as  a  single  layer  of  randomly  moving 
scatterers.  As  will  be  seen,  such  an  assumption  leads  to  the  development 
of  a  class  of  channels  known  as  time -selective  fading  channels. 
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To  begin  with,  assume  that  a  signal  SqCO  given  by  (2.1)  is 
transmitted  through  a  scattering  medium  that  can  be  modeled  as  a  single 


layer.  The  received  signal  is 


s(t)  =  Re[P(t-t0)uQ(t-t0)exp(j2  tt  f0(t-tQ))  ] 


(2.51) 


where  t^  is  the  propagation  time  between  the  transmitter  and  receiver  and  P(t) 
is  the  time-varying  transmission  coefficient  due  to  the  scattering  medium. 

From  the  comments  in  Section  2.1,  we  will  assume  8 (t)  is  a  sample  function 
from  a  zero-mean  complex  Gaussian  random  process.  Equation  (2.51)  may 
also  be  written  as 


where 


s(t)  =  Re[u(t-tQ)exp(j2  tt  fQ(t-t0))  ] 


u(t)  »  P(t)uQ(t) 


(2.52) 


(2.53) 


For  the  WSSUS  channel  model,  j3(t)  is  a  stationary  random  process. 

It  is  important  to  note  that  since  8(t)  is  a  complex-valued  process  it 
influences  both  the  amplitude  and  the  phase  of  the  transmitted  signal  Sg(t) 

Equations  (2.51)  or  (2.52)  and  (2.53)  together  with  the  condition 
that  P(t)  be  a  stationary  process  constitute  what  is  known  as  the  time- 
selective  fadir\  channel  model  [Bello  and  Nelin,  1963].  Other  adjectives 
often  used  to  describe  this  channel  model  include  channels  dispersive  only 
in  frequency  [Kennedy,  1969],  frequency-flat  fading  channels  [Bello  and 
Nelin,  1963],  or  Doppler-spread  channels  [Van  Trees,  1971].  The  basis  of 
this  terminology  will  become  clear  as  the  properties  of  time-selective 
fading  channels  are  discussed. 


4 
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Equation  (2.51)  could  just  as  easily  have  been  derived  directly  from 
the  doubly-spread  fading  channel  model  by  noting  that  (2.51)  is  identical 
in  form  to  (2.11),  except  for  the  argument  of  (3.  Hence  the  general 
fading  channel  model  (2.12)  reduces  to  (2.11)  where  the  range  variable 
t  of  (3  is  no  longer  important.  In  terms  of  covariance  functions  of  the 
channel,  the  channel  covariance  function  for  the  time-selective  fading 
channel  is 

p(T,t-s)  =  p  (0,t-s)6  (t)  .  (2.54) 

In  the  following,  we  will  use  (2.54)  for  the  channel  correlation  function 
for  time-selective  fading  channels  rather  than  define  a  new  time  auto¬ 
correlation  function  that  directly  characterizes  (3(t)  in  (2.51). 

Some  comments  on  the  form  of  (2.53)  are  in  order.  First,  since  Ug(t), 
the  transmitted  envelope,  is  being  multiplied  by  P(t),  which  is  independent 
of  frequency,  the  received  envelope  u(t)  in  (2.53)  is  easily  seen  to 
undergo  fading  that  is  independent  of  frequency:  i.e.,  the  various 
frequency  components  of  Ug(t)  fade  identically  (frequency-flat  fading). 
Nevertheless,  (3(t)  in  (2.53)  is  a  time-varying  function  and  as  such, 
acts  to  modulate  the  transmitted  envelope  Ug(t).  This  leads  to  spreading 
of  the  Fourier  transform  of  the  transmitted  envelope  in  the  frequency- 
domain;  hence,  the  origin  of  the  term  Doppler-spread  fading.  Finally, 
since  P(t)  in  (2.53)  is  independent  of  t,  the  uncorrelated  scattering 
condition  of  WSSUS  channels  is  not  really  necessary  in  characterizing 
time- selective  fading  channels. 
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Using  (2.54)  in  (2.33),  the  scattering  function  for  a  time- selective 
fading  channel  is 

00 

SDR(r,f)  -  SD(0,f)5(r)  -  6(t)J  e"j2TTftp(0,t)dt.  (2.55) 

.CO 

The  term  SD(0,f)  in  (2.55)  is  known  as  the  Doppler  scattering  function 
of  the  channel. 

From  (2.55)  and  (2.36),  we  see  that  the  mean  Doppler  shift  for  a 
time- selective  fading  channel  is 

00 

nu  =  J  fS  (0,f)df;  (2.56) 

2  a  -« 

similarly,  from  (2.55)  and  (2.37)  we  see  that  the  mean-square  Doppler  spread  is 

00 

B  =  J  f2S  (0,f)df  -  m?  (2.57) 

2a  -® 

D 

where 

2(Jb  =  J  sD(°»f>df-  (2-58> 

.09 

The  mean-square  Doppler  spread  characterizes  the  Doppler  spread  around  the 
mean  Doppler  shift  due  to  the  time- selective  properties  of  the  fading 
channe 1 . 

Using  (2.55)  in  (2.34)  and  (2.35),  it  is  readily  seen  that  the  mean 
delay  and  mean-square  delay  spread  for  a  time -selective  fading  channel 
are  identically  zero: 


L  =  0. 


(2.59) 


Thus  the  time -selective  fading  channel  exhibits  spreading  in  frequency 
but  not  in  delay  and  is  therefore  often  called  a  singly-spread  channel. 
Another  singly-spread  channel  will  be  encountered  in  the  next  section. 

To  clarify  the  relationship  between  correlation  time,  the  channel 
correlation  function,  and  the  correlation  between  time  samples  of  the 
received  signal  envelope  for  a  time-selective  fading  channel,  consider  the 
following  example.  Note  first  that  by  using  (2.54)  in  (2.44),  the 
correlation  between  time  samples  of  the  received  signal  envelope  for  a 
time-selective  fading  channel  is 


RTED^tl5  *2^  *  P (°» ci“t2^u0^tl^u0^t2^ 


(2.6 


For  this  example,  let  the  Doppler  scattering  function  for  the  channel  be 
of  the  form 


sD(o,f)  = 


1‘  -  —  <  f  <  — 

2*2 


(2.6 


0;  elsewhere  . 


The  corresponding  channel  correlation  function  p(0,At)  is  therefore 


P(0, At)  »  B 


sirrnBAt 


rrBAt 


The  correlation  time  as  defined  by  (2.45)  is  given  by 


c  B  * 


(2.6 


provided 


uQ(t1)u^(t1  -  At)  4  0  for  0  <  At  <  - 


(2.6! 


Thus  time  samples  of  the  received  signal  separated  by  1/B  seconds  will 


be  uncorrelated. 
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Suppose,  however,  that  we  transmit  a  pulse  of  duration  T  through 
the  channel  where 

T  «  £  .  (2.64) 

In  this  case,  it  is  the  pulse  itself  that  determines  and  the 
correlation  time  is  equal  to  the  duration  of  the  pulse.  From  (2.60) 
and  this  example,  we  may  conclude  that,  provided  (2.64)  holds,  time 
samples  of  the  received  envelope  separated  fay  less  than  the  pulse 
duration  T  will  be  correlated.  In  Section  2.7  we  will  see  that  (2.64) 
is  a  necessary  condition  for  the  development  of  the  non-dispersive  fading 
channel  model. 

2.6.  Frequency-Selective  Fading  Channels 

The  second  subclass  of  doubly-spread  fading  channels  that  we  will 
look  at  are  known  as  frequency-selective  fading  channels  [Bello  and  Nelin, 
1963].  As  is  the  case  for  time-selective  fading  channels,  frequency- 
selective  fading  channels  result  from  doubly-spread  channels  by  making 
different  assumptions  about  characteristics  of  the  scattering  medium.  For 
frequency-selective  fading  channels,  we  assume  that  the  scattering  medium 
may  be  modeled  as  a  fixed  (non-moving)  volume  consisting  of  differential 
layers.  For  a  transmitted  signal  Sg(t)  given  by  (2.1),  the  received 


signal  scattered  from  a  single  layer  is 
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s ( t ,t )  =  Re[P  (T)u0(t-tQ-T)exp(j2  TT  f0(t-tQ-T))  ]dT  (2.65) 


where  is  the  nominal  propagation  delay  between  the  transmitter  and 
receiver,  r  is  the  additional  delay  due  to  the  scattering  layer  and  P(t) 
is  the  transmission  coefficient  for  a  wave  scattering  off  a  layer  whose 
propagation  delay  is  t.  Note  in  particular  that  P(t)  is  a  random  variable 
that  is  indexed  by  the  additional  delay  t.  We  will  assume  P(t)  is  a 
zero-mean  complex  random  variable. 

The  total  received  signal  is  a  superposition  of  all  the  responses  due 
to  the  individual  layers: 


s(t)  =  Re 


J 


■  i 

P  (T)u0(t-t0-r)exp(j2  rr  f0(t“t0'T))dT! 


(2.66) 


where  (3(t)  is  assumed  to  be  zero  except  for  the  finite  region  where 
scatters  exist.  Since  P(t)  is  a  zero-mean  complex  Gaussian  random 
variable  with  a  Rayleigh  distributed  amplitude  and  a  uniformly  distributed 
phase,  the  factor  exp(- j2tTf  )  may  be  absorbed  into  P  (t )  .  In  this  case 
(2.66)  becomes 


s(t)  =  Re 


)  P(T)U  (t-t  - 

•CD 


r )exp(j2  tt  fQ(t-t0))dT| 


or  equivalently 

s(t)  =  Re[u(t-tQ)exp(j2  tt  fQ(t-t0))  ] 


(2.67) 


(2.68) 


where 


U  (t)  =J  P(t)u  (t-T)dT  *J  P(t-T)u  (T)dT  . 


(2.69) 
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Equations  (2.67)  or  (2.68)  and  (2.69)  together  with  the  condition  that  the 
individual  scattering  layers  are  uncorrelated  constitute  the  frequency- 
selective  fading  channel  model.  Observe  that  we  could  have  derived  (2.67) 
directly  from  (2.12)  by  simply  dropping  the  time  dependence  of  (3  (t  ,  t)  . 

Thus  frequency  selective  fading  channels  are  properly  designated  a  subclass 
of  doubly-spread  channels.  Note  that  (2.69)  describes  the  input-output 
characteristics  for  a  linear  time -invariant  filter  with  complex  impulse 
response  P(t).  For  reasons  that  will  become  apparent,  frequency-selective 
fading  channels  are  also  known  as  time-flat  fading  channels  [Bello  and 
Nelin,  19631,  channels  dispersive  only  in  time  [Kennedy,  1969],  and  range 
or  delay-spread  fading  channels  [Van  Trees, 1971]. 

The  channel  covariance  function  for  a  frequency-selective  fading 
channel  is 

g(T)  =  P(T,0)  .  (2.70) 

The  term  g(r)  is  also  known  as  the  range  scattering  function  for  a 
frequency-selective  fading  channel. 

The  channel  scattering  function  for  a  frequency-selective  fading 
channel  may  be  found  by  using  (2.70)  in  (2.33): 

SDR(r,f)  =  g(r)6(f)  .  (2.71) 

Using  (2.71),  the  mean  delay  for  a  frequency-selective  fading  channel  is 

00 

iru  =  J  t  g  (t  )  dT  (2.72) 

la  -« 

D 

while  the  mean-square  delay  spread  is  given  by 
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ti  r 


t  1  p  2  ,  . .  2 

L  =  — 2  J  T  g(T)dT  - 


(2.73) 


where 


2ct  *  «  J  g  (t  )  dT  . 


(2.74) 


Using  (2.71)  in  (2.36)  and  (2.37),  the  mean  Doppler  shift  and  the 
mean-square  Doppler  spread  for  a  frequency-selective  fading  channel  are 
found  to  be  identically  zero. 

We  will  now  demonstrate,  by  way  of  example,  the  relationship  between  the 
two-frequency  correlation  function,  the  transmitted  signal  envelope,  and  the 
correlation  between  frequency  samples  of  the  received  signal  envelope.  To 
do  this  we  first  need  to  compute  the  correlation  between  frequency  samples 
of  the  Fourier  transform  of  the  received  signal  envelope  when  a  bandlimited 
signal  is  transmitted.  Using  (2.70)  in  (2.47),  the  frequency  correlation 
function  RFE(fpf£)  becomes 


^FE (f  1* f2^  “  JJJg<T)uo(ti‘T^uo^t2'T^exp^2  n  Cf2t2"fit:i^dTdtidt2  *  (2'75> 


Equation  (2.75)  may  be  rewritten  as 


RFE(f1,f2)  =  UQ(f1)U*(f2)  J  g(T)exp[-j2  ttt  (f1-f2)  ]dT 


(2.76) 


where  Ug(f)  is  the  Fourier  transform  of  u^(t)  as  defined  in  (2.2).  Note 


that  the  quantity 


R^f.O)  =  J  g(T)exp(-j2  HTf)dT 


(2.77) 


appearing  in  (2.76)  is  the  two- frequency  correlation  function  R^E(f,t) 
for  zero  time  separation  as  defined  by  (2.49).  Using  (2.77),  (2.76)  may 


be  written  as 
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RFE(fl,f2^  ~  U0('fl)U0('f2^RR<'fl"f2,0) 


(2.78) 


Now  let  g (t  )  be  of  the  form 


/  1  L  ^  .  L 

1  -  —  <  T  <  — 

2  “  2 


g(T)  = 


(2.79) 


0,  elsewhere  • 

Using  (2.77),  the  two-frequency  correlation  function  for  the  channel 
covariance  function  given  by  (2.79)  is 

sin  TTL(f1-  f2) 


^R(f  l”f2  ’  °>  L  tt  L(f1-f2) 


(2.80) 


Assume  a  signal  is  transmitted  whose  envelope  funccion  has  a  Fourier 
transform  given  by 


Vf)  = 


fi  -  S£<  f  <s 

2  “  2 


0,  elsewhere 


(2.81) 


Using  (2.80)  and  (2.81)  in  (2.78),  it  is  readily  seen  that  if 

„>i 


(2.82) 


then  the  frequency  components  of  Ug(f)  separated  by  multiples  of  1/L  Hz 
will  be  uncorrelated.  From  (2.48),  we  see  that  the  correlation  bandwidth 
of  the  channel  for  this  case  is 


W 


c  L 


(2.83) 


If,  however,  instead  of  condition  (2.82)  we  have  the  condition 

W  «  i,  (2.84) 

then  using  (2.80)  and  (2.81)  in  (2.78),  we  see  that  samples  of  the 
Fourier  transform  of  the  envelope  of  the  received  signal  are  correlated 
for  all  frequencies  in  the  bandwidth  of  the  signal.  In  Section  2.7  we 
will  see  that  (2.84)  is  a  necessary  condition  for  the  development  of  the 
non-dispersive  fading  channel  model. 

Before  leaving  this  section,  we  make  mention  of  a  concept  often 
encountered  in  the  literature  on  fading  channels,  which  is  that  of 
time- frequency  duality  ([Bello,  1964],  [Van  Trees,  Sec.  12.3,  1971]). 

The  reader  may  have  noticed  a  certain  similarity  between  the  equations 
characterizing  a  time-selective  fading  channel  and  the  equations 
characterizing  a  frequency-selective  fading  channel.  As  an  example 
of  this,  compare  the  expressions  for  the  correlation  between  time 
samples  of  the  received  signal  envelope  for  a  time- selective  fading 
channel  given  by  (2.60)  and  the  correlation  between  frequency  samples 
of  the  Fourier  transform  of  the  received  signal  envelope  for  a 
frequency- selective  fading  channel  given  by  (2,78).  Except  for  the  fact 
that  time  appears  in  the  arguments  of  (2.60)  while  frequency  appears  in 
the  arguments  of  (2.78)  and  that  Fourier  transforms  of  the  terms  of  (2.60) 


are  used  in  (2.78),  the  two  expressions  are  identical  in  form.  To 
emphasize  this  "duality"  concept,  (2.69)  may  be  written  as 

U (f )  =  H(f)U0(f)  (2.85) 

where  H(f)  is  the  Fourier  transform  of  a  sample  function  of  P(t): 

CO 

H(f)  =  J  (3  (t  )exp(-j2  n  f  t  )dT  .  (2.86) 

.OP 

Comparing  (2.53)  with  (2.85),  we  see  that  except  for  the  fact  that 
Fourier  transforms  are  used  in  (2.85),  the  defining  equations  for  the 
two  channel  models  are  identical  in  form.  As  was  mentioned  in  Section 
2.3,  even  the  conditions  of  wide-sense  stationarity  for  the  time-selective 
fading  channel  model  and  uncorrelated  scattering  for  the  frequency- 
selective  fading  channel  model  are  dual  concepts.  Bello  (1969)  has 
thoroughly  discussed  the  various  concepts  associated  with  time- 
frequency  duality;  the  reader  is  urged  to  consult  [Bello,  1969]  for  more 
details.  Using  Bello's  definitions,  Van  Trees  (1971)  proves  that  a 
time-selective  fading  channel  is  the  dual  of  a  frequency- selective 
fading  channel.  The  principal  advantages  of  recognizing  this  duality 
appear  to  be  in  deriving  equivalent  circuit  models  of  fading  channels 
[Bello,  1963],  designing  optimal  receivers  for  the  two  types  of 
channels  [Van  Trees,  1971],  and  in  evaluating  the  performance  of 
communication  systems  over  time- selective  and  frequency-selective 
fading  channels  ([Bello  and  Nelin,  1962],  [Bello  and  Nelin,  1963]). 


In  the  previous  two  sections,  we  have  discussed  two  subclasses  of 
the  doubly-spread  fading  channel  model.  To  derive  the  channel  models 
for  these  subclasses  we  have  assumed  either  that  the  scattering  medium 


is  moving  but  may  be  modeled  as  a  single  layer  or  that  the  scattering 
medium  is  fixed  but  consists  of  differential  layers.  In  this  section 
we  assume  that  the  scattering  medium  is  both  fixed  (non-moving)  and  may 
be  modeled  as  a  single  layer. 

For  a  transmitted  narrowband  bandpass  signal  given  by  (2.1),  the 
received  signal,  after  passing  through  a  fixed,  single  layer  scattering 
medium,  is  given  by 

s(t)  =  Re[Pu0(t-t0)exp(j2nfQ(t-t0))]  (2.87) 

where  t^  is  the  propagation  time  between  the  transmitter  and  receiver 
and  (3  is  assumed  to  be  a  zero-mean  complex  Gaussian  random  variable. 
Equation  (2.87)  constitutes  the  fading  channel  model  for  a  nondispersive 
fading  channel  [Kennedy,  1969].  Note  that  (2.87)  could  have  been  derived 
directly  from  the  channel  model  for  a  doubly-spread  fading  channel  (2.12) 
by  dropping  the  time  dependency  of  P(T,t)  and  by  noting  that  scattering 
is  due  to  a  single  layer.  A  more  precise  derivation  of  the  nondispersive 
fading  channel  model  from  the  doubly-spread  fading  channel  model  is  given 


below. 


The  channel  covariance  function  for  a  nondispersive  fading  channel 
is  given  by 

p(T,At)  =  p  (0,0)6  (t)  , -®  <  t  <  ®,-®  <  At  <  ®.  (2.8 

Using  (2.88)  in  (2.33),  the  scattering  function  for  the  nondispersive 
fading  channel  is  given  by 

(T,f)  =  p  (0,0)6  (r)6(f).  (2.8 

From  (2.89)  and  (2 .34) - (2 .37) ,  the  mean  delay,  mean-square  delay  spread, 
mean  Doppler  shift,  and  mean-square  Doppler  shift  for  a  nondispersive 
fading  channel  are  all  identically  zero.  The  nondispersive  fading 
channel  thus  exhibits  no  spreading  in  frequency  or  time;  hence  the  name 
nondispersive.  Another  term  often  used  to  describe  the  nondispersive 
fading  channel  is  flat-flat  fading  channel  [Bello  and  Nelin,  1963]. 

Since  P  in  (2.87)  has  a  Rayleigh  distributed  amplitude  and  a  uniformly 
distributed  phase,  the  nondispersive  fading  channel  is  also  known  as 
simply  a  Rayleigh  fading  channel.  In  this  study,  however,  by  a 
Rayleigh  fading  channel  we  will  mean  any  doubly-spread  fading  channel 
or  any  of  its  subclasses  which  do  not  contain  specular  components. 

As  mentioned  above,  a  nondispersive  fading  channel  is  a  special 
case  of  a  doubly-spread  channel.  More  generally,  it  is  a  special 
case  of  both  time -selective  fading  and  frequency-selective  fading 
channels  (which  are,  in  turn,  special  cases  of  the  doubly-spread 
channel).  In  Section  2.5  we  saw  that,  for  a  time-selective  fading 
channel  and  a  signal  satisfying  (2.64),  time  samples  of  the  received 
envelope  would  be  correlated  for  the  duration  of  the  transmitted  envelope 


In  Section  2.6  we  saw  that,  for  a  frequency-selective  fading  channel  and 
a  signal  satisfying  (2.84),  frequency  samples  of  the  received  envelope 
would  be  correlated  for  all  frequencies  in  the  bandwidth  of  the  trans¬ 
mitted  signal.  Since  a  nondispersive  fading  channel  exhibits  no  spreading 
in  either  time  or  frequency,  (2.64)  and  (2.84)  combined  must  be  the 
condition  for  a  doubly-spread  channel  to  exhibit  nondispersive  fading: 

BL  «  ^  .  (2.90) 

Since  the  time -bandwidth  product  for  any  signal  must  be  greater  than 
unity  (c.f.  [Papoulis,  Sec.  4-4,  4-5,  1962]) 

TW  >  1,  (2.91) 

(2.91)  combined  with  (2.90)  indicates  that  nondispersive  fading  will 
occur  provided 

BL  «  1.  (2.92) 

From  (2.42),  it  is  seen  that  nondispersive  fading  will  occur  in  a  doubly- 
spread  channel  provided  that  the  channel  is  underspread. 

Finally  we  mention  a  special  case  of  the  nondispersive  fading 
channel  known  either  as  the  "slow  and  flat"  fading  channel  model 
[Nesenbergs,  1967]  or  the  pure  fading  channel  model  [Wozencraft  and  Jacobs, 
1965].  For  this  case,  8  appearing  in  (2.87)  is  assumed  to  be  a  real 
Rayleigh  (or  Rician,  if  a  specular  component  is  present)  distributed 
random  variable.  This  channel  model  will  result  from  the  nondispersive 
fading  channel  model  if  the  phase  of  P  may  be  modeled  as  a  real, 
nonrandom  variable. 

Figure  3  summarizes  the  fading  channel  models  that  we  have 
considered  up  to  this  point. 


GENERAL 

TIME-VARYING 

CHANNEL 


Summary  of  fading  channel  models.  Numbers  in  parentheses 
equations  for  the  channel. 
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CHAPTER  3 

AVERAGE  SNR  ANALYSIS  OF  DS/SSMA  COMMUNICATIONS  VIA  FADING  CHANNELS 
One  useful  measure  of  performance  of  a  DS/SSMA  system  is  the  average 
signal-to-noise  ratio  SNR.  In  this  chapter,  we  analyze  the  average 
signal-to-noise  ratio  performance  of  DS/SSMA  communications  for  a  general 
class  of  fading  channels.  In  making  this  analysis  we  shall  use  the  DS/SSMA 
model  from  Chapter  1  together  with  the  WSSUS  fading  channel  model  of 
Chapter  2.  The  results  presented  represent  a  generalization  of  the 
performance  analysis  of  [Pursley,  1977]  which  considered  only  additive 
white  Gaussian  noise  (AWGN)  channels. 

The  type  of  fading  considered  in  this  chapter  is  Rician  or  specular- 
plus-Rayleigh  fading  [Stein,  p.  372,  1966].  That  is,  for  a  single  trans¬ 
mitted  signal  given  by  (2.1),  the  received  signal  (2.18)  consists  of  a 
replica  of  the  transmitted  signal  plus  a  weaker  Rayleigh-faded  version  of 
this  signal.  This  will  be  the  situation  whenever  the  transmission  medium 
is  such  that  there  is  a  strong  stable  path  and  a  number  of  weak  paths.  The 
component  of  the  received  signal  that  is  produced  by  the  strong  stable  path 
is  called  the  specular  component  or  the  desired  signal  component;  this 
component  can  be  coherently  demodulated.  However,  unless  restrictions 
are  imposed  on  the  fading  rate,  the  Rayleigh-fading  signal  that  arises 
from  the  weak  paths  cannot  in  general  be  coherently  demodulated  [Stein, 
pp.  406-407,  1966],  and  hence  we  shall  treat  it  as  interference.  An 
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example  of  a  situation  in  which  this  type  of  fading  occurs  is  the  example 
given  in  Section  2.1  of  communication  between  an  aircraft  and  a  satellite. 
In  this  example  there  exists  a  direct  communication  path  between  the 
aircraft  and  the  satellite  in  addition  to  one  or  more  communication  paths 
due  to  reflection  off  the  earth. 

If  the  channel  exhibits  Rayleigh  fading  rather  than  Rician  fading 
(i.e.,  if  there  is  no  strong  stable  component  in  the  received  signal), 
then  the  direct-sequence  form  of  spread- spectrum  communication  would  not 
generally  be  suitable,  since  it  could  not  be  coherently  demodulated  unless 
the  fading  is  sufficiently  slow.  This  follows  from  the  results  of  Viterbi 
(1965)  who  showed  that  for  a  correlation  receiver  with  a  phase  reference 
having  a  uniform  probability  density  function,  as  would  be  the  case  for 
a  nondispersive  Rayleigh- fading  channel,  the  probability  of  error  for  an 
antipodal  signaling  set  is  %.  For  this  chapter,  it  is  assumed  that  the 
specular  component  is  of  sufficient  amplitude  that  the  effects  of  fading 
on  the  performance  of  the  DS/SSMA  synchronization  subsystem  may  be 
neglected.  The  performance  of  synchronization  subsystems  in  the  presence 
of  fading  communication  channels  has  been  analyzed  elsewhere  (e.g.,  see 
[Weber,  1976]). 

3.1.  DS/SSMA  System  Performance  for  Doubly-Spread  Fading  Channels 

In  this  section,  we  shall  develop  the  results  for  the  average 
signal- to-noise  ratio  at  the  output  of  a  correlation  receiver  for 
DS/SSMA  communications  via  a  doubly-spread  channel.  In  Sections  3.2 
and  3.3  we  will  specialize  the  results  obtained  to  the  singly-spread 
time-selective  and  frequency- selective  channels,  respectively. 
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The  system  model  that  we  will  consider  is  shown  in  Fig.  4  for  K 

users.  The  system  model  as  shown  is  assumed  to  be  time-asynchronous, 

i.e.  there  exists  no  common  timing  reference  among  each  of  the  K-users. 

In  Fig.  4,  the  fading  channel  is  assumed  to  be  able  to  be  modeled  by  a 

linear  time-varying  filter  whose  transfer  function  is,  in  general, 

nondeterministic .  With  this  one  exception,  the  system  model  of  Fig.  4  is 

identical  to  the  system  model  of  Fig.  1  for  the  AWGN  channel. 

For  each  k  (1  <  k  ^  K)  the  k-th  user's  data  signal  b^(t)  is  a 

sequence  of  statistically  independent,  unit  amplitude,  positive  and 

negative,  rectangular  pulses  of  duration  T.  The  data  signal  for  the  k-th 

user  is  therefore  given  by  (1.1).  As  was  the  case  for  the  AWGN  channel 

model,  each  user  is  assigned  a  code  waveform  a^t)  which  consists  of  a 

periodic  sequence  of  unit  amplitude,  positive  and  negative,  rectangular 

pulses  of  duration  T  .  The  code  waveform  for  the  k-th  user  is  therefore 

c 

given  by  (1.2).  We  assume  that  each  signature  sequence  has  period 
N  =  T/Tc  so  that  there  is  one  code  period  per  data  symbol. 

The  data  signal  h^(t)  is  modulated  onto  the  phase-coded  carrier  c^(t), 
which  is  given  in  complex  envelope  form  by 

ck(t)  *  ,/2P  Re[ak(t)exp(jouct  +  j©^)  ] ,  (3.1) 

so  that  the  transmitted  signal  for  the  k-th  user  is 

sk(t)  =  Jl?  Re[ak(t)bk(t)exp(jcuct  +  j©k)  ]  .  (3.2) 

In  the  above  expressions,  P  is  the  common  signal  power,  9k  is  the  phase 
of  the  k-th  carrier,  and  is  the  common  carrier  frequency. 
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Figure  4.  DS/SSMA  system  model  with  a  fading  channel. 
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In  the  absence  of  fading,  the  received  signal  at  the  input  to  a 
receiver  is  given  by 

K 

r(t)  =  n(t)  +  Z  s,  (t-r  ),  (3.3) 

k=l  *  K 

where  n(t)  is  the  AWGN  term  and  accounts  for  the  nominal  propagation 
time  for  the  k-th  signal  as  well  as  for  the  time  asynchronism  between  the 
k-th  transmitter  and  the  other  K-l  transmitters  in  the  system. 

In  the  presence  of  fading,  however,  the  expression  for  the  received 
signal  at  the  input  to  a  receiver  becomes  more  complex.  From  Chapter  2, 
for  a  transmitted  signal  given  by  (2.1),  the  received  signal  at  the  output 
of  a  doubly-spread  fading  channel  is  given  by  (2.12)  or  (2.13)  and  (2.14). 
As  noted  in  the  introduction  to  this  chapter,  the  type  of  fading  being 
considered  is  specular-plus-Rayleigh  fading.  Thus  (2.18)  must  be  used  to 
describe  the  fading  process.  Finally,  in  addition  to  the  fading  present 
in  the  channel,  we  will  also  assume  that  additive  white  Gaussian  noise  is 
also  present.  Thus  if  x^(t)cos (<u  t+9k)  *-s  t*ie  Input  to  the  channel  from 
the  k-th  transmitter,  then  the  corresponding  output  yk(t)  Is  given  by 

yk(t)  *  ReCuk(t-Tk)exp[ juJc(t-Tk)  +  j©k] }  +  n(t),  (3.4) 

where  n(t)  is  the  AWGN  term  and 

00 

uk(t)  =  Yk  J*  &k(T>t)xk(t-T)dT  +  x  (t)  .  (3.5) 

•  0D 

The  first  term  of  (3.5)  is  the  faded  portion  of  the  received  signal  and 
the  second  term  is  the  specular  component.  The  non-negative  real 
parameter  Yk  is  the  transmission  coefficient  for  the  fading  channel  as 


seen  by  the  k-th  signal;  and  B^(r,t)  is  a  zero  mean  unit-energy  (i.e., 

2 

0^  =  h  in  (2.38))  complex  Gaussian  random  process  which  will  be  referred 
to  as  the  fading  process.  As  pointed  out  in  Section  2.2,  0k(r,t),  as  it 
appears  in  (3.5),  represents  the  equivalent  low-pass  time-varying  impulse 
response  for  the  faded  portion  of  the  fading  channel.  The  covariance 
function  for  the  fading  process  in  a  doubly-spread  channel  is  given  by 
(2.32): 

Ak(T,a;t,s)  =  ^E[Pk(T,t)8*(a,s)} 

=  pk(T,t-s)5  (t  -  o' )  .  (3.6) 

As  will  be  seen  below  (see  (3.21)),  for  practical  WSSUS  channels,  pk(T,t-s) 
is  a  real-valued  function. 

From  equations  (3. 2) -(3. 4),  it  follows  that  for  a  Rician  doubly-spread 
fading  channel  with  AWGN,  x^t)  *  \fl P  ak(t)bk(t),  so  that  the  received 
signal  at  the  input  to  a  receiver  is 

K 

r(t)  »  n(t)  +  I  Re[uk(t  -  Tk)exp{j  [u^t  +  cpk]}]  (3.7) 

k=l 

where  cpk  =  9^  -  tucT k .  If  the  received  signal  r(t)  is  the  input  to  a 
correlation  receiver  matched  to  s^(t),  the  corresponding  output  is 
T 

Z  -  f  r(t)Re{a  (t  -  t  )exp(ju>  t  +  jcp  )]dt.  (3.8) 

i  0  CL 

Since  we  are  interested  in  phase  angles  modulo  2tt  and  time  delays  modulo  T, 
we  shall  (without  loss  of  generality)  assume  =  0  and  9^  =  0  (and 
hence  cp^  =  0).  Therefore,  as  in  [Pursley,  1977],  the  phase  angles  cpk 
and  9k  and  the  time  delay  (for  k  ^  i)  are  all  measured  relative  to 
the  phase  and  time  delay  of  the  i-th  signal.  From  (3.7)  and  (3.8),  it 


follows  that,  except  for  double- frequency  terms, 


T  K  T  jcp 

Z.  =  f  n(t)a.  (t)cos  uu  t  dt  +  I  %  f  Re{u,  (t  -  t  )a.  (t)e  jdt,  (3.9 

1  ^  »  X  C  .  .  ^  rC  tC  1 

0  k=l  0 


where  we  have  used  the  identity  (2.10).  The  double  frequency  terms  which 
are  omitted  from  (3.9)  may  be  ignored  for  practical  implementations  of 
SSMA . 

As  in  [Pursley,  1977],  the  phase  angles,  time  delays,  and  data 

symbols  for  the  k-th  signal  (k  /  i)  are  modeled  as  mutually  independent 

random  variables  which  are  uniformly  distributed  on  [0,2rr],  [0,T],  and 

£-+-1,-1} ,  respectively.  The  average  signal-to-noise  ratio  SNR^  at  the 

output  of  the  i-th  correlation  receiver  is  then  defined  in  terms  of 

probabilistic  averages  (i.e.,  expected  values)  conditioned  upon 

b^(t)  for  0  <  t  <  T,  where,  without  loss  of  generality,  we  may  assume 

b. (t)  =  +1  over  this  interval, 
i 

For  convenience,  the  following  notation  will  be  introduced.  Let 


hk  i(T;t)  =  \(t  -  T>bk(t  -  T)aL(t) 


(3.1 


and 


fu  „  .  (T,a;t,s)  =  h,  .  (t  ;  t)h  .(a;s) 


k,n,i 


k,i 


n,  l 


(3.1 


for  1  <  k  <  K,  1  <  n  <  K,  1  <  i  <  K;  0  <  t  <  T,  0  <  s  <  T;  and  -®  <  t  ® 


a>  <  o  <  ®.  The  function  h.  .  is  denoted  by  h.  ,  and  f.  .  .  is  denoted  by 

i  i  '  i  '  ill  ' 


1,1,1 


f^.  In  terms  of  this  notation,  the  i-th  specular  component,  the  i-th 


(complex)  faded  signal  component,  the  AWGN  component,  and  the  complex 


interference  component  are  defined  by 


D  =  V?P  j  b^(t)dt 


T  ® 


F  =  J*  j  Y-Pi(T,t)h  (T  ;t)dTdt 
1  0  -® 


N.  =  j  n(t)a.(t)cos  a)  t  dt 

1  Jo  1 


K 


i .  -  E  (y,  i,  .  +  i;  .  > 

1  k=l  k  k,:L 


k^i 


where 


*k,i  S  JQ  J,  V1-'  -  Tk)hk,i(Tk  +  T'.t)exp(jcpk)dTdt 


^k.i  ’  I  hk,i(Tk;t)e*I’(%)dt 
0 


If  we  then  let  =  Re(p\]  and  1^  =  Re{]L},  we  see  from  (3.9)  that 


Z.  =  N.  +  D.  +  F.  +  I.  . 

i  x  i  i  i 


Since  EftL  =  EI^  =  EF^  =  0  and  ED^  =  *]% P  J  b^(t)dt  =  Ts/^P,  then 


EZ^  =  x/5p.  Assuming  that  n(t)  has  a  two-sided  spectral  density  %Nq, 


the  variance  of  the  noise  component  Z^  is 


Var  =  fcNQTf 


while  the  conditional  variance  of  the  i-th  specular  signal  component, 
conditioned  upon  b^(t)  =  +1  for  0  ^  t  <  T,  is  zero.  Since  the  fading 
process  has  zero  mean,  the  variance  of  the  faded  portion  of  i he  i-th 
signal  is  given  by 


Var  F.  =  E[F.]"  =  E[Re  F.]2  =  %E[Re[F2}  ]  +  J>E[F*F.],  (3.20 

which  follows  from  an  application  of  (2.10)  with  x  =  y  =  F  ,  and  from 

~  i  ~  i 2 

the  fact  F^F^  =  |F^|  is  real.  The  expressions  for  the  two  terms  in 
(3.20)  will  involve  the  covariance  function  A^(T,a;t,s)  given  by  (3.6)  and 
the  function  A^(r,<J;t,s)  given  by  (2.21).  In  particular 


XX®® 


e|f  |2  *  py2  f  r  J  j  a  (T,<J;t,s)f  (T ,CT;t,s)dTdCTdtds. 

1  A  A  ^ 


(3.21 


0  0  -®  -® 


The  expression  for  E[Re[f\}]  is  the  right-hand  side  of  (3.21)  with  A^ 
replaced  by  A^.  However  as  noted  in  Section  2.2,  we  will  assume  that 
the  WSSUS  fading  process  is  such  that  (2.23)  holds  so  that  eCf^  =  0 


and  thus 


Var  F.  =  %E|F.|2 


(3.22 


From  (3.21),  (3.22),  and  (3.6),  it  can  be  seen  that 


p  2  T  T  “ 

Var  F  =tV.  J  J  J  p  (T,t-s)f  (t  ,t  ;  t,  s)dTdtds  . 
1  z  x  0  0  -®  1  1 


(3.23 


By  a  similar  analysis  of  the  interference  component  1^  defined  by 
(3 . 15) -(3 . 17)  it  is  clear  that  since  the  fading  process  has  zero  mean 
and  the  data  signals  are  uncorrelated  and  have  zero  mean,  then 


Var  I  -  %E|I,|2  -  «■  I  (V?  v.  .  +  v-  .) 


(3.24 
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A  ^  2  A  ~  2 

where  =  E 1 ±  l  and  v£  ±  =  Ejl^J  are  given  by 


T  T  « 

EJ  I  I  Pk(T,t-s)fkjkji(T  +  Tk,T  +  Tk;t,s)dTdtdS 


(3.25) 


0  0  -eo 


vk,i  ’  E  ■fjJ’o  fk,k,i(Tk’Tk:t'5)dtds 


(3.26) 


In  general,  the  expression  for  v,  .  in  (3.25)  cannot  be  reduced  further 

K,  1 

without  additional  constraints  on  the  covariance  function  pk(r,t-s)  for 

the  fading  process.  The  expression  v'  .  in  (3.26),  however,  does  not 

k  ,  i 

depend  on  the  fading  process.  Since  the  assumptions  used  to  derive  v,  . 

K,1 

are  identical  to  those  used  to  derive  the  variance  of  the  k-th  interference 
component  for  the  AWGN  channel  of  Chapter  1,  (3.26)  can  be  simplified  as 
in  [Purs ley,  1977]  to  give 


'U - 3  rk  i> 

’  6N  ’ 


(3.27) 


where  r,  .  is  Che  interference  parameter  given  by  (1.16). 
k,  i 

The  signal-to-noise  ratio  SNR^  at  the  output  of  the  i-th  correlation 
receiver  is  defined  by 

A  -k 

SNR.  =  (ED , ) [ Var  F.  +  Var  I.  +  Var  N. ]  ^ 
i  i  i  1  i 


T  Jh?  [Var  Fi  +  Var  I  +  ^NQT ] .  (3.28) 


Notice  that  for  SNR^  defined  in  this  manner,  the  faded  portion  of  the 
i-th  signal  is  considered  to  be  an  interference  component  of  (the 
correlation  receiver  cannot  make  use  of  this  component) . 


Several  observations  about  the  results  obtained  up  to  this  point  will 
now  be  made.  First,  throughout  the  analysis  presented  above,  nowhere  has 
there  been  made  any  assumptions  about  the  independence  of  the  fading 
processes  3^(T,t)  for  1  ^  k  <  K,  apart  from  the  implicit  assumption  that 
P^(T,t)  is  independent  of  the  random  variables  b^(t),  cp^,  and  T  . 
Consequently,  the  above  results  are  valid  for  the  important  cases  where 
(1)  all  the  fading  processes  are  independent  and  (2)  all  the  fading  processes 
are  identical,  i.e.,  P^(r,t)  =  £  (t  ,  t)  for  all  k,  1  <  k  <  K.  Case  (1) 
would  arise,  for  example,  in  the  situation  where  multiple  mobile  users 
are  communicating  via  the  ionosphere,  using  DS/SSMA  communications.  For 
this  situation,  it  would  be  expected  that  each  of  the  transmitted  paths 
would  involve  different  sets  of  scatterers  and  hence,  following  an 
argument  similar  to  that  used  to  develop  the  uncorrelated  scattering 
channel  model  of  Chapter  2,  the  K  fading  processes  would  be  independent. 

Case  (2),  on  the  other  hand,  would  arise  if  a  satellite  is  transmitting 
multiple-user  data  to  a  single  aircraft  using  DS/SSMA  communications.  In 
this  case,  all  of  the  K  users  would  undergo  the  same  fading  resulting  in 
identical  fading  processes.  Thus  the  results  above  are  general  enough  to 
allow  both  of  these  fading  situations. 

Second,  from  the  results  above  (viz.,  (3.23),  (3.24),  (3.25),  and 
(3.27)),  we  see  that  SNR^  depends  upon  the  signature  sequences  through  the 
functions  f^  in  (3.23)  and  f^  ^  ^  in  (3.25)  and  through  the  interference 
parameters  r^  ^  in  (3.27).  In  addition,  SNR^  depends  on  the  channel 
covariance  function  pfc,  1  <  k  <  K,  via  (3.23)  and  (3.25). 
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Finally,  throughout  the  analysis  presented  above  we  have  used  a 
correlation  receiver  structure  for  the  receiver.  This  is  not  to  imply 
that  a  correlation  receiver  is  the  optimum  receiver  for  DS/SSMA  communica¬ 
tions  via  Rician  fading  channels.  In  fact,  even  for  DS/SSMA  communications 
via  AWGN  channels,  the  correlation  receiver  is  not  the  optimum  receiver. 

For  while  a  correlation  receiver  i£  optimum  for  an  antipodal  signalling 
set  and  an  AWGN  channel  (e.g.  see  [Viterbi,  1965] )>  a  rate  1  convolutional 
decoder  using  the  Viterbi  algorithm  is  optimal  for  DS/SSMA  communications 
via  an  AWGN  channel  [Schneider,  1979].  For  a  nondispersive  Rician  fading 
channel,  it  has  been  shown  ([Viterbi,  1965],  [Turin,  1958])  that  a  linear 
combination  of  the  optimal  coherent  and  optimal  noncoherent  receiver  is 
optimum  for  a  single  user  system.  Other  schemes  have  been  devised  to 
provide  low  bit  error  rates  in  the  presence  of  "slowly"  fading  channels 
[Monsen,  1973].  Thus,  it  is  evident  that  a  correlation  receiver  is  not 
optimum  for  DS/SSMA  communications  via  fading  channels  (however,  in 
general  the  optimal  receiver  is  not  known).  Nonetheless,  the  analysis 
presented  above  is  valuable  in  a  practical  sense,  since  existing 
DS/ spread- spectrum  receivers  do  use  simple  correlation  receiver  structures 
to  simplify  receiver  design  and  construction  [Cahn,  1973] .  Therefore  in 
the  analysis  that  follows  we  will  continue  to  use  correlation  receivers 
matched  to  the  i-th  (1  ^  i  <  K)  user's  code  waveform. 

3.2.  DS/SSMA  System  Performance  for  Time-Selective  Fading  Channels 
Analysis .  In  this  section,  consideration  is  given  to  one  subclass  of 
doubly-spread  channels  known  as  time-selective  fading  channels.  The 
channel  covariance  function  for  a  time-selective  fading  channel  is  given 
by  (2.54): 
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Pi(T,t-s)  =  pjL(0,t-s)6(T). 


(3.29) 


For  convenience,  p^(0,t-s)  is  denoted  by  r^(t-s).  It  follows  from  (3.23) 
and  (3.29)  that  the  variance  of  the  faded  portion  of  the  i-th  signal  is 


T  T 


Var  Fi  “  2"  Yi  i*  JQ  ri(t-s)bi(t)bi(s)dtds  . 


(3.30) 


The  expression  for  the  data  signal  is  given  by  (1.1).  Making  a  change 

of  variables,  u  =  t-s,  and  using  (1.1),  (3.30)  becomes 

p  2  T  T  0  u+T 

Var  F.  =  -r  V ft  P  f  r.  (u)dtdu  +  P  f  r.  (u)dtdu] 

1  ^  1  0  u  1  -T  O1 


38  I  Y.  [ J  r  (u)  (T-u)du  +  P  r  (u)  (T+u)du] 
0  -T 


“  P  Yi  i"  ri^  (T_u)du 


(3.31) 


In  a  similar  fashion,  the  variance  of  the  interference  component 
1^  may  be  found  by  using  (3.29)  in  (3.25).  Interchanging  the  order  of 
integration  and  expectation  in  (3.25),  we  find 


Vi  'LL  rk(t-s)ECtk,k,i(vvt-s)}dtds’ 


(3.32) 


where  the  expectation  is  over  b^(t)  and  T^.  Note  that  given  t^, 

is  a  semi-random  binary  process  [Papoulis,  1965]  with  an  autocorrelation 

function  for  t^,  t,  and  s  in  [0,T] 


j  1,  Tk  <  min(t,s) 
E[bu(t-Tk)bk(s-Tk)|Tk)  =  *  1,  T,  £  max (t, s) 


0,  otherwise. 


(3.33) 


L*X*J 


Evaluating  the  expectation  in  (3.32)  using  (3.33),  we  find  that  for  t  and  s 
in  [0,1] 


min(t,s) 


Et£k,k,i<Tk'Tk;l:>s)1  ■  Vt)at.<S>  T  [J„  ak(t-Tk)ak(s-Tk)dTk 


0 


+  i  .V'-W-V^kl 

max ( t , s ) 


1  p 


T+min(0,v) 


=  ai(t)ai(s)  -J  ak(u)ak(u-v)du 

max(0,v)  ^ 


(3.34) 


where  u  =  t-T^,  v  *  t-s  and  we  have  used  the  fact  that  a^(t)  is  periodic 
with  period  T.  For  v  positive  we  find 


E^fk,k,i^Tk,Tk;t,S^  =  T  ^c(t”s)ai(t)ai(s)  '>  t-s  >  °»  (3.35) 


while  for  v  negative  we  find 


E^fk,k,i(Tk,Tk;t,s)^  =  T  \(|t-s|)a..  (^a^s);  t-s  >  0,  (3.36) 


where  R^(T)  is  the  continuous -time  partial  autocorrelation  function 
defined  in  Chapter  1.  Combining  (3.35)  and  (3.36),  the  expectation  in 
(3.32)  reduces  to 


Et£k,k.i(vTk;t's)J  -Iviwiwi 


(3.37) 


Using  (3.37)  in  (3.32)  together  with  the  change  of  variables  u 
(3.32)  reduces  to 

T 


t-s, 


c,t  =  rk<u)Rk(u)Rt(u)du. 


(3.38) 
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From  (3.24),  (3.32),  and  (3.27),  the  variance  of  the  interference 
component  1^  is 


Var  I. 


K 

=  Z 
k=l 


rk(u)Rk(u)Ri(u)du  + 


PT 


12N 


3  rk,i) 


(3.39) 


k^i 

Using  (3.31)  and  (3.39)  in  (3.28),  the  signal-to-noise  ratio  at  the  output 
of  the  i-th  correlation  receiver  for  a  DS/SSMA  system  operating  over  a 
time-selective  fading  channel  is 


N  2y2  T 
i  *  ^2S  +  ~~2  I  r.(u)(T-u)du 


SNR 


T  0 
T 


K  T 

+  2  v?  J  r  (u)R^(u)R  (u)du  +  —  r,  .]}  (3.40) 

k=l  k  0  *  k  1  6NJ  k’1 


k^i 


where  3  =  PT  is  the  energy  per  data  bit. 

Numerical  Results.  The  general  result  (3-40)  is  now  specialized  to  channels 
having  specific  channel  covariance  functions  rk(t)  for  1  <  k  ^  K. 

Throughout  the  remainder  of  this  section  it  is  assumed  that 


r  (t)  -  0  for  |t|  >  XT  (3.41) 

where  X  =  (tvff3),  n  is  a  positive  integer  less  than  N,  P  is  in  the  range 
0  <  P  <  1,  and  X  is  a  positive  number  less  than  N.  The  quantity  XT£  is 
called  the  correlation  duration  for  the  fading  channel.  Note  that  from 
(2.45)  and  (2.60),  the  correlation  duration  for  the  time-selective  fading 


channel  is  equal  to  the  correlation  time  T  for  the  channel  as  defined 
in  Section  2.4,  provided  we  are  using  signals  satisfying  condition  (2.63). 
For  convenience,  the  following  notation  is  introduced  at  this  point. 

Let  A^Oi)  =  C^(£+l)  -  and  tt^  ^(i.)  =  A^(A)C^(A)  .  The  function  ^ 

will  be  denoted  by  tt^ . 

From  (1.10)  and  (3.41),  (3.40)  is  seen  to  imply 


NQ  2y.  T 

SNR  =  [rg  +  — 1/  r  (u)  (T-u)du 

T  o 


K  2  .2  n 


+  s  [  3  \  E  J* 

k-i  r  a=o  k  x 

kH 


W+P,)TC  2 

T  rk(u)du 

n 

c 


+  ["k,i«>  +  "llkct)1  J- 


<i+VTc 

JLT 


Tc (u-£Tc)rk(u)du 


I 


<***»'  ,  ..  v  2 


IT 


(u-XT  )  r  (u)du}  +  — -  r  ]} 
c  K  6NJ  K,1 


(3.42) 


where  0,  *  1  for  l  <  n  and  0  =0.  Observe  that  the  evaluation  of  (3.42) 

it  n, 

requires  only  knowledge  of  the  discrete  aperiodic  autocorrelation  functions 
of  the  K  signature  sequences  as  well  as  knowledge  of  the  channel 
covariance  function  r  (t) .  In  this  chapter,  we  shall  consider  two  types 

K 

of  channel  covariance  functions  to  illustrate  the  general  results.  These 
are  a  triangular  covariance  function  and  an  exponential  covariance  function 


By  a  triangular  covariance  function,  we  shall  mean  a  function  of  the  form 
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A(T-v  t  ) ,  t  <  XT 


rk1}^  -  i 


,  0,  |t|  > 


(3-43) 


where  v  *  (XT£)  T.  A  truncated  exponential  covariance  function  is 


defined  by 


t  <  XT 


r‘2)(t>  - 


(3.44) 


0,  |t|  >  XTc  . 


It  should  be  noted  that  r£  (t)  is  not  a  valid  autocovariance  function, 

but  it  is  an  approximation  to  the  valid  autocovariance  r(t)»  B  exp(-a|t|), 

(2)  (2) 

_a>  <  t  <  “.  It  is  true  that  r^  (t)  =  (-t)  which  is  the  only  property 

of  an  autocovariance  function  that  is  needed  in  what  follows. 

For  the  time- selective  fading  channel,  the  unit  variance  constraint 
on  the  covariance  function  of  the  channel  fading  process  results  in  the 


conditions 


r£1}(0)  -  AT  -  1,  r£2)(0)  -  B  -  1. 


Using  (3.43)  and  (3.45)  in  (3.42),  SNIL  for  a  DS/SSMA  system  over  a 
time- selective  fading  channel  with  a  triangular  covariance  function 
is  found  to  be 


(3.45) 
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2V2X  ,  ,2 

smj  '  -  C-j-  U  -  K'H-D  5  +  ^1 

IN  3N^ 


K  2\  n  p  P1  2 

+  h[~$  ^0c^ck(X)Gi (l)rx<*> +  A  [TTk,i^> + 


*  3  AkW)A1«»'1(t)}]  +  £  A  rk  +  §  }** 

k=l  6N 


(3.46) 


where  r .(C)  =  N-v(j&+cp  ).  In  an  identical  fashion,  it  can  be  shown  that 

X/  Xj 

SNR^  for  a  DS/SSMA  system  over  a  time- selective  fading  channel  with 

covariance  function  r,^2"*  is 

k 


lij>  ■  C-S  K1-'"’')  +  53  <*'“V+i)-i>) 


K  2\.  n  .  -qP  . 

+  2  [  Z  e"aA[C  U)C.U)(l-e  Z) 

k=l  N ct  t=0 
ktt 

-«P*  -aP- 

+  [nkji(£)  +  TTijk(i) ][%[(«£  +  1)  -e  *(a(X  +  PA)  +  1)]  -  X(l-e  *)] 

+  \U)A.a)CU2  +  ^  +  -§  ]-e  X[(*+PA)2  +  f(Jt+Px)  +  1 

a  a 

+  ~  [-(erf  +  1)  +  e  *(cr(A  +  PA)+1]  +  JT(l-e  )}3  ] 

i  K  N  , 

+  ~  2  r.  •  +  yh  . 

6N3  k-1  k>1  ^ 


(3.47) 
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In  Fig.  5,  we  have  plotted  SNR  =  min{SNR^}  for  K  =  6,  using 

equations  (3.46)  and  (3.47),  as  a  function  of  the  correlation  duration  of 

the  channel  correlation  function.  In  evaluating  (3.46)  and  (3.47)  it  has  been 
2  2 

assumed  that  y  =  y  for  all  k,  1  <  k  <  6  and  C~~)  =  10  dB.  The  correlation 
k  Nq 

duration  used  in  evaluating  (3.46)  was  that  value  XT^  satisfying  the 

(2)  -5 

equation  r'  (XT  )  =  e  (i.e.,  aXT  =5).  This  correlation  duration  was 

K  C  C 

selected  numerically  to  minimize  the  effect  of  ignoring  the  tail  of  an 

(2) 

exponential  correlation  function  (see  (3.41)  above).  Thus  r^  is  a  good 
model  for  exponential  correlation.  The  sequences  used  to  evaluate  (3.46) 
and  (3.47)  consisted  of  6  maximal- length  shift  register  sequences  (m-sequences) 
of  period  127  which  collectively  form  a  maximal-connected  set  [Gold  and 
Kopitzke,  1965].  The  shift-register  loading  selected  for  each  m-sequence 
was  the  least-sidelobe  energy/auto-optimal  (LSE/AO)  phase  for  that 
particular  m-sequence  [Pursley  and  Roefs,  1979].  Further  details  on  the 
actual  loadings  used  may  be  found  in  the  first  six  entries  of  Fig.  B.l  of 
[Pursley  and  Roefs,  1979].  In  evaluating  (3.46)  and  (3.47)  we  have  used 


IN 

only  values  of  correlation  duration  satisfying  \  The  reason  for 

N 

considering  this  limited  range  of  correlation  duration  is  that  for  X  >  ~  , 
the  effects  of  adjacent  bits  must  be  taken  into  account  in  the  performance 
analysis  of  time-selective  fading  channels  (note  that  r^^  (t)  and  r^\t) 
are  double- sided  functions  (see  (3.43)  and  (3.44)). 


Correlation  Duration  (Multiples  of  Tr) 

r  V,  FP-6515 

SNR  vs.  correlation  duration  for  a  time-selective  fading  channel  with  triangular  and 
truncated-exponential  covariance  functions  for  two  values  of  y*  (N  =  127,  K  =  6, 

<?/N„  =  lOdB) . 


Note  in  Fig.  5  that  as  the  correlation  duration  increases,  the  DS/SSMA 

system  performance  is  degraded  significantly.  An  explanation  for  this 

degradation  is  as  follows:  as  the  correlation  duration  of  either  r£^ (t) 

(2) 

or  r^  ' (t)  increases,  the  time-selective  fading  channel  model  reduces  to  the 
nondispersive  fading  channel  model  (see  (2.64))  of  Section  2.7,  i.e.,  the  fading 
portion  of  the  transmitted  signal  is  being  multiplied  by  a  complex  random 
variable  as  opposed  to  being  multiplied  by  a  time-varying  random  process .  Since 
the  receiver  structure  being  used  is  a  correlation  receiver,  the  effects  of  the 
smoothing  by  the  integrator  in  the  correlation  receiver  on  the  diffuse 
portion  of  the  received  signals  is  reduced  as  the  channel  correlation 
duration  increases.  Heuristically  speaking,  if  the  fading  process  is  a 
zero-mean,  stationary  ergodic  process  so  that  time  averages  and  probabilistic 
averages  (expectations)  may  be  interchanged,  then  for  small  values  of  correla¬ 
tion  duration,  the  diffuse  portion  of  the  received  signal  is  averaged  over  the 
bit  duration  by  the  receiver  integrator  so  that  its  contribution  to  is 
essentially  zero.  For  larger  values  of  correlation  duration,  the  receiver 
integrator  is  integrating  over  a  time -independent  random  variable,  so  that  the 
contribution  to  of  the  diffuse  portion  of  the  received  signals  becomes 
significant.  As  a  result,  there  is  a  corresponding  decrease  in  SNR^ .  The 
fact  that  a  DS/SSMA  communication  system  performs  better  over  a  time-selective 
fading  channel  with  an  exponential  channel  correlation  function  than  one  with  a 
triangular  channel  correlation  function  may  be  explained  by  examining  (3.40). 
Since  for  a  given  correlation  duration,  the  tail  of  an  exponential  channel 
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correlation  function  falls  off  faster  than  that  of  a  triangular  channel 
correlation  function,  using  (3.45)  we  see  from  (3.40)  that  for  a  given 
value  of  \Tc,  the  contribution  to  of  the  diffuse  portion  of  the 
received  signal  is  smaller  for  the  exponential  channel  correlation 
function. 

Figure  6  illustrates  the  dependence  of  SNR  on  the  transmission 

coefficient  Y  for  a  time-selective  fading  channel  with  a  triangular 

channel  correlation  function.  The  same  signature  sequences  as  described 

above  are  used.  As  might  be  expected  from  (3.46),  SNR  is  degraded  in 

proportion  to  the  intensity  of  the  fading  process.  From  Figs.  5  and  6, 

we  may  ascertain  that  for  a  given  minimum  SNR  at  the  output  of  a 

correlation  receiver,  there  exists  an  admissible  range  of  y  and  \Tc  for 

which  DS/SSMA  communications  is  possible.  Suppose,  for  example,  we  use 

the  LSE/AO  m-sequences  described  above  for  DS/SSMA  communications  over  a 

time- selective  fading  channel  with  a  triangular  channel  correlation 

2 

function.  A  desired  minimum  SNR  of  10  dB  is  possible  provided  y  <0.2 

and  the  correlation  duration  of  the  channel  is  less  than  20  T  . 

c 

If  a  receiver  SNR  higher  than  that  obtainable  from  the  fading  channel 

for  a  given  5/Nq  is  required,  one  would  expect  that  by  increasing  $/Nq, 

the  desired  performance  could  be  obtained.  That  this  approach  ma/  be 

fruitless  is  illustrated  by  Fig.  7  in  which  SNR  is  shown  as  a  function 

2 

of  <S /NQ  for  a  time-selective  fading  channel  for  three  values  of  y  .  In 
Fig.  7,  the  correlation  duration  is  10  T£  and  the  set  of  signature 


lective  fading  channel  with  a  triangular  covariance 
elation  duration  (N  =  127.  K  =  6,  3/ =  lOdB) . 


lective  fading  channel  with  a  triangular  covariance  function 
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sequences  is  the  same  as  used  for  Fig.  1.  The  observed  leveling-off 
effect  of  the  SNR  vs.  <S/N q  curve  may  be  expected  for  any  communication 
system  in  which  the  expression  for  Var  contains  a  sufficiently 
large  interference  component  whose  power  is  directly  proportional  to  the 
transmitted  signal  power  (see  (3.31)  and  (3.39)).  Hence,  this  effect  will 
occur  for  values  of  K  and  N  in  the  range  of  interest;  even  for  SSMA 
communication  via  an  AWGN  channel. 

In  the  three  figures  immediately  above  we  have  chosen  to  use  a 
maximal-connected  set  of  m-sequences  with  LSE/AO  phases  as  the  signature 
sequences  without  giving  any  justification  for  this  choice  of  signature 
sequences.  Some  justification  for  this  choice  will  now  be  given.  First, 
we  have  used  a  maximal-connected  set  of  m-sequences  because  such  a  choice 
yields  the  largest  possible  set  of  m-sequences  of  period  N  for  which  any 
two  sequences  in  the  set  have  a  preferred  three-valued  crosscorrelation 
function  [Gold  and  Kopitzke,  1965].  Hence,  by  using  such  a  set  we 
automatically  are  using  a  set  whose  sequences  have  bounded  pairwise 
periodic  crosscorrelation  values,  a  desirable  feature  for  synchronization 
purposes  in  DS/SSMA  systems  and  for  communication  in  synchronous  DS/SSMA 
systems.  Second,  we  have  used  the  LSE/AO  phases  for  the  individual 
m-sequences,  primarily  as  a  sieve  over  the  phases  of  the  m-sequences,  but 
secondarily  for  the  reason  that  such  a  choice  offers  promise  for  good  per¬ 
formance,  at  least  over  AWGN  channels  [Golay,  1972].  From  [Pursley  and  Sarwate, 
1977b],  we  note  that  the  interference  parameter  r,  .  is  bounded  by 

K,  1 


4N2  -  6(S(a(k))S(a(i))}^  <  r  <  4N2  +  6{s (a(k) )S (a(i) ) 

K,  1 


(3.48) 
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where  S (u)  is  the  sidelobe  energy  of  sequence  u  defined  by 

N-l  2 

S  (u)  =  Z  C  CJL).  (3.49) 

£-1 

Hence  by  bounding  the  sidelobe  energy  of  each  of  the  signature  sequences 

through  the  use  of  LSE/AO  phases,  we  are  bounding  the  value  of  r, 

K  ,  1 

thereby  bounding  the  performance  of  a  DS/SSMA  system  over  an  AWGN 
channel  (see  Eq.  (1.19)).  In  fact,  the  choice  of  LSE/AO  phases 
is  not  necessarily  critical  for  fading  channels.  For  example,  if 
AO/LSE  phases  [Pursley  and  Roefs,  1979]  were  used,  the  results  obtained 
for  Fig.  5  would  differ  from  the  results  obtained  using  AO/LSE  phases 
by  less  than  0.004  dB.  In  Section  3.5  we  will  see  that  maximal- 
connected  sets  of  m-sequences  with  LSE/AO  phases  perform  very  much  like 
"typical1'  sequences.  Thus  it  is  expected  that  the  numerical  results 
presented  above  in  Figs.  5-7  are  typical  results  that  would  be  obtained 
using  any  reasonable  set  of  6  signature  sequences  of  length  127.  However 
as  will  be  shown  later,  bad  choices  of  signature  sequences  will  give  far 
worse  performance. 

3.3.  DS/SSMA  System  Performance  for  Frequency-Selective  Fading  Channels 
Analysis.  The  dual  of  the  time-selective  fading  channel  is  another 
special  case  of  a  WSSUS  channel  known  as  a  frequency-selective  fading 
channel.  In  this  section  the  effect  of  frequency-selective  fading  on  the 
performance  of  DS/SSMA  systems  is  considered. 

The  channel  correlation  function  for  a  frequency-selective  fading 
channel  is  defined  (as  in  (2.70))  by 


g.(T)  =  Pi(T,0). 


(3.50) 
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An  additional  assumption  about  the  selectivity  of  the  channel  is  necessary 
in  order  to  obtain  useful  results  regarding  the  performance  of  DS/SSMA 
systems  over  frequency-selective  fading  channels.  As  can  be  seen  from 
the  equations  that  follow,  a  frequency-selective  fading  channel  exhibits 
memory  and  therefore  introduces  inters,  mbol  interference  into  the  received 
signal.  We  will  assume  in  this  section  that  the  selectivity  of  the  channel 
is  such  that,  in  the  detection  of  a  given  data  symbol,  we  need  only  be 
concerned  with  the  intersymbol  interference  due  to  the  two  adjacent  data 
symbols.  This  condition  is  equivalent  to  assuming  that 

gi(T)  m  0  for  |t|  >  T.  (3.51) 

Channels  exhibiting  higher  degrees  of  selectivity  will  require  the  inclusion 
of  the  effects  of  more  than  the  two  adjacent  bits  in  the  analysis  of  the 
performance  of  the  system,  in  which  case  the  analysis  that  follows  can  be 
modified  in  a  straightforward  manner.  However,  as  noted  in  [Bello  and 
Nelin,  1963],  the  need  to  include  more  than  two  adjacent  bits  in  the 
system  analysis  is  frequently  an  indication  that  the  channel  is  too 
frequency- selective  for  use  in  practical  communication  systems. 

Using  (3.50)  and  (3.51)  in  (3.23),  we  find  that  the  variance  of  the 
faded  portion  of  the  i-th  signal  for  a  frequency-selective  fading  channel 
is 

p  2  T  T  T 

Var  F  =  t  y  J*  g  (t)  J  J  f .  (t,t  ;t,s)dtdsdT  .  (3.52) 

-T  00 

Notice  from  (3.11)  that  the  double  integral  with  respect  to  t  and  s  in 
(3.52)  reduces  to  the  product  of  two  single  integrals,  each  of  which  is 
given  by 


80 


(  bi,oVT)  +  bi,-iVT)’  T  s  ° 


I  h.  (T  ;  t)dt  = 
0  1 


(3.53) 


V  b.^R.C.T)  +  b.^R.C-T),  T  <  0 


where  R^(t)  is  Che  continuous- time  partial  autocorrelation  function 
defined  in  Chapter  1.  Upon  substitution  of  (3.53)  into  (3.52),  we 
obtain 


Var  F.  =  PyJ  J  g£ CO  [R? (t )+[b.  ^ Qb.  # _1+b .  , Qb.  ^  ^R.  (T)R.(r)+Rr (t)  ]dr  .  (3.54) 

Averaging  (3.54)  over  all  data  patterns,  it  is  easy  to  see  that  the 
variance  of  the  faded  portion  of  the  i-th  signal  is 

T 


Var  F.  =  PY.  J  g.  (t)[R.  (t)  +  R  (T)]dx. 

x  x  0  i  i 


(3.55) 


The  variance  of  the  interference  component  1^  may  be  evaluated  for 
a  frequency-selective  fading  channel  by  using  (3.50)  and  (3.51)  in  (3.25) 
to  obtain 

T  T  T 


V.  .  -  J  J  J  g .(r)a  (t)a  (s)  X(t,s,T)dTdtds 
0  0  -T  K 


where 


X(t,s,T)  =  ^  J’o%(t“Tk"T)ak(s"Tk'T)E^bk(t:_Tk"T)bk(s"Tk"T)^dTk'  (3.57) 

Letting  u  =  t-T  and  v  =  s-T  in  (3.57),  we  note  from  (3.33)  that 


(3.56) 


ai(t)ai(s)X(t,s,T)  =  E^fk}kji(T^k»T+Tk>t’s^» 


(3.58) 
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and  hence  using  (3.37),  (3.57)  reduces  to 


X(t,s,T)  =  ^  R_k(  |  t-s  | ) 


(3.59) 


Using  (3.59),  (3.56)  may  now  be  rewritten  as 

T  T  T 

..  :  j 

•T  0  0 

After  a  change  of  variables  T  =  t-s,  (3.60)  reduces  to 
T  T 


1  T  ,T  ,T 

vk,i  =  T  >f  sk(T)dT  J"  I  ai(t)ai(s)Rk(|t-s|)dtds. 


vk,i  =  t{t  gk(T)dT  J*0  VT)VT)dT 


(3.60) 


2T3  T 

=  It  rk.i  J0 


(3.61) 


where  the  second  step  follows  from  [Pursley  and  Sarwate,  Section  V,  1977a] 
and  the  last  step  follows  from  [Pursley,  eqs .  (11)-(12),  1977].  From 
(3.24),  (3.61),  and  (3.27),  the  variance  of  the  interference  component 


I.  is 
l 


Var 


3 

K  T  T  2 

h  '  p  s  rk  i[vklT  I  gk(T)dT  +  ! 

i  k*l  0  12Nd 


(3.62) 


k^i 


Using  (3.55)  and  (3.62)  in  (3.28),  the  signal-to-noise  ratio  at  the 
output  of  the  i-th  correlation  receiver  for  a  DS/SSMA  system  operating 
over  a  frequency- selective  fading  channel  is 
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n  2y7  r 
S®i  ■  C  25+^I0 


K 


i:  CO[R2(T)  +  R2(T)]dT 

J.  1  1 


T3  T 


2  E  rk  i[\  3T  J*  §k(T)dT  +  - 3 

k=l  ’  *  0  12N 


+  -4  Sr. 
T 

k^i 


-h 


(3.63) 


Numerical  Results .  In  this  subsection  (3.63)  is  specialized  to  the  case 
of  a  specific  channel  covariance  function  g^(T).  In  all  that  follows  it 
is  assumed  that  g,  (t)  satisfies  (3.41);  that  is  g^CO  is  non-zero  only 
over  the  interval  |t|  <  XT .  Analogous  to  the  time -selective  fading 
channel  case,  r,e  define  XT  to  be  the  correlation  duration  of  the  frequency- 
selective  fading  channel.  From  (2.77)  and  (2.78)  we  see  that,  provided 
(2.82)  is  satisfied,  the  first  zero-crossing  of  the  Fourier  transform  of 
g.  (T)  is  defined  to  be  the  correlation  bandwidth  of  the  frequency- selective 
fading  channel.  In  practice,  however,  the  correlation  bandwidth  of  the 
channel  is  usually  defined  in  this  manner  without  taking  into  account  the 
bandwidth  of  the  transmitted  signal  ([Monsen,  1971],  [Bello  and  Nelin, 
1963]).  Thus  from  the  theory  of  Fourier  transforms,  if  a  frequency- 
selective  fading  channel  has  a  correlation  duration  XTc>  the  correlation 
bandwidth  of  the  channel  is  proportional  to  (XT£)  ^  (see  [Monsen,  1971]) . 
Consequently,  the  correlation  duration  of  the  channel  is  also  a  measure 
of  the  correlation  bandwidth  of  the  channel. 

Using  (1.9),  (1.10),  and  (3.41)  in  (3.63),  we  find  that  for  a 
frequency- selective  fading  channel  with  a  specified  channel  covariance 
function  which  satisfies  (3.41),  SNR^  is 
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c  ly  T 

:  rk  i[”  J*  8k(T>dT  +  3} 

=1  3N  0  K  6tT 


<1+2\>rk,i 


(3.66) 


Triangular  and  truncated  exponential  covariance  functions  for  a 
frequency-selective  fading  channel  are  defined  in  a  manner  analogous  to 
tne  definition  of  the  corresponding  channel  covariance  functions  for  a 
time- selective  fading  channel  (i.e.,  g^  is  defined  as  in  (3.43)  and  (3.44)). 
For  this  channel  model  however,  (3.65)  implies  that  the  normalization 
factors  A  and  B  appearing  in  (3.43)  and  (3.44)  are  given  implicitly  by 


84 


2  -a\T 

.  T  2B  ri  ci  -  i 

A  —  =  —  [1  -  e  ]  -  1, 

v  a 


(3.67) 


where  X  is  defined  in  (3.41).  For  a  frequency-selective  fading  channel 

with  a  triangular  channel  covariance  function,  SNR^  is  found  to  be 

/in  n  .  o  -I 

SNR.U;  =  C— 3  2  [P ^  (cf (f -N)  +  C^a)}rA(i) 

YN  x=o 

g3 

+  p]Ctt.(ji-n)+tt  .(f) }rx(-|)  +  -j  CaJu-n)  +  Aj(X)}rx(|)] 


1  K  ?  N  , 

+  A  2  (1+2Vrk  i  +  zl  }  * 

6N3  k=l  k  k,i  20 


(3.68) 


The  corresponding  result  for  a  frequency-selective  fading  channel  with 
a  truncated-exponential  channel  covariance  function  is 


SNRi2)  -  +  T' '^"ai  '  E  e'aA[(cJ(£-N)  +  C^CO)(l-e  ^ l) 

1  2<S  N2[l-e'aX]  A=0 

1  "«P  ,  -°P . 

+  2(n.  U-N)  +  it  (i))(i[aX+l  -e  VU+0  ,)+l)  }-l  (1-e  *)) 

1  i.  U  Xj 

2  2  2  21  2  2  3(£+p.)  o 

+  (A.a-Nj+Ajumir  +  ~  +  -fi-e  ^tu+p^r  +  — +  -ji 

a  a 

2 ,  2  "aPf 

+  “ [e  *(a(£+Px)+l)  -  «a+l)]  +  JT(l-e  *))] 


6N  k*l 


(3.69) 
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Equations  (3.68)  and  (3.69)  were  evaluated  as  a  function  of  the 
correlation  duration  using  the  same  code  sequences  as  for  Fig.  5,  and 
the  results  are  shown  in  Fig.  8.  As  in  the  time- selective  fading  case, 
it  is  assumed  that  Yk  =  Y  for  k,  1  <  k  <  6,  and  <S/Nq  =  10  dB  .  as 
before,  the  correlation  duration  for  a  triangular  channel  covariance 
function  is  defined  to  be  XTc-  However,  for  a  truncated  exponential 
channel  covariance  function,  we  define  the  correlation  duration  to  be 
the  value  of  \Tc  satisfying  the  equation  g^(XT^)  =  10  ^  (i.e.,  a  satisfies 
-a\T  .  -aXT  c 

§a[l-e  ]  *e  =  10  ).  Once  again,  this  particular  definition  of 

correlation  duration  was  selected  to  minimize  the  effect  of  truncating  &  (t) 
From  Fig.  8,  we  note  that  for  small  values  of  XT^,  SNR  for  a  frequency- 
selective  fading  channel  is  degraded  significantly.  Because  the  frequency- 
selective  fading  channel  is  the  dual  of  the  time-selective  fading  channel, 
as  was  pointed  out  in  Section  2.6,  such  a  result  could  have  been  predicted 
directly  from  Fig.  5  using  a  time- frequency  duality  argument.  Physically, 
what  is  occurring  is  that  as  the  correlation  duration  of  the  channel  goes 
to  zero,  the  correlation  bandwidth  of  the  channel  approaches  infinity.  From 
(2.84),  we  see  that  for  small  values  of  correlation  duration,  the  frequency- 
selective  fading  channel  model  may  therefore  be  replaced  by  the  non- 
dispersive  fading  channel  model  of  Section  2.7.  Following  an  argument 
identical  to  that  used  for  this  limiting  situation  in  the  time-selective 
fading  channel  case,  we  see  that  small  values  of  correlation  duration 
lead  to  a  reduction  of  the  receiver  integrating  effect  on  the  diffuse 
portion  of  the  i-th  signal.  Note  that  the  contribution  to  of  the 
diffuse  portion  of  the  K-l  interfering  signals  is  unaffected  by  the 
correlation  duration  of  the  channel,  as  can  be  seen  from  (3.66). 
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Figure  9  illustrates  the  dependence  of  the  signal- to-noise  ratio 
on  the  square  of  the  transmission  coefficient  y  for  a  frequency-selective 
fading  channel  with  a  triangular  covariance  function  using  the  same 
signature  sequences  as  used  for  Fig.  8.  As  might  be  expected,  the 
signal-to-noise  ratio  is  degraded  in  proportion  to  the  intensity  of  the 
fading  process-  From  Figs.  8  and  9,  we  may  ascertain  that  for  a  given 
signal-to-noise  ratio  to  be  achievable  the  parameters  y  and  XT  must  be 
constrained  to  a  certain  range. 

The  graph  of  SNR  vs.  S/N^  for  a  frequency-selective  fading  channel 
with  a  triangular  channel  covariance  function  is  shown  in  Fig.  10.  As 
in  Fig.  7,  the  SNR  curve  "f lattens-out"  as  <S/Nq  increases  for  the  same 
reasons  as  given  for  the  time-selective  fading  channel. 

The  comments  made  in  the  last  section  about  the  choice  of  signature 
signatures  used  to  evaluate  expressions  (3.68)  and  (3.69)  are  valid  for 
the  frequency-selective  fading  case,  also.  In  particular,  had  the  AO/LSE 
phases  of  the  m-sequences  been  used  instead  of  the  LSE/AO  phases,  the 
results  obtained  would  differ  from  the  results  presented  in  Figs.  8-10 
by  less  than  0.006  dB . 

In  the  discussion  up  to  this  point,  we  have  not  indicated  which  of 
the  two  singly-spread  channels  is  a  better  model  for  actual  radio 
channels,  especially  channels  over  which  DS/SSMA  communication  systems 
are  to  be  used.  We  will  now  correct  this  omission.  In  practice,  most 
fading  channels  exhibit  some  degree  of  frequency-selectivity;  in 
particular,  multipath  channels  are  frequency-selective  fading  channels 
(e.g.,  see  [Stein,  pp .  351-355,  1966]).  In  addition,  because  the  bandwidth 
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for  a  frequency-selective  fading  channel  with  a  triangular  covariance  function 
values  of  correlation  duration  (N  =  127,  K  =  6,  8/ =  lOdB) . 


occupied  by  spread-spectrum  systems  (including  DS/SSMA  systems)  is 
typically  quite  large,  the  channels  used  by  spread-spectrum  systems  often 
appear  frequency- selective  to  the  spread-spectrum  signal,  even  though 
the  same  channel  would  not  be  frequency-selective  to  a  non-spread-spectrum 
signal  transmitting  the  same  data  (see  (2.84)). 


From  the  above  we  may  conclude  that  while  the  time- selective  fading 
channel  may  be  a  good  model  for  some  radio  channels,  the  frequency- 
selective  fading  channel  model  is  of  more  practical  value  when  discussing 
the  performance  of  direct-sequence/spread-spectrum  communications  via 
multipath  channels.  Hence,  in  the  remainder  of  this  chapter,  emphasis 
will  be  placed  on  frequency- selective  fading  channels. 

3.4  Comparison  of  Single-User  PSK  and  DS/SS  Communications  via 
Fading  Channels 

In  Chapter  1  it  was  noted  that  spread-spectrum  techniques  could,  with 
the  proper  choice  of  modulation  technique,  be  effective  in  combating 
the  effects  of  fading  encountered  by  conventional  modulation  systems 
(e.g.,  PSK,  FSK)  over  communications  channels.  In  surveying  the 
literature,  we  find  many  references  which  note  in  particular  that  use  of 
direct-sequence/spread-spectrum  (DS/SS)  modulation  is  effective  in 
combating  multipath  and  the  associated  frequency-selectivity  (see 
Section  3.3)  present  on  many  channels  (e.g.,  see  [Cahn,  1974], 


[Jacobs,  1975],  [Massey  and  Uhran,  1975],  and  [Viterbi,  1979]).  Although 
several  "heuristic"  proofs  of  this  conjecture  exist  ([Cahn,  1974],  [Jacobs, 
1975]),  to  date  no  rigorous  demonstration  of  this  conjecture  has  been  given 
In  the  following,  we  will  give  such  a  demonstration.  In  particular,  we 
will  compare  the  performance  of  single-user  PSK  and  DS/SS  communication 
systems  over  singly-spread  Rician  fading  channels.  Throughout  this 
comparison,  average  signal-to-noise  ratio  at  the  output  of  a  correlation 
receiver  will  be  used  as  the  performance  measure. 

From  (3.40),  we  see  that  for  a  single-user  DS/SS  system  operating 
over  a  time- selective  fading  channel,  the  average  signal-to-noise  ratio 
at  the  output  of  a  correlation  receiver  is  given  by 

N  2  T 

SNR  »  r (u)  (T-u)du}”^  .  (3.70) 

T  O 

Since  the  expression  for  average  SNR  given  by  (3.70)  does  not  depend  on 
the  signature  sequence  used  by  the  DS/SS  system,  we  may  conclude  that 
equation  (3.70)  is  also  the  expression  for  average  SNR  for  a  biphase 
PSK  system  operating  over  a  time- selective  fading  channel.  For  if  the 
signature  sequence  (a^^3  is  identically  one  for  -®  <  n  <  ®  in  (1.2),  the 
DS/SSMA  system  model  shown  in  Fig.  1  reduces  to  a  PSK  system,  when  K  =  1. 
From  the  above  discussion  we  may  conclude  that  so  far  as  single-user 
systems  are  concerned,  the  performance  of  PSK  and  DS/SS  systems  are 
identical  over  time -selective  fading  channels. 

For  frequency- selective  fading  channels  however,  the  result  of  this 
comparison  of  the  two  systems  is  quite  different.  From  (3.63),  the  average 
SNR  at  the  output  of  a  correlation  receiver  for  a  single-user  DS/SS 
system  operating  over  a  frequency-selective  fading  channel  is 


(3.71 


SNR 


_0 

2(5 


g(T)[R:(T)+Rf(T)]dr}‘2 

T  O 


The  average  SNR  for  a  PSK  system  operating  over  a  frequency-selective 
fading  channel  may  be  found  either  directly  from  (3.52)  by  letting  a^(t)  = 
in  (3.11)  or  by  noting  that  for  a  PSK  system,  R^(t)  =  t  and  R^(r)  =  (T-r), 
and  using  (3.71).  Both  of  these  results  follow  from  the  comments  made  in 
the  previous  paragraph  about  the  reduction  of  the  DS/SS  system  model  to 
the  PSK  system  model.  Using  either  method,  the  average  SNR  for  a  PSK 
system  is  found  to  be 


SNR 


g (T ) [T2+(T-T)2]dT  . 


(3.72 


In  Fig.  11  we  have  evaluated  eqs.  (3.71)  and  (3.72)  as  a  function  of 

2  ' 

the  correlation  duration  of  the  channel  for  two  values  of  y  and  for 

S 

—  =  10  dB .  In  evaluating  (3.71),  we  have  used  an  m-sequence  of  length 
N0 

127  with  its  LSE/AO  phase  as  the  signature  sequence  for  the  DS/SS  system. 

Thus  for  the  PSK  system,  the  correlation  duration  scale  is  interpreted  as 

being  marked  in  intervals  of  yjy  of  a  bit.  The  relatively  large  advantage 

2 

of  DS/SS  over  PSK  (e.g.,  6.04  dB  for  X  =  40,  y  =0.2)  may  be  thought  of 
as  frequency  diversity  of  spread-spectrum  communication.  Most  of  this 
improvement  results  from  the  multipath  rejection  capability  of  a  properly 
designed  spread-spectrum  system.  It  should  be  noted  that  this  multipath 
rejection  is  fully  achieved  only  if  the  signature  sequences  have  good 
aperiodic  autocorrelation  properties.  As  will  be  seen  below,  poor  choices 
of  signature  sequences  will  result  in  a  signal-to-noise  ratio  that  is  as 
much  as  3  dB  below  the  DS/SS  curve  of  Tig.  11. 


SNR  vs.  correlation  duration  for  a  frequency-selective  fading  channel  with  a 
triangular  covariance  function  for  single-user  DS/spread-spectrum  and  biphase 


We  will  now  examine  Fig.  11  in  greater  detail.  First,  from  (3.41) 
and  (3.65),  for  small  values  of  correlation  duration,  g(T)  approaches 

A 

a  delta  function  of  unit  area.  Since  R^(0)  =  0  and  EL  (0)  =  T  from  (1.7) 
and  (1.8),  we  see  that  for  a  correlation  duration  of  value  zero,  (3.71) 
and  (3.72)  both  reduce  to 

Nn  2  jl 

SNR  =  +  Y  }"*  (3.73) 

i.e.,  the  performance  of  the  two  systems  is  identical  for  XT  =0.  From 

c 

(3.73),  we  see  that  for  XT£  =  0,  SNR  =  8.239  for  =  0*1  and  SNR  =  6.020  dB 
2 

for  y  =  0.2  for  both  systems.  This  convergence  of  performance  at  X  =  0 
as  well  as  the  performance  of  both  types  of  systems  at  other  values  of  XTc 
can  be  explained  physically  by  considering  the  autocorrelation  function  of 
the  signature  sequence  together  with  the  concept  of  multipath.  To  see 
this,  note  that  the  continuous -time  periodic  autocorrelation  function 
^(t)  of  a  code  waveform  a^(t)  is  defined  as 

A  T 

0.(T)iJ'  a.  (t-T)a.(t)dt  (3.74) 

-  R. (t)  +  R. (T),  (3.75) 

where  the  last  line  follows  from  (1.7)  and  (1.8).  Using  (3.75)  in  (3.71), 
we  see  that  SNR  for  a  DS/SS  system  operating  over  a  frequency-selective 
fading  channel  becomes 

N  2  T 

SNR  =  fcS  +  ^J*  g(T)[92(T)  -  2R  (T)R  (T)]dT}“* 

i  J.  i. 


(3.76) 


which  is  a  function  of  both  the  aperiodic  and  periodic  autocorrelation 
functions  of  the  code  waveform.  From  the  discussion  in  Section  2.6, 
g(T)  appearing  in  (3.76)  may  be  considered  to  be  a  weighted  "window- 
function"  on  the  allowable  range  of  delay-spread  of  the  multipath  of  the 
fading  channel.  Thus  for  a  given  value  of  correlation  duration,  g(f) 
"passes"  only  multipath  having  range  delays  less  than  X.T .  Consequently, 
multipath  present  on  the  channel  will  degrade  the  performance  of  a  DS/SS 
system,  the  extent  of  degradation  being  determined  by  the  exact  form  of 
the  continuous- time  autocorrelation  functions  of  the  code  waveform. 

To  give  a  concrete  example  of  how  the  shape  of  the  autocorrelation 
functions  of  the  code  waveform  affect  the  performance  of  the  DS/SS  system, 
assume  that  the  DS/SS  system  employs  the  N  3  120  sequence  given  by 
{a^3  3  {l,  1, 1,0, 0,0, . . .) .  The  magnitude  of  the  continuous -time  periodic 
autocorrelation  function  of  a^t)  is  shown  in  Fig.  12.  From  Fig.  12  and 
(3.76)  we  note  that  multipath  signals  having  delays  which  are  odd  multiples 
of  3Tc/2  will  be  rejected  by  the  DS/SS  receiver.  On  the  other  hand, 
multipath  signals  having  delays  which  are  even  multiples  of  3T Jl  will  not 
be  rejected  by  the  receiver,  thus  causing  degradation  of  the  DS/SS  system 
performance.  Multipath  signals  having  delays  of  other  values  will  also 
tend  to  degrade  the  system  performance,  the  extent  of  degradation  being 
determined  by  the  value  of  |9(t)|  for  a  particular  value  of  delay  t.  Now 
consider  the  performance  of  the  DS/SS  system  when  an  m-sequence  of  period  127 
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T 
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Figure  13.  Magnitude  of  the  autocorrelation  function  of  an  m-sequence 
(N  =  127). 
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is  used  for  the  signature  sequence.  The  magnitude  of  the  continuous -time 

periodic  autocorrelation  function  of  this  signature  sequence  is  shown  in 

Fig.  13  [Stiffler,  1971].  In  this  case,  because  of  the  "noise-like"  periodic 

autocorrelation  function  of  the  m-sequence,  all  multipath  signals 

arriving  with  delays  greater  than  one  chip  duration  will  tend  to  be 

rejected  by  the  receiver,  as  can  be  seen  from  (3.76).  Because  of  its 

good  periodic  autocorrelation  function  an  m-sequence  is  expected  to  be  a 

better  candidate  for  a  signature  sequence  than  the  "periodic"  signature 

sequence  discussed  above,  when  system  performance  is  being  evaluated  for 

a  frequency- selective  fading  channel.  In  Fig.  14  we  have  evaluated  (3.76) 

for  three  "periodic"  signature  sequences  of  length  120  having  periods  of 

2,  6,  and  12,  respectively,  and  for  one  m-sequence  of  period  127  for  a 

&  2 

triangular  channel  covariance  function  and  for  —  =  10  dB,  y  =0.2.  From 
this  figure  we  see  that  our  conjectures  about  the  performance  of  the  two  types 
of  signature  sequences  are  valid.  For  values  of  X  greater  than  two,  the 
m-sequence  outperforms  all  other  signature  sequences,  primarily  due  to  the 
periodicity  of  the  periodic  autocorrelation  functions  of  the  other  signature  se¬ 
quences.  Note  that  for  X  »  20,  the  difference  in  performance  of  the  two  types  of 
signature  sequences  is  2.87  dB .  The  reason  that  all  of  the  SNR  curves  exhibit  a 
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general  upward  trend  is  due  to  the  unit  variance  constraint  (3.65)  on  the 
channel  covariance  function  (t )  .  For  large  values  of  correlation 
duration,  g^ (t )  is  "spread-out"  over  the  range  0  ^  t  <  XTc  and  since  the 
area  under  g^(r)  is  constant,  the  magnitude  of  g^(T)  decreases  for 
increasing  values  of  XT^,  resulting  in  a  smaller  contribution  of  the 
integral  term  appearing  in  (3.76). 

The  above  explanation  has  centered  upon  the  periodic  autocorrelation 
function  of  the  code  waveform  as  the  basis  for  DS/SS  system  performance, 
without  taking  into  account  the  aperiodic  autocorrelation  terms  which 
also  appear  in  (3.76).  In  addition,  we  have  neglected  the  "smoothing" 
effects  of  the  integral  appearing  in  (3.76).  Thus  the  explanadon  given 
above  should  be  treated  as  a  "first-order"  model  of  what  is  physically 
occurring  in  the  DS/SS  system.  Nevertheless,  this  first-order  model  is 
accurate  enough  to  predict  the  performance  of  a  PSK  system  over  a 
frequency-selective  fading  channel.  Treating  the  PSK  system  as  simply 
a  DS/SS  system  with  a  constant-code  waveform,  the  periodic  autocorrelation 
function  of  the  PSK  code  waveform  is  a  constant  function.  Thus  the  PSK 
system  does  not  reject  any  multipath  signals,  resulting  in  its  inferior 
performance  compared  to  a  DS/SS  system  using  an  m-sequence  for  its 
signature  sequence. 

3.5.  Performance  Evaluation  for  Random  Signature  Sequences 

In  this  section  approximate  expressions  are  obtained  for  SNR^  for 
both  time- selective  and  frequency-selective  fading  channels  with  a 
specific  channel  correlation  function.  These  approximations  are  the 
expected  values  of  SNR^  when  the  signal  re  sequences  are  random  binary 


sequences.  By  random  binary  sequences  we  mean  a  sec  of  K  statistically 

independent  binary  sequences  1  <  k  <  K.  Each  sequence  {a^}  is 

(k)  r  (k)  i 

a  sequence  of  independent  random  variables  aK  '  for  which  Pr{a '  =  +1]  = 

n  n 

Pr(a^k)  s  -l}  =  The  expected  value  of  SNR^  is  just  the  average 
signal-to-noise  ratio  SNR,  where  the  averaging  is  over  all  possible  sets 
of  K  binary  sequences  of  period  N.  One  reason  for  computing  SNR  is  that, 
in  some  sense,  SNR  is  a  measure  of  the  asymptotic  performance  of  a  DS/SSMA 
system  for  which  the  signature  sequence  length  N  is  very  large  [Roefs,  1977] 
A  second  reason  for  computing  SNR  is  that,  as  we  will  show  below,  SNR  is 
a  close  approximation  (for  a  reasonable  set  of  signature  sequence)  to 
SNR^,  yet  evaluation  of  SNR  does  not  require  computation  of  the  aperiodic 
correlation  functions  of  any  specific  sets  of  signature  sequences. 

First  we  note  that  for  random  binary  sequences  the  moments  of  the 
aperiodic  autocorrelation  functions  are  given  by 


E(CkU)}  =  { 


E{Ck(£)}  = 


)  0,4*0 
N ,  4  =0 


N-|jt|  ,4*0 
2 


N 


,4=0 


E(cka)ci0o3 


(*0,4*0,  k*i 


N  ,  4=0 


(3.77a) 


(3.77b) 


(3.77c) 


E{Ck(4)Ci(4+l)3  =  0  Vk,i  . 


(3 . 77d) 


Equations  (3 . 77a) - (3 . 77d)  were  previously  derived  by  Roefs  (1977)  except 
for  the  case  l  =  0  in  (3.77a)  through  (3.77c).  From  the  definitions  of 
C^OO  and  random  binary  sequences,  for  i  =  0  we  obtain 


N-l 

E(c,  (0)}  =  e(  z 

j=0 


[ajk)]2 


N-l 

e[i  i}  =  n. 

j=0 


(3.78) 


Similarly  we  find  that 


N-l 

(0)]  =  E{  Z 

3=0 


N-l 

Z 

h=0 


I,"']2}  -  »2 


Vk,i  . 


(3.79) 


Using  (3.77),  we  may  evaluate  the  expected  values  of  the  code-dependent 
terms  in  (3.46)  and  (3.47).  Let  SNR  =  T  V%P  [E[var  Ft  +  Var  1^  +  ^NqT]"^, 
where  Var  F^  and  Var  1^  are  now  random  variables  because  the  signature 
sequences  [a^)]  are  random.  Using  the  above  results,  we  find  that  for  a 
time- selective  fading  channel  with  a  triangular  channel  covariance 
function,  SNR  is 


Hi-  U  -  1|  1  »2-^3  +  ^  +  S  I'4 


(3.80a) 


for  X  <  1  and 
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(3.80b) 


for  X  2  1.  The  corresponding  result  for  a  time-selective  fading  channel 
with  a  truncated-exponential  channel  covariance  function  is 


SNR  =  C  ^  [f  N  +  ^  (4  -  f)  -  -  +  (n-1)  (3  -  ^ 
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(3.82c) 


for  2  <  \  <  N-l.  The  corresponding  result  for  a  frequency-selective 
fading  channel  with  a  truncated-exponential  channel  covariance  function  is 


SNR  =  c - 1—r  [(l-e"“P)  -  l  (l-e'^Cctf  +  1)) 
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(3.83a) 


for  X  <  1, 
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for  1  ^  X  <  2,  and 


SNR  =  C - — -  [N2(l-e_Qf  -  z  +  z  e~a (a  +  1)) 
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for  2  <  X  <  N-1.  Note  that  expressions  (3.80)  through  (3.83)  for  SNR 

2 

are  functions  only  of  the  channel  covariance  function,  the  parameter  y  , 
the  number  of  users  K,  the  length  of  the  signature  sequence,  and  the  bit 
energy  to  noise  spectral  density. 


In  order  to  ascertain  the  accuracy  of  the  above  results  as  an 

approximation  to  (3.46),  (3.47),  (3.68),  and  (3.69),  SNR  was  evaluated 

2 

for  various  values  of  V  and  XT  for  N  =  127,  K  =  6,  and  these  results 

c 

were  compared  with  those  obtained  in  Sections  3.2  and  3.3  for  the 

corresponding  channel  type  and  channel  correlation  function.  For  the 

ranges  of  y^ ,  XTc,  and  6/N q  used  in  Figs.  5  and  8,  the  expression  for  SNR 

differed  from  the  exact  expression  for  SNR^  by  less  than  0.05  dB,  for  the 

specific  signature  sequences  used  to  obtain  Figs.  5  and  8.  As  a  further 

check  on  the  accuracy  of  the  above  results,  (3.46)  was  evaluated  for  a 

maximal- connected  set  of  m-sequences  using  LSE/AO  phases  [Pursley  and 

Roefs,  1979]  for  N  =  31  and  N  =  63  and  the  results  were  compared  with  the 

2 

corresponding  results  obtained  from  (3.80).  Once  again,  for  y  in  the 
range  0.05  ^  0.2  and  X  5:  0.1N  the  difference  between  SNR  and  SNR.^  was 
found  to  be  less  than  0.05  dB.  Thus,  based  on  this  comparison  of 
numerical  results,  we  conclude  that  SNR  is  a  very  good  approximation  to 
SNR^  for  the  sets  of  signature  sequences  used  above.  Furthermore,  from 
the  comments  in  Section  3.2,  it  follows  that  SNR  is  also  a  good 
approximation  to  SNR^  for  maximal-connected  sets  of  m-sequences  of  length 
127  with  AO/LSE  phases.  Since  the  expressions  for  SNR  are  considerably 
easier  to  evaluate  than  the  corresponding  exact  expressions,  the  quantity 
SNR  is  a  useful  approximation  for  the  preliminary  design  of  a  DS/SS  system 
operating  over  a  Rician  fading  channel.  Examples  of  this  application  of 
SNR  for  the  AWGN  channel  may  be  found  in  [Pursley  and  Roefs,  1979].  It 
should  be  noted  that  in  the  absence  of  fading,  expressions  (3.80) 
through  (3.83)  all  reduce  to  the  expression  obtained  by  Roefs  and  Pursley 
(1977)  for  SNR  for  the  AWGN  channel  by  setting  y^  =  0  in  (3 . 80) - (3 .83) . 


The  accuracy  of  the  expressions  for  SNR  as  approximations  to  SNR^ 
for  the  sets  of  signature  sequences  studied  above  (i.e.,  maximal-connected 
sets  of  m-sequences  with  LSE/AO  or  AO/LSE  loadings)  does  not  imply  that 
any  arbitrary  collection  of  m-sequences  will  provide  the  same  performance. 
As  an  example  of  this,  (3.68)  was  evaluated  for  N  =  31,  K  =  3,  using  two 
different  collections  of  m-sequences.  One  set  maximizes  the  parameters 
r,  .  and  the  other  minimizes  r,  .  .  These  two  sets  were  found  by  an 

1C  y  2.  tC  y  X 

exhaustive  search  of  sets  of  3  m-sequences  of  length  31  having  every 

possible  phase  [Garber,  1978].  The  comparison  of  (3.82)  with  (3.68)  for 

2 

these  two  sets  is  shown  in  Fig.  15  for  y  =0.05.  Note  that  there  is  a 
significant  difference  (>  0.7  dB)  in  performance  between  the  best  set 
(curve  1)  and  the  worst  set  (curve  3)  of  signature  sequences.  This 

2 

difference  is  larger  for  larger  values  of  S/N q  or  smaller  values  of  y  . 

2 

For  very  large  tf/N^  and  very  small  y  ,  this  difference  is  given  by 

/  \ 
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(3.84: 


which,  for  these  two  sets  of  signature  sequences,  is  2.353  dB. 
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CHAPTER  4 

PROBABILITY  OF  ERROR  BOUNDS  FOR  DS/SSMA  COMMUNICATIONS 
VIA  FADING  CHANNELS 

In  Chapter  1  we  noted  that  for  DS/SSMA  communications  via  an  AWGN 
channel,  three  measures  of  system  performance  are  available:  average 
signal- to-noise  ratio  (SNR),  worst  case  performance,  and  probability  of 
error  (P  ) .  In  the  previous  cha  ter  we  have  analyzed  the  performance  of 
DS/SSMA  communications  via  fading  channels  using  SNR  as  the  system 
performance  measure.  In  this  chapter  we  analyze  the  performance  of 
DS/SSMA  communications  via  fading  channels  using  probability  of  error  at 
the  output  of  a  correlation  receiver  as  the  system  performance  measure. 

The  results  obtained  here  represent  a  generalization  of  the  performance 
analysis  of  [Yao,  1977]  which  considered  only  AWGN  channels. 

Although  we  shall  use  probability  of  error  as  the  system  performance 
measure,  as  noted  in  Section  4.1,  exact  evaluation  of  the  probability  of 
error  at  the  output  of  a  correlation  receiver  for  a  DS/SSMA  system 
operating  over  a  fading  channel  is  very  difficult.  Hence  we  shall 
concentrate  our  effort  on  obtaining  bounds  on  Pg .  Specifically  we  shall 
bound  Pg  through  the  use  of  an  isomorphism  theorem  from  the  theory  of 
moment  spaces  which  provides  relationships  among  arbitrary  moments  of 
a  random  variable. 

Throughout  this  chapter  we  shall  concentrate  on  the  two  singly-spread 
channels  considered  in  Chapter  3:  time-selective  fading  channels  and 
frequency- selective  fading  channels.  With  appropriate  changes  in  the 
equations  that  follow,  the  results  presented  below  may  be  generalized 
to  other  classes  of  WSSUS  fading  channels. 
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4.1.  Moment  Space  Bounds  on  Probability  of  Error 

In  this  section  we  shall  present  the  general  problem  of  obtaining 
moment  space  bounds  on  the  probability  of  error  at  the  output  of  a 
correlation  receiver  for  the  DS/SSMA  system  model  considered  in  Chapter  3. 
In  Sections  4.2  and  4.3  we  shall  consider  specific  examples  of  the  moment 
space  bounding  technique  applied  to  time-selective  and  frequency- selective 
fading  channels  for  second  and  N-th  moment  bounds,  respectively. 

From  Chapter  3,  we  found  that  for  a  DS/SSMA  system  operating  over  a 
doubly-spread  fading  channel,  the  output  of  the  correlation  receiver 
matched  to  the  i-th  user’s  phase-coded  waveform  is  given  by  (see  3.18) 

Z.  =  N.  +  D.  +  F.  +  I.  (4.1) 

i  n.  i  i  x. 

where  and  are  given  by  (3.14)  and  (3.12),  respectively  and  F^  and  1^ 
are  the  real  parts  of  the  complex  quantities  F^  and  1^  defined  by  (3.13) 
and  (3.15),  respectively.  Assuming  that  the  data  bit  b.  n  takes  on  the 
values  {+l,-l}  with  equal  probability,  the  probability  of  error  at  the 
output  of  the  i-th  correlation  at  the  decision  instant  t  =  T  is  given  by 

Pe  =  \  Pr{error|bi  Q=  -l]  +  \  PrCerror^  Q  =  +l] 

«  h  Pr{z.  >  0|b.>0  =  -1)  +  h  Pr[z.  <  0ib.j0  =  +lj 

-  \  El>r[Z.  >  0|b.(t),bi  0  =  -1]3  +  *§  E{Pr[Z.  <  0jb.(t),b.>0  -  +1]) 

(4.2) 

Note  that  in  (4.2)  we  are  simply  evaluating  Pg  conditioned  upon  b^(t) 
and  then  averaging  over  b^(t)  to  obtain  an  expression  for  P  . 


Then  the  first 


For  convenience,  define  h  =  Vp/2  T  and  z  -  F.  +  I. . 
term  inside  the  curly  brackets  of  (4.2)  may  be  rewritten  as 


PrCz.  >  0|b.(t),b.j0  =  -1}  =  E{Q(^)|b.(t),b.}0  =  -1} 


(4.3) 


where  Q(x)  =  1  -  § (x)  and  a  =  V^NqT  is  the  standard  deviation  of 
(see  3.19).  In  (4.3),  the  expectation  over  z  conditioned  upon  b^(t)  and 
b.  n  =  -1  denotes  expectation  over  all  the  random  variables  b  (t) (k  #  i) , 
Tk’  ^k*  anc*  ^k^T,t^  aPPeari-nS  i-n  2*  Since  given  b^(t),  z  is  symmetric 
about  zero,  (4.3)  may  also  be  written  as 


Pr(zi  >  0|bt(t),blj0-  -1}  -  Efe)<^)  +  %QCiSgS)|b1<«.b1>0-  -1]  (4.4) 


In  a  similar  fashion,  the  second  term  inside  the  curly  brackets  of  (4.2) 
may  be  expressed  by 

Pr{z.  <  0|bi(t),bi>0=  +1}  =  EfcQ (^)  +  %QC~ 5-)|bi(t),bij0=  +l}  (4.5) 

From  (4.2),  (4.4),  and  (4.5)  we  find  that 

Pr[error|b.j0=  -l}=  JjE[e£Q(^)+Q(^)  |b.  (t)  ,b.j()=  -l}|b.>()=  -1]  (4.6) 

and 

PrUrror|b.>0=  +l}=  %E[E{Q(^)-HK^)  |bi(t)  ,b.  +l3|b.j0=  +1]  .  (4.7) 

Hence,  if  the  expectations  in  (4.6)  and  (4.7)  can  be  evaluated,  the  probability 
of  error  at  the  output  of  the  i-th  correlation  receiver  may  be  obtained. 


Ill 


B 

K 


Several  approaches  are  available  for  evaluating  (4.2).  The  first  is 
the  direct  approach,  in  which  the  expectations  in  (4.6)  and  (4.7)  are 
directly  evaluated.  In  general,  because  of  the  complicated  form  of  z, 
however,  P^  is  very  difficult  to  evaluate  either  analytically  or 
numerically.  To  give  an  indication  of  the  amount  of  difficulty  required 
to  evaluate  (4.2)  numerically,  a  recent  attempt  to  evaluate  the  expecta¬ 
tions  in  (4.6)  over  only  and  b^  required  1.4  minutes  of  CPU  time  on  a 
Digital  Equipment  Corporation  DEC-10  computer  system  for  a  two-user 
DS/SSMA  system  employing  a  code  sequence  of  length  127  for  a  single 
value  of  <$/Nq.  Clearly  the  analytical  or  numerical  approach  to  evaluating 
(4.2)  becomes  unwieldy  for  larger  numbers  of  users,  longer  code  sequences, 
or  a  large  number  of  <S/Nq  values. 

A  second  approach  to  evaluation  of  (4.2)  is  by  simulation  and  use  of 
the  Monte  Carlo  method.  Such  an  approach  was  taken  by  Orr  (1977)  in 
evaluating  the  probability  of  error  of  a  DS/SSMA  system  using  a  slowlv- 
fading  (nondispersive)  fading  channel.  The  major  drawbacks  of  the 
simulation  method  of  evaluating  Pg  are  that  the  technique  requires  a 
considerable  amount  of  computation  time  and  that  it  provides  very  little 
insight  into  the  design  aspect  of  the  problem. 

A  third  technique  to  evaluating  (4.2)  and  the  one  employed  in  the 
remainder  of  this  chapter  is  the  moment” space  bounding  method.  In  order 
to  understand  the  theoretical  background  for  the  moment- space  bounding 
method,  we  will  now  state,  without  proof,  an  isomorphism  theorem  [Yao,  1977] 
originally  developed  in  the  theory  of  games  ([Dresher,  et.  al.,  1950], 
[Dresher,  1953])  which  provides  relationships  among  arbitrary  moments  of 


a  random  variable. 


Isomorphism  Theorem:  Let  Y  be  a  random  variable  with  a  probability 
distribution  function  G^(y)  defined  over  a  finite  closed  interval 
I  =  [a,b].  Let  k^(y) .k^Cy) , . . . ,k^(y)  be  a  set  of  N  continuous  functions 
defined  on  I.  The  generalized  moment  of  the  random  variable  Y  induced 
by  the  function  k^(y)  is 

mi  =  JIki(y)  dGy(y)  =  EyCk^y)},  i  =  1,...,N.  (4.8) 

We  denote  the  N-th  moment  space  %  by 

b 

^  =  [m  =  (m1,...,mN)€  RN|mi  =  J  (y)dGy(y) , 1  <  i  <  N,Gy€  P(I)}  (4.9) 

where  P(I)  is  the  set  of  probability  distributions  defined  on  I  =  [a,b] 

N 

and  R  denotes  N-dimensional  Euclidean  space.  Then  7/1  is  a  closed,  bounded, 

N 

and  convex  set.  Now  let <3  be  the  curve  £  =  (r^,...,r^)  traced  out  in  R 
by  r^  =  k^(y)  for  y  in  I.  Let  X  be  the  convex  hull  of  (3.  Then 

JC  -  7n.  (4.10) 

A  brief  proof  of  the  above  theorem  may  be  found  in  [Yao  and  Tobin,  1976]; 
a  detailed  discussion  of  the  proof  is  given  in  [Dresher,  et.  al.,  1950]. 

To  demonstrate  how  the  isomorphism  theorem  can  be  used  to  solve 
the  problem  at  hand,  let  N  =  2  and  let  k2(z)  be  equal  to  the  expression 
inside  the  curly  brackets  of  (4.6).  Let  k^(z)  be  some  continuous  function 
of  z  whose  generalized  moment  m^  given  by  (4.8)  may  be  readily  evaluated. 
Now  consider  a  plot  of  k^(z)  versus  k^(z)  and  denote  the  convex  hull  of 
the  resulting  figure  by  JC.  From  the  isomorphism  theorem,  JC  =  >7?  so  that 
knowledge  of  ra^  =  E[k^(z) | b^ (t) ,b^  q  =  -l]  enables  us  to  obtain  bounds  on 
m2  =  E[k2(z) |bi(t) ,bi  0  =  -l}.  But  by  our  choice  of  k2(z), 
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Pr{error|bi(t) ,bi  Q  =  -l}  =  m2 


(4.11) 


so  that  we  have,  in  fact  obtained  bounds  on  the  conditional  probability 
of  error,  conditioned  upon  b^(t)  and  b_^  ^  =  -1.  Repeating  this  process 
for  Pr{error | b. (t) ,b.  ^  =  +l}  and  averaging  the  sum  of  these  two  bounds 
over  b.(t),  from  (4.2)  we  see  that  we  have  bounded  the  probability  of  error 
(P^)  through  use  of  the  isomorphism  theorem. 

In  practice,  we  would  like  to  choose  k^(z)  so  that  the  convex  hull 

of  the  plot  of  k^(z)  versus  k^(z)  is  "thin"  and  hence  the  Pg  bounds 

obtained  are  tight.  Ideally  we  would  like  to  choose  k^(z)  =  k2(z),  in 

which  case  the  convex  hull  K  is  infinitely  thin  and  the  upper  bound 

P^U  equals  the  lower  bound  P  However,  this  would  require 

knowledge  of  tn^  =  E[k2(z)}  which  we  assumed  could  not  be  directly 

evaluated.  Thus  in  choosing  k^(z)  we  must,  in  general,  trade-off  the 

thickness  of  K  with  the  ease  of  evaluating  m^.  Based  on  previous  results 

in  the  published  literature  ([Yan,  1975], [Yao  and  Tobin,  1976], [Yao,  1977]) 

N 

a  promising  choice  for  k^(z)  is  k^(z)  =  z  ,  where  N  is  even.  The  error 
bounds  resulting  from  such  a  choice  are  known  as  N-th  moment  bounds. 

For  a  given  choice  of  k^(z)  and  k2(z),  three  methods  are  available 
for  evaluating  Pe^  and  Pe^  through  the  use  of  the  isomorphism  theorem: 
graphical  techniques ,  numerical  techniques ,  and  analytical  techniques . 

The  graphical  technique  is  attractive  from  the  viewpoint  that  in  order 
to  employ  it,  we  need  only  plot  k2(z)  versus  k^(z)  over  the  range  of  z 
and  complete  the  convex  hull  3C  by  means  of  a  straightedge.  Assuming  that 
m^  is  known,  the  conditional  upper  and  lower  error  bounds,  conditioned 


upon  b^(t),  may  be  read  directly  off  the  graph.  Repeating  this  process 
for  all  data  bit  patterns  b^(t)  of  interest  and  averaging  the  conditional 
error  bounds  over  b.(t),  we  obtain  bounds  on  Pe^  and  Pe^.  However  as 

3. 

pointed  out  in  [Yan,  1975],  for  sufficiently  complex  plots  of  k^(z)  versus 
k^(z),  determination  of  K  by  eye  becomes  quite  difficult.  In  these  cases, 
it  becomes  necessary  to  find  3C  by  numerical  or  analytical  means.  As 
might  be  expected,  the  numerical  technique  involves  finding  the  convex  U 
and  convex  H  regions  of  C-  through  numerical  evaluation  of  the  second 
derivation  of  k^(z)  with  respect  to  k^(z)  over  the  range  of  z.  The 
convex  hull  of  (3-  is  then  found  directly  using  numerical  techniques.  The 
analytical  approach  to  finding  K  is  by  writing  equations  for  K  based  on 
knowledge  of  the  convex  regions  of  C.  From  a  design  viewpoint,  the 
analytical  technique  of  finding  K  is  the  preferred  method  since  the 
resulting  equations  describing  X  are  explicit  functions  of  the  DS/SSMA 
system  parameters  and  the  fading  channel  characteristics.  Because  bounds 
on  PeU  and  PgL  are  determined  by  X  and  m^,  DS/SSMA  system  performance  can 
therefore  be  optimized  for  specific  fading  channels.  In  practice,  as 
noted  in  Sections  4.2  and  4.3,  the  actual  method  used  in  evaluating  X  is 
a  combination  of  both  numerical  and  analytical  techniques. 

The  above  techniques  of  bounding  P^  through  the  use  of  the  isomorphism 
theorem  have  already  been  used  successfully  in  the  evaluation  of  the 
performance  of  several  other  types  of  communication  systems  including 
binary  pulse-amplitude  modulation  (PAM)  systems  with  intersymbol  inter¬ 
ference  ([Yao  and  Tobin,  1975], [Yan,  1975], [Yao  and  Tobin,  1976]), 
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coherent  phase  shift  keyed  (CPSK)  systems  with  cochannel  interference 
[Tobin  and  Yao,  1977],  DS/SSMA  systems  over  AWGN  channels  ([Yao,  1976a], 
[Yao,  1976b], [Yao,  1977]),  and  PSK  systems  employing  bandpass  limiters 

[Yao  and  Milstein,  1978].  In  all  of  these  studies,  either  Nth  moment 
N 

(i.e.,  k^(z)  =  z  )  or  exponential  moment  (i.e.,  k^(z)  =  exp(c(h+z)) 
where  c  is  a  constant  chosen  to  optimize  the  bounds)  bounds  were  used  to 
bound  the  probability  of  error.  In  the  next  two  sections  we  shall  apply 
Nth  moment  space  bounds  to  bound  the  probability  of  error  of  DS/SSMA 
systems  operating  over  specific  classes  of  fading  channels. 

4.2.  Second  Moment  Bounds 

In  this  section  we  shall  evaluate  the  second  moment  space  bound  on 
the  probability  of  error  for  a  DS/SSMA  system  operating  over  either  a 
time-selective  or  a  frequency- selective  fading  channel.  In  the  following 
section  we  shall  generalize  the  results  obtained  to  Nth  moment  bounds. 

In  order  to  clearly  present  our  results  we  shall  first  consider  the  time- 
selective  fading  channel  and  later,  the  frequency-selective  fading  channel. 
At  the  end  of  this  section  we  will  present  various  criteria  that  are  useful 
in  selecting  the  range  of  z  denoted  by  I  in  the  isomorphism  theorem  of  the 
previous  section. 

4.2.1.  Second  Moment  Error  Bounds  for  Time -Selective  Fading  Channels 
To  simplify  the  derivation  of  moment  space  bounds  for  the  probability 
of  error  of  a  DS/SSMA  system  with  a  time -selective  fading  channel,  instead 
of  deriving  results  directly  from  the  doubly-spread  fading  channel  model 
as  was  done  in  Section  3.2,  we  shall  use  the  simpler  (but  equivalent) 


model  of  a  time-selective  fading  channel  developed  in  Section  2.5. 
Hence,  from  (2.52)  and  (2.53),  the  output  of  the  i-th  correlation 
receiver  for  a  DS/SSMA  system  with  a  time-selective  fading  channel  is 


given  by  (4.1),  where  and  are  given  by  (3.12)  and  (3.14), 
respectively  and  =  Re[f\]  and  1^  =  where 

T 

F.-v'isPj'  Y,  P,(t)b  (t)dt 

o 

h  -  -  Rel^  A  Vk,i} 

k=l 

k^i 

I!  =  Re  [I!}  =  ReUkP  El'} 

1  k=l  it>:L 

k^i 

and 

T 

\,i 4  J"0  pk(t-Tk)hk,i<Tkit)e>cp(j'Pk)dt 
*k,i  ’  I0  hk,i(Tk;';)eXpU'Pk)dt- 

In  (4.12)  and  (4.16)  0^(t)  is  a  zero-mean,  unit  variance  complex 
Gaussian  random  process  with  covariance  functions 

EUPk(t)P*(s)}  =  rkCt-s) 


(4.12) 

(4.13) 

(4.14) 

(4.15) 

(4.16) 

(4.17a) 


E .  ( t) p.  (s ) }  =  0  . 


(4.17b) 


From  Section  4.1  we  note  that  in  order  to  determine  second  moment 

bounds  on  P  ,  we  must  first  evaluate  the  moments 
e 

mx  =  E[z2|b.(t)}  (4.18) 

where  z  is  defined  as  in  Section  4.1  and  the  expectation  is  over  all 
b^(t)  for  k  4  i,  cp^  and  in  2  conditioned  upon  all  values  of  b^(t) 

of  interest.  As  in  Chapter  3,  we  shall  model  the  phase  angles,  time 
delays,  and  data  symbols  for  the  k-th  signal  as  mutually  independent 
random  variables  which  are  uniformly  distributed  on  [0,2rr],  [0,T],  and 
{+l,-l},  respectively.  From  (4.12)  through  (4.16)  we  find  that 

E[z2|b.(t)3  =  E(F.2|b.(t)}  +  E{l.2|b.(t)3 

-  E{F.2|b.(t)}  +  E{I.2}  (4.19) 

where  we  have  used  fact  that  1^  is  independent  of  b^(t)  and  that  the  cp^ 
are  independent  random  variables.  Because  F^  and  1^  are  identical  to 
the  terms  F^  and  1^  appearing  in  Section  3.2,  noting  that  1^  is  a  zero 
mean  random  variable,  the  term  e£X^2]  appearing  in  (4.19)  is  identical 
to  Var  1^,  previously  evaluated  in  (3.39).  Furthermore,  using  (1.1) 
in  (4.12),  we  find  that 

5t  -  ■A*  bl>0  . 

Hence  using  (2.10)  and  (4.17),  we  find,  as  in  (3.31),  that 

T  T 

ECFi2 I bi Ct)}  =  f  Y.2  J  J  r  (t-s)dtds 
L  L  0  0 


(4.20) 


which  is  independent  of  fc^(t).  Thus  from  (4.18),  (3.39),  and  (4.21), 

E£z2|bi(t)}  =  e[f.2}  +  e{i.2} 

-  E[z2}  =  ,  (4.22) 

a  result  independent  of  the  choice  of  b^(t). 

Once  the  finite  closed  interval  I  of  z  is  known,  as  required  by  the 
isomorphism  theorem,  we  may  evaluate  second  moment  error  bounds  for  a 
DS/SSMA  system  with  a  time-selective  fading  channel.  Let  the  distortion 
D  be  defined  as  the  maximum  value  of  the  random  variable  z,  i.e. 

D  *  max|z|.  (4.23a) 

z 

By  the  symmetry  of  z  about  zero,  the  interval  I  defined  in  the 
isomorphism  theorem  may  therefore  be  taken  to  be  I  «  [-D,D] .  However 
by  the  choice  of  k^(z)  made  in  Section  4.1  (see  (4.6)),  we  may  reduce  the 
actual  range  I  to  the  interval  I  =  [0,D]  because  of  the  symmetry  of  k^(z) 
about  zero.  In  Section  4.2.3  we  will  present  various  criteria  for  the 
actual  selection  of  D. 

In  Appendix  B,  a  presentation  of  the  methods  used  to  evaluate  the 
convex  hull  of  the  plot  of  k^(z)  versus  k^(z)  over  I  is  given  along  with 
development  of  the  conditional  upper  and  lower  error  bounds,  conditioned 
upon  b^(t),  through  the  use  of  the  isomorphism  theorem.  Once  the 
conditional  error  bounds  are  known,  from  (4.4),  (4.6),  and  (4.7)  we  may 
obtain  bounds  on  P  ^  and  P 
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In  Figure  16  we  have  plotted  upper  and  lower  second  moment  space 
bounds  on  the  probability  of  error  for  a  DS/SSMA  system  with  a  time- 
selective  fading  channel  with  a  triangular  channel  covariance  function  as 
a  function  of  the  bit  energy  to  noise  spectral  density.  In  plotting 
Figure  16  we  have  used  a  value  of  normalized  distortion  of  0.5,  where  the 
normalized  distortion  D'  is  defined  by 

D'  =>  £  .  (4.23b) 

h 

In  addition  to  the  upper  and  lower  error  bounds,  we  have  also  plotted 
the  Gaussian  approximation  to  the  probability  of  error 

PgG  =  1  -  §(SNR  )  (4.24) 

where  SNR^  is  given  by  (3.40)  for  a  time -selective  fading  channel.  The 
Gaussian  approximation  is  exact  when  F^  and  1^  are  Gaussian  random 
variables  [see  the  definition  of  SNR^,  eq.  (3.28)].  Note  that  for  small 
values  of  <S/Nq,  the  bounds  on  Pg  are  quite  tight.  However  for  larger 
values  of  <S/Nq,  the  upper  bound  on  P^  becomes  practically  useless. 
Furthermore  for  larger  values  of  D',  larger  numbers  of  users,  and 
shorter  signature  sequences,  even  the  bounds  for  smaller  values  of 
<?/Nq  become  loose.  This  phenomenon  of  second  order  bounds  was  also 
noted  by  Yao  (1977)  for  DS/SSMA  systems  with  AWGN  channels.  To  improve 
the  probability  of  error  bounds,  in  Section  4.3  we  will  consider  Nth 


moment  bounds  for  n  2  4. 
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4.2.2.  Second  Moment  Error  Boands  for  Frequency-Selective 


Fading  Channels 

Following  the  evaluation  of  second  moment  bounds  for  the  time- 
selective  fading  channel  as  presented  in  the  previous  subsection,  we 
shall  now  evaluate  the  second  moment  error  bounds  for  a  DS/SSMA  system 
employing  a  frequency- selective  fading  channel.  Once  again  we  shall  use 
the  singly-spread  channel  characterization  of  the  frequency-selective 
fading  channel  developed  in  Section  2.6  rather  than  employ  the  equivalent 
method  of  developing  the  frequency-selective  fading  channel  model  from 
the  doubly-spread  channel  model,  as  was  done  in  Section  3.3.  For  a 
DS/SSMA  system  with  a  frequency-selective  fading  channel,  the  output  of 
the  i-th  correlation  receiver  is  once  again  given  by  (4.1),  where  IL 
and  are  given  by  (3.14)  and  (3.12),  respectively.  In  this  case, 
however  F.  and  I.  are  given  by  F.  =  RefF.}  and  I.  -I.  +  i!  where 

i  i  °  J  i  l  ill 


P,(T)f  h  (Tjt)dtdT 


(4.25 


h  -  aeth!  - .c.  Vk,i] 

k=l 


(4.26 


I!  =  Re{l!}  *  ReCv^P  E  I’  } 
1 


(4.27 


00  T 

ik,i  "  I  Pk(T)J  hk>i(Tk  +  T;t)exp(JcPk)dtdT  <4*28 


122 


I^i  -  J*  hk>i(Tk;t)exp(j9k)dt  .  (4.29) 

In  (4.24)  and  (4.27)  Pk(T)  is  a  zero  mean,  unit  energy  (i.e., 

2 

=  \  in  (2.74))  complex  Gaussian  random  process  with  covariance 
func tions 


E[%Pk(T)p*(a)3  =  gk(T)6(T  -  a)  (4.30a) 

and 

EC^k(T)Pk(a)3  =  0  .  (4.30b) 


In  order  to  obtain  moment  space  error  bounds  for  the  frequency- 

selective  fading  channel  case,  we  must  again  obtain  bounds  on  the 

innermost  conditional  expectation  of  (4.6)  and  (4.7)  through  the  evaluation 
A  2 

of  m^  =  E(z  jb^t)}.  For  the  time-selective  fading  channel  case,  it  was 


only.  Hence  conditional  error  bounds  were  needed  for  only  two  data 
bit  values:  b^  ^  *  +1  and  b^  q  =  -1.  From  the  form  of  (4.25),  however, 
for  the  frequency-selective  fading  channel  case,  conditional  error  bounds 
are  needed  for  an  inf ini te  number  of  data  bit  patterns  because  the 
frequency-selective  fading  channel  exhibits  memory.  To  display  this 
dependence  of  F^  on  all  preceeding  and  successive  values  of  data  bits, 
using  (1.1)  and  (3.10)  in  (4.25),  we  find  that 


F. 

a. 


oo  X  ao 

-®  0  1= 


b.  , 


®  (j+l)T 

Jh?  y.  Z  J  P .  (t )  [ b 

j—  jT 


PT(t-T-iT)ai(t-T)ai(t)dtdT 

-j-lRi(T"jT)  +  bi-A(T-jT)]dT, 


(4.31) 
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where  R^(r)  and  R^(r)  are  defined  by  (1.7)  and  (1.8),  respectively. 
Clearly,  we  cannot  evaluate  an  infinite  number  of  moment  space  bounds 
as  required  by  (4.6),  (4.7),  and  (4.31).  Hence  we  must  truncate  the 
number  of  preceeding  and  successive  data  bits  required  to  evaluate  (4.6) 
and  (4.7)  in  some  sensible  manner.  Note  that  this  corresponds  to  the 
problem  of  evaluation  of  the  probability  of  error  of  a  baseband  PAM 
system  with  intersymbol  interference  (e.g.,  see  [Lucky,  et.  al.,  1968]). 

In  Section  3.3  we  assumed  the  selectivity  of  the  channel  is  such  that  only 
the  two  adjacent  symbols  need  be  taken  into  account  in  evaluating  the 
performance  of  a  DS/SSMA  system  with  a  frequency-selective  fading 
channel.  If  we  also  apply  this  constraint  to  (4.31),  (4.31)  reduces  to 


T  T 

Fi  =  Yt  J  PiOO  J  hi(T;t)dtdT 


(4.32) 


-T 

Note  that  the  constraint  on  (3^(t)  implied  by  (4.32)  is  a  stronger 
constraint  than  that  made  in  (3.51),  since  the  former  is  actually  a 
constraint  on  the  range  spread  of  the  channel  [Van  Trees,  1971],  while 
the  latter  is  a  constraint  on  the  selectivity  of  the  channel.  In  a 
similar  manner,  (4.28)  becomes 

A  T  T 

Ik>i  ■  J  Pk(T)  J  hk>i(Tk+T  ;t)exp(jcpk)dtdT  .  (4.33) 

Thus  for  the  frequency-selective  fading  channel  case,  we  need  to 

evaluate  the  conditional  bounds  on  the  probability  of  error  for  8 

different  sets  of  the  bit  pattern  [b.  l5b.  ^jb.  Each  evaluation 

r  l,-1  i,0  i,+l 

2 

of  the  conditional  error  bounds  requires  evaluation  of  m^  =  E{z  |b^(t)3 


i 

r  - 

I 


for  one  of  the  eight  bit  patterns.  Following  the  same  reasoning  given 
for  (4.18),  E{z2|b.(t)3  is  given  by  (4.18),  where  for  the  frequency- 

r  2n 

selective  fading  case,  E(I^  i  is  given  by  Var  I  in  (3.62).  Using  (2.10) 

2 

and  (4.30)  we  find  that  e[f^  |b^(t)}  is  given  by  (3.54).  Thus  we  can 
evaluate  m^  for  each  of  the  eight  sets  of  bit  patterns. 

Using  the  methods  of  Appendix  B,  we  have  evaluated  the  second  moment 
space  bounds  for  the  probability  of  error  of  a  DS/SSMA  system  with  a 
frequency-selective  fading  channel  having  a  triangular  channel  covariance 
function  and  have  plotted  the  results  as  a  function  of  <S/Nq  in  Figure  17. 
In  plotting  these  results  we  have  used  K  =  2  users  employing  signature 
sequences  of  length  127  and  have  assumed  a  normalized  distortion  of 
D'  =0.5.  In  Figure  17  we  have  also  plotted  the  Gaussian  approximation 
(4.24)  to  the  probability  of  error,  where  SNR^  for  the  frequency- 
selective  fading  channel  is  given  by  (3.63).  The  comments  made  in  the 
previous  section  on  the  looseness  of  the  second  order  moment  space  bounds 
apply  here  also.  In  Section  4.3  we  will  consider  higher  order  moment 
space  bounds  in  an  effort  to  tighten  the  bounds  on  the  probability  of 
error . 

4.2.3.  Selection  of  the  Normalized  Distortion 

Up  to  this  point  we  have  not  given  any  consideration  toward  selecting 


the  value  of  the  normalized  distortion  D'  to  be  used  in  evaluating  the 
moment  space  bounds  for  a  DS/SSMA  system  operating  over  a  given  fading 


.26 


channel.  From  Che  isomorphism  theorem  of  Section  4.1  and  the  determina¬ 
tion  of  moment  space  bounds  given  in  Appendix  B,  such  a  selection  of  a 
particular  value  of  D'  is  required  before  moment  space  bounds  can  be 
obtained.  In  this  subsection  three  possible  criteria  are  given  for 
selection  of  D'  . 

At  first  examination,  it  may  appear  that  determination  of  D'  is 
straightforward  from  the  definition  of  D1  given  by  (4.23).  Further 
examination  of  z  however  reveals  that  because  we  have  assumed  t)  to 

be  a  complex  Gaussian  random  process,  the  "tails"  of  Re[p^(T;t)}  are 
infinite  and  hence  D'  *  ®.  The  isomorphism  theorem  requires,  however, 
that  D'  be  finite.  Our  goal  then  is  to  sensibly  choose  a  finite  D' 
in  such  a  manner  that  truncation  of  D'  will  only  marginally  affect  the 
actual  bounds. 

The  first  approach  will  be  termed  the  distribution  function  method. 
Basically  this  method  assume  a  priori  knowledge  of  the  distribution 
function  of  z  in  order  to  sensibly  truncate  the  "tails"  of  z.  One  method 
of  truncation  is  to  simply  assume  that  the  tails  of  z  may  be  neglected  if 

PrCz>D}<YgPe  (4.34) 

where  Pg  is  the  probability  of  error.  For  example,  if  z  is  resumed  to 

2 

be  a  zero  mean  Gaussian  random  variable  with  variance  a  ,  for  on  the 
order  of  10  ,  D  should  be  selected  so  as  to  satisfy 

D  s  4.753423  a  (4.35a) 

or 

D'  £  4.753423  £  . 


(4.35b) 
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Of  course  the  major  drawback  of  this  method  is  that  we  have  also  assumed 
Pg  could  not  be  evaluated  directly,  which  would  be  the  implication  if 
the  distribution  function  of  z  were  known  (see  (4.1)).  Still,  this 
method  merits  consideration  for  cases  in  which  the  moment  space  bounds 
are  exceptionally  tight  and  an  approximating  distribution  function  for  z 
yields  results  for  which  lie  between  the  upper  and  lower  bounds. 

A  second  approach  utilizes  Chebyshev's  inequality  to  obtain  a  finite 
value  for  D'.  By  Chebychev's  inequality  we  have  that 

PrC|z|  >  D}  <  .  (4.36) 

Since  we  have  assumed  that  the  moments  E{zn}  can  be  evaluated,  (4.36) 
may  be  used  in  a  manner  identical  to  (4.34)  to  bound  D'.  As  an  example 
of  this,  for  the  second  moment  case  for  values  of  Pe  on  the  order  of  10  , 
we  choose  D  to  satisfy 

,  2 

10"°  5  (4.37) 

D 

or  D'  ^  10  —  ,  where  a  =  E[z  }.  Assuming  that  higher  order  moments  of 

z  can  be  evaluated,  we  may  then  use  (4.36)  to  further  reduce  this 
initial  choice  of  D' 

A  third  approach  to  choosing  D'  is  based  on  the  actual  DS/SSMA 
system  dynamic  range.  Because  any  physical  system  must  have  a  finite 
dynamic  range,  D  is  automatically  limited  to  the  dynamic  range  of  the 
DS/SSMA  system.  If  we  use  this  approach  to  limit  the  range  of  D',  we 
must  have  prior  knowledge  of  the  actual  system  dynamic  range. 


Ideally  we  would,  of  course,  like  to  choose  D'  as  large  as  possible 
so  that  the  effects  on  the  probability  of  error  bounds  of  truncating  D' 


can  be  minimized.  Fortunately,  as  will  be  seen  in  the  next  section, 

|  for  Nth  moment  space  bounds,  the  probability  of  error  bounds  are 

relatively  insensitive  to  the  choice  of  D1,  for  values  of  D'  greater 
than  a  certain  minimum  value.  For  the  remainder  of  this  chapter  we 
|  will  compute  moment  space  error  bounds  with  the  assumption  that  D'  has 

been  computed  using  one  of  the  three  methods  outlined  above. 

4.3.  Nth  Moment  Bounds 

<  In  this  section  we  shall  consider  Nth  moment  bounds  for  n  ^  2.  Our 

goal  here  is  to  obtain  tighter  moment  space  bounds  on  the  probability  of 
error  of  a  DS/SSMA  system  with  a  fading  channel  than  could  be  obtained 
[  using  the  second  moment  bounds  of  the  previous  section. 

4.3.1.  Fourth  Moment  Bounds 

As  noted  in  Section  4.1,  in  order  to  evaluate  the  Nth  moment  space 
bounds  on  the  probability  of  error  of  a  DS/SSMA  system  with  a  fading 
channel,  we  first  need  to  evaluate  the  moment 

m1  =  ECzn|bi(t)} 

=  E[[F.  +  I.]n|b.(t)}  (4.38} 

for  all  data  bit  patterns  b^(t)  of  interest.  For  the  fourth  moment  space 
bound,  we  therefore  need  to  evaluate 


E[[F.+I.  ]4jb.  (t)  )  =  ECF.‘++4F.JI.+6F.2I.2+4F.I.'5+I.tf|b.  (t)} 
i  l  1  1  3.  11  11  11  i  i 


(4.39 
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In  the  following  we  shall  assume  that  the  fading  processes  P.(T;t) 
and  P^(T;t)  are  independent  for  j  ^  k.  Then  from  (3.13),  (3.15),  and 
(3.16),  (4.39)  reduces  to 

EC[F.+I.]4|b.(t)}  =  E[F.4+6F.2Ii2+Ii4jbi(t)} 

=  E[Fi4|bi(t)}+6E{Fi2|bi(t)}E{li2}+E[li4}  (4.40) 

where  we  have  used  the  fact  that  1^  is  independent  of  b^(t)  and  that 
P^(T;t)  is  a  zero  mean  process.  The  middle  term  of  (4.40)  has  been 
previously  evaluated  in  Sections  4.2.1  and  4.2.2  for  time-selective 
fading  channels  and  frequency-selective  fading  channels,  respectively. 
Hence  in  this  section  we  shall  concentrate  on  evaluating  the  first  and 
last  terms  of  (4.40).  Recognizing  that  F^  =  Re{F\}  and  1^  =  ReDL}, 
in  order  to  evaluate  the  fourth  moments  of  these  quantities,  we  shall 
need  the  following  identity  which  may  be  derived  directly  from  (2.10): 

8Re(w)Re(x)Re(y)Re(z)  =  Re(wxyz)+Re(w*x*yz)+Re(wxy*z) 

+  Re (w*x*y*z)+Re (w*xyz)+Re (wx*yz)+Re(w*xy*z) 

+  Re(wx*y*z),  (4.41) 

for  complex  numbers  w,  x,  y,  z. 

From  (4.40),  (4.41),  (3.15),  and  (3.16),  we  note  that  the  fourth 
moments  of  F.  and  I.  will  be  a  function  of  the  fourth  moments  of  0.  (T;t), 

K 

for  1  <  k  <  K.  To  evaluate  the  fourth  moment  of  P^(T;t)  we  shall  make 
use  of  the  following  theorem  due  to  Reed  (1962)  and  Miller  (1968)  . 


Theorem:  Let  x(t)  be  a  zero-mean,  complex  Gaussian  random  process  with 
covariance  functions 


E[x*(t)x(s)}  =  R(t,s) 
E[x(t)x(s)}  =  0 


(4.42) 

(4.43) 


Suppose  x^  =  x(t^)  f°r  n  €  {l,2,...,N}  are  samples  from  x(t) . 


a.)  If  s  ^  t,  then 


E{x*  x*  ...x*  x  x  ...x  }  =  0 

m.  m-  m  n.  nn  n 

12  s  i  2  t 


where  m^  and  n^  are  integers  from  the  set  {1,2,3, ... ,n] . 


b.)  If  s  =  t,  then 

e{x*  x*  ...x*  x  x 

m^  n^  ^ 


•  x  ]  = 
nt 


£  E{x*  x  }e[x*  x  }...E{x*  x  } 

TT  mTT(l)  ni  mTT(2)  n2  mTT(t)  nt 

where  tt  is  a  permutation  of  the  set  of  integers  [l,2,3 . t} . 


(4.44) 


(4.45) 


A  proof  of  this  theorem  for  stationary  random  processes  was  originally  given 
in  [Reed,  1962];  Miller  (1968)  subsequently  generalized  the  theorem  to 
nonstationary,  non-zero  mean  random  processes.  A  discussion  of  both  of 
these  results  may  be  found  in  [McGee,  1971]. 

Hence  from  (4.44)  and  (4.45),  the  fourth  moments  of  a  zero-mean 
complex  Gaussian  process  8(t)  satisfying  (4.42)  and  (4.43)  are 


ECP*(t)P*(s)P(t,)P(s')}  =  R(t,s)R(t',s)  +  R(t,s,)R(t',s) 
E{p(t)P(s)p(t*)P(s')  -  E{P*(t)P(s)P(t,)P(s')} 

=  E{B*(t)P*(s)P*(t’)e(s’)}  =  0  • 


(4.46) 


(4.47) 
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For  the  time-selective  fading  case,  from  (4.41),  (4.46),  (4.47), 
and  (4.20)  the  conditional  fourth  moment  of  F^  conditioned  upon  b^(t) 
is  given  by 

T  T  T  T 

E £f . ^ | b . C t) }  =  t  J  J  J  r. (t-s)r. (u-v)dtdsdudv} 

o  o  o  o  1  x 

=  3[Var  Ft]2  (4.48) 


where  Var  F^  is  given  by  (3.31)  and  we  have  used  (4.17). 

In  a  similar  fashion,  for  the  frequency-selective  fading  case 
E{F.4|b.(t)3  is  given  by 


EtF  >  (c)}  .  |  p\  4Re{J’  J-  g  (T)g  (T')J  J-  J  J 

“  *  rn  rr»  A  A  n 


T  T 


T  T  T  T 


-T  -T 


0  0  0  0 


(T  >T  >  t,s) f^ (t ' ,t ' ; t' ,s')dtdsdt'ds'dTdT ' } 

=  3[Var  F± ]2  (4.49) 

where  Var  F^  is  given  by  (3.54).  In  deriving  (4.49)  from  (4.32),  we 

have  used  (4.41),  (4.46),  (4.47),  and  (4.30). 

r  4i 

In  order  to  evaluate  E {.I  }  appearing  in  (4.40),  assume,  without  loss 
of  generality,  that  user  1  is  the  i-th  user  (i.e.,  i  =  1).  Letting 

I.  =  i.  +  I!  (4.50) 

ill 

A 

where  1^  and  1^  are  given  by  (4.13)  and  (4.14),  respectively,  for  the 
time -selective  fading  channel  and  by  (4.26)  and  (4.27),  respectively, 
for  the  frequency-selective  fading  channel,  we  find  that 


32 


to 


m 


E[I.4] 


Eli.4  +  4i. 3i!  +  6x.2(ii)2  +  4i. d:)3  +  (x : )4J 


11  11 


1  1 


E[i.4  +  6I.2(ip2  +  (ip4"} 


t  \4-| 


(4.51) 


where  we  have  used  the  zero-mean  condition  on  P^(T;t)  for  the  second  step. 
Noting  that  the  last  term  of  (4.51)  has  been  previously  evaluated  in 
[Yao,  eq.  (14)-(18),  1977],  we  shall  concentrate  on  evaluating  the  first 
two  terms  of  (4.51).  Let 


\,i  s^Vetit,i! 


where  I  .  is  defined  by  (3.16).  Then 

k,  i 

K 


(4.52) 


I.  =  2  I,  .  , 


k=2 


k,i 


(4.53) 


which  follows  from  (3.15).  Using  (4.53)  in  (4.51),  we  see  that 


E[X,4}  =  E([  2  I,. 


K 


k=2 


k,i 


K  K  K  K 

=  E[  2  2  2  2  I  I  I  I  } 

h=2  j=2  k=2  Jj=2  J’  k’ 

K  .  K-l  K  9  9 

-  EC  2  I  4}+6  2  2  E[I  2I  2} 

k=2  ,L  k=2  j=k+l  ’  3,1 

K  .  K-l  K  „  , 

=•  2  e(I,  ;  +  6  2  S  Etl  23e{I.  .3 
k=2  k,i  k=2  j=k+l  k,:L  J’1 


(4.54) 


where  we  have  used  the  independence  of  Pjc(T>c)  f°r  2  <  k  5=  K  and  the 
fourth  moment  properties  of  p,  (t ; t) . 


For  the  time-selective  fading  case,  from  (4.52),  (4.15),  (4.17), 
(4.41),  (4.46),  and  (4.47),  E[l,  .  is  given  by 

K,1 


Etlk,i“J  =  4  p2\  !g  J"Q  Sg  -f0  tk(t'S),:k(u-V)EC£k,k,i<Tk-Tk;':>S) 


T  T  T  T 


,fk,k,i^Tk’Tk;u,V^dtdsdudv’ 


(4.55) 


where  the  expectation  is  over  bk(t)  and  t  .  From  (3.10)  we  see  that  in 
order  to  evaluate  this  expectation,  we  need  to  evaluate 


T  J“0  ak(t-Tk)ak(S-Tk)ak(u-Tk)ak<'"Tk)Etbk<t-Tk)bk(s-Tk)bk(u-Tk)bk(v-Tk>lir 

(4.56) 

Note  that  (4.56)  is  a  function  of  t,  s,  u,  and  v  and  consequently  may 
only  be  evaluated  numerically  for  given  values  of  the  four  parameters. 
Furthermore,  unlike  previous  "moments"  of  a^t)  that  we  have  encountered 
up  to  this  point,  evaluation  of  (4.56)  requires  knowledge  of  the  actual 
sequence  a^t)  rather  than  just  its  aperiodic  correlation  functions. 
Consequently,  the  objective  of  applying  higher  order  moment  space  bounds 
to  bound  Pg  becomes  questionable  unless  (4.56)  can  be  evaluated  in  some 
manner.  One  such  approach  is  to  simply  bound  the  expectation  appearing 
in  (4.55)  so  that  bounds  on  E{.I,  .  }  may  be  obtained.  Noting  that  a 

"worst  case"  bounding  approach  is  to  assume  that  the  expectation  appearing 
in  (4.55)  is  upper  bounded  by  1  and  lower  bounded  by  0,  (4.55)  becomes, 
using  (3.30)  and  (3.31), 


>  -  < 


0,  best  case 


3[Var  F^]  ,  worst  case 


(4.57) 


where  Var  is  given  by  (3.31).  In  a  similar  manner,  the  second  term  of 
(4.51)  may  also  be  bounded,  and  hence,  from  (4.40),  (4.48),  (4.51),  (4.54), 
and  (4.57),  upper  and  lower  bounds  on  m^  appearing  in  (4.38)  may  be 
obtained.  If  we  use  these  upper  and  lower  bounds  on  to  evaluate  upper 
and  lower  moment  space  bounds  on  P^  through  use  of  the  methods  described 
in  Appendix  B,  however,  we  find  that  the  results  obtained  are  actually 
slightly  looser  than  those  obtained  using  the  second  moment  space  bounds 
described  in  the  previous  section.  This  looseness  of  the  higher  order 
bound  may  be  attributed  to  the  crude  bounds  on  the  expectation  appearing 
in  (4.55). 

For  the  frequency- selective  fading  case,  from  (4.52),  (4.30),  (4.33), 

(4.41),  (4.46),  and  (4.47),  e{I.  is  given  by 

k,  i 


- 1  pV  r  rw>.f0 10 10  J-0EfFk,k,i<Tk^-Tk+rjt-s> 


T  T 


T  T  T  T 


•F^  k  ^(T^+cr,T^+CTjU,v)}dtdsdudvdTdCT , 


(4.58) 


As  was  the  case  for  the  time-selective  fading  channel,  the  expectation 
appearing  in  (4.58)  may  not  be  easily  evaluated  analytically.  Consequently, 
once  again  a  simple  bound  on  this  expectation  may  be  made  in  an  effort  to 
bound  and  subsequently  m^.  However,  it  was  found  that  the 

resulting  bounds  on  P  obtained  through  the  use  of  the  bounds  on  m  are 

6  JL 

actually  looser  than  the  second  order  bounds  obtained  in  Section  4.2. 
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4.3.2.  Higher  Moment  Bounds  for  DS/SS  Communications  via 
Fading  Channels 

In  the  previous  subsection  we  noted  that  for  the  cases  in  which 
multiple  users  are  present,  fourth  moment  space  bounds  on  the  probability 
of  error  of  DS/SSMA  communications  via  a  fading  channel  cannot  be  readily 
evaluated.  This  observation  may,  in  fact,  be  generalized  to  Nth  order 
bounds  for  n  >  2,  since  evaluation  of  higher  order  bounds  will  require 
evaluation  of  higher  order  "moments"  of  the  signature  sequences .  For  a 
single  user  system,  however,  the  above  comments  do  not  apply,  since,  as 
noted  in  the  previous  subsection,  we  are  able  to  evaluate  the  fourth 
moments  of  F^  for  both  time-selective  and  frequency-selective  fading 
channels.  In  this  subsection  we  shall  consider  moment  space  bounds  for 
single  user  DS/SS  systems  with  fading  channels. 

To  begin  with,  note  from  (4.20),  (4.31),  and  (4.32)  that  Fi  =  Re[E\} 
is  conditionally  Gaussian,  conditioned  upon  b^(t).  This  follows  from  the 
observation  that  the  real  part  of  a  complex  Gaussian  random  variable  is 
Gaussian.  Furthermore,  throughout  all  of  our  analysis  presented  to  this 
point,  we  have  assumed  that  the  fading  process  P(T;t)  is  independent  of 
the  additive  noise  process  n(t).  Hence  for  a  single  user  DS/SS  system 
with  either  a  time-selective  or  a  frequency-selective  fading  channel,  we 
may  write  an  exact  expression  for  the  probability  of  error.  For  a  single 
user  DS/SS  system,  the  output  of  the  i-th  user's  correlation  receiver 
is  given  by 

Z.  =  N.  +  D.  +  F.  . 

riii. 


(4.59) 


where  we  have  used  (3-54),  (4.2),  (4 . 59) - (4.63) ,  and  where  g.(T)  is  the 
covariance  function  of  the  i-th  user's  fading  process. 

Because  we  can  evaluate  the  probability  of  error  exactly  for  the 
single  user  DS/SS  case,  there  appears  to  be  no  need  to  evaluate  moment 
space  bounds  on  the  probability  of  error.  Nevertheless,  it  is  of  interest 
to  evaluate  moment  space  error  bounds  in  this  instance  since  the  results 
are  useful  in  evaluating  the  performance  of  the  bounds  themselves  for 
related  DS/SSMA  system  problems.  Two  such  examples  of  related  problems 
are  DS/SSMA  systems  with  faded  multiple-access  interference  and  with 
non-faded  multiple-access  interference.  In  the  following  we  will  evaluate 
Nth  order  moment  space  bounds  for  a  DS/SS  system  in  the  presence  of  a 
fading  channel. 

In  the  previous  subsection  we  presented  methods  for  obtaining  fourth 

moment  error  bounds  for  a  DS/SSMA  system  with  faded  multiple-access 

interference.  For  a  single  user  system,  these  same  methods  apply  except 

that  the  interference  term  1^  is  neglected  in  all  of  the  calculations. 

From  (4.38),  it  follows  that  in  order  to  evaluate  Nth  moment  bounds 

on  P  ,  we  need  to  evaluate 
e 

=*  EC[Fi]n|b.(t)}  .  (4.66) 

Two  approaches  are  available  for  evaluating  (4.66).  The  first  requires 
derivation  of  an  expression  for  the  product  of  real  parts  of  complex 
numbers  in  terms  of  sums  of  real  parts  of  products  of  complex  numbers 
and  then  applying  (4.44)  and  (4.45)  to  this  result.  This  is  the  approach 
taken  in  the  previous  subsection.  The  second  approach  makes  use  of  the 


fact  that  is  conditionally  Gaussian,  as  noted  above.  To  use  this 
second  approach,  we  will  need  the  following  result  from  [Wang  and 
Uhlenbeck,  1945]  regarding  moments  of  real  multivariate  Gaussian 
dis  tributions : 


Theorem:  Let  X  =  [x. ,x„,...,x  }  have  a  zero-mean  multivariate  Gaussian 
-  n  1  2  n 

distribution  with  covariance  e[x.x.}  =  R. .. 

i  j  ij 

(i)  Then  for  m  odd: 


x.  x. 

1  12  L3 


}  =  0. 


1 

m 


(4.67) 


(ii)  For  m  even: 


E '  x .  x.  x. 
X2  X2  X3 


•x.  3  =  .  .  R. 

Lm  J1J2  J3J4 


m- 1  m 


(4.68) 


where  the  sum  is  taken  over  all  possible  ways  of  dividing  the  m  points 
into  m/2  combinations  of  pairs.  The  number  of  terms  in  the  sum  is  equal 


to 

1*3  *5  •  *  *  (m-3)  (m-1)  . 


(4.69) 


From  (4.66),  (4.68),  (4.69),  and  the  observation  that  F^  is  conditionally 
Gaussian,  we  have  the  following  result:  for  n  even 

ml  =  EtlF.]n|b.(t)}  =  l-3'5-‘*(n-3)(n-l)[ECF.2tb. (t)}]n/2  .  (4.70) 

Hence  using  (4.70),  (4.2),  and  the  procedures  given  in  Appendix  B,  we 
may  evaluate  the  Nth  moment  space  bounds  on  the  probability  of  error  of 
a  single  user  DS/SS  system  with  a  fading  channel. 
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In  Figure  18  we  have  plotted  upper  and  lower  sixth  moment  space 

bounds  on  the  probability  of  error  as  a  function  of  <S/Nq  for  a  single 

user  DS/SS  system  with  a  time-selective  fading  channel.  In  plotting 

Figure  18  we  have  assumed  that  the  fading  channel  has  a  triangular 

2 

covariance  function  with  correlation  duration  10T„  and  y.  =  0.05.  Also 

C  l 

shown  in  Figure  18  is  the  exact  probability  of  error  for  this  system  which 

was  evaluated  using  (4.64).  For  the  purposes  of  this  plot  we  have 

assumed  a  normalized  distortion  of  D'  =  1.0  in  evaluating  the  moment  space 

bounds.  Note  that  the  moment  space  bounds  on  the  probability  of  error  are 

quite  tight  for  this  particular  set  of  channel  specifications.  From  the  form 

of  (4.64)  we  see  that  Pg  for  the  time -selective  fading  channel  is  independent 

of  the  choice  of  signature  sequences  used.  Following  a  line  of  reasoning 

identical  to  that  given  in  Section  3.4  in  discussing  the  time -selective 

fading  channel,  we  may  conclude  that  our  results  are  valid  also  for  a 

PSK  system  using  a  Rician  fading  channel. 

In  our  discussion  in  Section  3.3  we  noted  the  relative  importance 

of  frequency- selective  fading  channels  as  compared  to  time -selective 

fading  channels  for  SSMA  systems.  For  the  remainder  of  this  subsection 

we  will  present  results  only  for  a  single  user  DS/SS  system  with  a 

frequency-selective  fading  channel.  In  Figures  19  and  20  we  have 

plotted  eighth  moment  bounds  on  P^  versus  <S/Nq  for  a  DS/SS  system  with 

a  frequency- selective  fading  channel  having  a  triangular  channel 

2 

covariance  function  for  two  values  of  y^  •  In  plotting  these  figures 
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we  have  assumed  a  correlation  duration  of  20T^  and  a  normalized  distortion 

of  D'  =  1.0  and  have  used  an  m-sequence  of  length  127  with  its  LSE/AO 

2 

shift  register  loading.  Note  that  as  is  increased  from  0.1  to  0.5 

the  moment  space  bounds  become  less  tight.  This  looseness  in  the  bounds 

2 

is  to  be  expected  for  large  values  of  y^  because  of  the  strong  dependence 

2 

of  m^  on  Var  (see  (4.69))  and  consequently  on  the  parameter  y^ 

(see  (3.31)  and  (3.54)). 


In  order  to  compare  the  tightness  of  various  orders  of  bounds,  in 

Figure  21  we  have  plotted  fourth  moment,  sixth  moment,  and  eighth  moment 

space  bounds  on  the  probability  of  error  as  a  function  of  S/Nq  for  the 

same  DS/SS  system  parameters  and  channel  characteristics  used  in  plotting 

Figure  19.  Generally  speaking,  for  values  of  <S/Nq  less  than  10  dB, 

fourth  moment  bounds  are  tighter  than  the  higher  order  bounds .  For  larger 

values  of  <S/Nq,  though,  only  the  eighth  order  bounds  do  not  "level  off." 

This  phenomenon,  first  observed  by  Yan  (1975),  may  be  simply  explained 

using  the  results  of  Appendix  B.  In  Appendix  B  we  note  that  for  values 

of  <S/Nq  less  than  a  certain  critical  value  SNR^,  the  plot  of  k^(z) 

versus  k.,(z)  is  convex  fl  for  any  value  of  D'.  To  evaluate  P  U  and  P  L 
1  }  e  e 

for  SNR  <  SNR„,  we  need  to  evaluate  ran  =  e(f. n| b . (t) 3 ,  take  the  nth  root 
C  1  l  1  i 

of  m^,  and  then  use  this  result  to  find  the  upper  and  lower  bounds  on  Pg 


by  the  equations 


peu  -  ECk2(mi1/n)3 

k  (D1)  -  k  CO) 

Pe  =  E^k1(D')  -  k2(0)  (ml)  +  k2(0)^ 


(4.71) 


(4.72) 
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2 

of  fourth,  sixth,  and  eighth  moment  bounds  (N  =  127,  Y  =  0.1,  X  =  20,  D’  =  1.0, 

LSE /AO  sequence) . 
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where  the  expectation  is  over  b^(t).  From  (4.69),  however  we  note  that 

1  <  <  ffl5  <  $105  (4.73) 

and  consequently  that 

m  <  m  <  m  <  m..  (4.74) 

2  i4  i6  8 

where 

ml  =  •  (4.75) 

n  i 

Using  (4.74)  in  (4.71)  and  (4.72),  we  see  that  for  SNR  <  SNR  ,  the 
lower  order  bounds  on  are  actually  tighter.  We  may  conclude  that 
unless  a  "leveling-off"  effect  is  observed,  n-2  order  bounds  are  tighter 
than  n  order  bounds. 

In  the  previous  section  we  commented  that  for  fading  channels,  Nth 

moment  space  bounds  are  somewhat  insensitive  to  the  exact  value  of  the 

normalized  distortion  being  used  to  evaluate  the  bounds.  To  illustrate 

this  insensitivity  of  the  bounds  on  D' ,  eighth  order  moment  space  bounds 

were  evaluated  for  a  DS/SS  system  with  a  frequency-selective  fading  channel 

2 

for  a  correlation  duration  of  20Tr,  Y  =0.05,  and  for  three  values  of  D'. 
As  D'  was  varied  from  0.5  to  10.0,  the  error  bounds  varied  by  less  than 
5  X  10  ^  %,  for  the  case  considered.  Thus  in  most  cases  of  interest, 
we  may  apply  any  of  the  methods  suggested  in  the  previous  section  to 
bound  D'  and  still  obtain  approximately  equivalent  error  bounds. 


CHAPTER  5 


SUMMARY  AND  CONCLUSIONS 

This  study  has  investigated  the  problem  of  analyzing  the  performance 
of  a  biphase  DS/SSMA  system  operating  over  a  fading  channel.  The 
receiver  model  used  throughout  this  analysis  has  been  a  correlation 
receiver;  a  receiver  structure  which  is  not  optimal  for  signaling  in  the 
presence  of  multi-user  interference  and  fading  but  one  that  is  typically 
used  in  implementing  DS/SSMA  systems.  Two  measures  of  system  performance 
were  considered:  average  signal- to-noise  ratio  and  average  probability  of 
error.  A  third  system  performance  measure  considered  for  DS/SSMA  systems 
operating  over  AWGN  channels  --  worst  case  performance  --  was  not  treated 
here  since  it  is  easily  shown  that  for  a  channel  undergoing  any  degree 
of  fading,  the  worst  case  probability  of  error  is  always  \  for  equal 
a  priori  probabilities  of  the  transmitted  data  bits. 

The  study  began  with  a  review  of  various  models  of  fading  channels. 
From  physical  and  phenomenological  considerations  we  have  shown  how 
fading  may  be  modeled  as  a  complex  Gaussian  random  process.  With  this 
assumption,  a  general  fading  channel  model  was  developed  which  could  be 
modeled  as  a  time-varying  linear  filter  whose  time-varying  impulse 
response  is  the  fading  process.  Statistical  phenomena  observed 
for  actual  radio  channels  led  us  to  consider  a  specialization  of  the 
general  fading  channel  model  known  as  the  WSSUS  channel.  This  channel 
model  was  characterized  by  the  fact  that  it  is  the  simplest  channel 
model  of  physical  interest  exhibiting  both  time-  and  frequency-selective 
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behavior.  A  review  of  various  techniques  for  characterizing  the  WSSUS 
channel  in  terms  of  observable  effects  was  presented  next.  We  then  noted 
that  by  making  certain  additional  assumptions  about  the  covariance 
functions  of  the  WSSUS  doubly-spread  channel,  the  WSSUS  channel  further 
simplified  to  two  types  of  singly-spread  channels  --  time -selective 
fading  channels  an.'  frequency-selective  fading  channels  --  and  to  a  channel 
known  as  a  nondispersive-f ading  channel.  In  discussing  these  various 
channel  models  we  noted  the  importance  of  assuming  a  specular-plus-diffuse 
fading  channel  model  as  a  typical  fading  channel  model  over  which  DS/SSMA 
communications  could  be  conducted.  As  was  noted  in  Chapter  3,  a  channel 
exhibiting  only  diffuse  fast  fading  renders  a  DS/SSMA  system  useless  due  to 
the  lack  of  phase-reference  needed  for  coherent  communications . 

After  a  review  of  fading  channel  models  we  then  considered  analysis 
of  the  average  signal- to-noise  ratio  of  a  DS/SSMA  system  communicating 
over  a  WSSUS  Rician  fading  channel.  It  was  shown  that  with  such  a 
channel,  the  output  of  the  correlation  receiver  matched  to  the  i-th  user's 
code  waveform  consisted  of  four  components:  the  i-th  user's  direct  (non- 
faded)  signal,  a  faded  version  of  the  i-th  user's  transmitted  signal,  a 
multiple-user  interference  term  consisting  of  direct  and  faded  signal 
components,  and  an  AWGN  term.  Treating  the  data  symbols,  phase  shifts, 
and  time  delays  of  the  other  K-l  users  as  random  variables  with  given 
distribution  functions,  we  then  defined  the  average  signal-to-noise  ratio 
in  terms  of  expectations  of  the  four  components  at  the  output  of  the 
correlation  receiver.  In  general  it  was  found  that  the  resulting 


expectations  could  not  be  evaluated  directly  in  terms  of  the  correlation 
functions  of  the  fading  channel  unless  certain  additional  assumptions 
about  the  covariance  function  of  the  channel  were  made.  Hence,  we 
considered  the  two  singly-spread  channels  --  time-selective  fading 
channels  and  frequency-selective  fading  channels.  For  both  of  these 
subclasses  of  WSSUS  channels  we  were  able  to  evaluate  expressions  for 
the  average  signal-to-noise  ratio  in  terms  of  the  covariance  function  of 
the  fading  channel,  the  continuous -time  partial  crosscorrelation  functions 
of  the  K  users'  code  waveforms,  and  the  AWGN  spectral  density.  For 
specific  channel  covariance  functions,  the  resulting  expressions  become 
functions  of  the  discrete  aperiodic  autocorrelation  functions  of  the 
code  signature  sequences.  A  parameter  of  the  channel  covariance  function 
known  as  the  correlation  duration  was  defined  and  was  then  used  as  the 
independent  variable  in  plotting  SNI^  for  both  time-  and  frequency- 
selective  fading  channels.  It  was  found  that  the  general  trends  of  the 
plots  could  be  explained  by  exploiting  the  duality  existing  between  the 
two  singly-spread  channels  as  well  as  by  noting  the  effect  of  the 
integrator  in  the  correlation  receiver  on  the  performance  of  these  two 
channels.  We  noted  that  most  fading  channels  occurring  in  practice  exhibit 
some  degree  of  frequency-selectivity  and  that  this  phenomenon,  together  with 
the  observation  that  SSMA  systems  occupy  very  large  bandwidths,  led  us  to 
conclude  that  a  frequency-selective  fading  channel  is  a  more  realistic  fading 
channel  model  for  DS/SSMA  systems  than  a  time-selective  fading  channel  model. 
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We  then  examined  the  performance  of  single  user  DS/SS  systems  over 
both  classes  of  singly-spread  channels  and  compared  this  performance 
to  that  of  a  PSK  system  operating  over  the  channel.  We  found  the  DS/SS 
system  performance  identical  to  that  for  a  PSK  system  for  a  time-selective 
fading  channel,  but  the  DS/SS  system  outperformed  the  PSK  system  for  a 
frequency-selective  fading  channel.  This  performance  difference  in  the 
latter  channel  may  be  attributed  to  the  frequency-diversity  achieved 
by  the  DS/SS  system  through  the  use  of  a  signalling  set  with  an  improved 
autocorrelation  function. 

By  randomizing  the  code  signature  sequences  in  the  DS/SSMA  system 
we  were  able  to  obtain  approximate  expressions  for  the  average  SNR  which 
were  independent  of  the  actual  choice  of  signature  sequences  used.  These 
expressions  are  therefore  useful  in  the  preliminary  design  of  a  DS/SSMA 
system  using  a  fading  channel.  It  was  found  that  for  maximal  connected 
sets  of  m-sequences  employing  either  LSE/AO  or  AO/LSE  shift  register 
loadings,  the  difference  between  the  actual  SNR  and  the  approximate  value 
of  SNR  was  negligible. 

After  considering  SNR  as  a  performance  measure  we  next  considered 
average  probability  of  error  as  a  measure  of  system  performance.  At  the 
outset,  it  was  noted  that  due  to  the  complicated  form  of  the  resulting 
expression  for  the  average  probability  of  error,  it  would  be  desirable  to 
obtain  bounds  on  the  probability  of  error.  Because  of  the  success  of 
its  use  in  bounding  error  probability  in  other  types  of  digital  communi¬ 
cation  systems,  the  moment  space  bounding  technique  was  used  to  obtain 
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bounds  for  DS/SSMA  communication  via  fading  channels.  Based  on  results 
in  the  published  literature,  we  chose  to  use  Nth  moment  space  bounds, 
bounds  which  require  evaluation  of  the  Nth  moment  of  the  sum  of  the  faded 
version  of  the  i-th  user's  transmitted  signal  plus  the  multiple-access 
interference.  For  the  second  moment  bound,  we  were  able  to  evaluate  the 
probability  of  error  bounds  using  results  from  the  analysis  of  the  signal- 
to-noise  ratio  at  the  output  of  the  i-th  user's  correlation  receiver.  It 
was  noted  that,  in  general,  the  second  moment  error  bounds  were  too  loose 
for  practical  application.  Hence,  we  then  considered  higher  order  moment 
bounds  in  an  effort  to  obtain  tighter  error  bounds.  For  n  >  2,  however, 
we  were  unable  to  directly  evaluate  the  nth  moment  of  the  sum  of  the  faded 
version  of  the  i-th  user's  transmitted  signal  plus  the  multiple-access 
interference  and  consequently  we  were  unable  to  obtain  tighter  error  bounds. 
It  should  be  noted  that  much  of  the  difficulty  in  applying  the  moment  space 
bounding  technique  to  DS/SSMA  communications  via  fading  channels  is  due  to 
the  relative  complexity  of  the  singly-spread  fading  channel  model  being  used 
throughout  this  study.  We  conjecture  that  for  simpler  channel  models,  such 
as  the  recently  developed  three  path  fading  channel  model  mentioned  in 
Chapter  2,  it  would  be  possible  to  evaluate  moment  space  bounds  on  the 
probability  of  error.  Of  course,  such  a  channel  model  is  less  general  than 
the  singly- spread  channel  model. 

In  contrast  to  the  results  obtained  for  the  multiple-user  case,  we 
were  able  to  directly  evaluate  the  probability  of  error  for  a  single- 
user  DS/spread-spectrum  system  for  either  a  time-selective  or  a  frequency- 
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selective  fading  channel.  To  investigate  the  performance  of  Nth  moment 
space  error  bounds  for  fading  channels,  we  considered  the  single-user 
DS/spread-spectrum  system  and  evaluated  the  fourth,  sixth,  and  eighth  moment 
error  bounds.  We  found  that  for  small  values  of  SNR,  the  lower  moment  error 
bounds  yield  the  tightest  bounds,  but  for  larger  values  of  SNR  and  D',  only 
the  higher  moment  bounds  yield  good  results. 
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APPENDIX  A 

AVERAGE  SIGNAL- TO-NOISE  RATIO  RESULTS  FOR  LARGE  VALUES  OF  THE 
FADING  POWER  TRANSMISSION  COEFFICIENT 
In  Chapter  2  we  made  mention,  in  passing,  of  the  recent  development 
of  a  three-path  fading  channel  model  for  line-of- sight  microwave  channels 
([Rummler,  1978],  [Rummler,  1979]).  In  the  literature  describing  this 
channel  model,  mention  is  made  of  several  actual  channels  exhibiting 

2 

fading  to  the  extent  that  the  fading  power  transmission  coefficient  V  for 

2 

these  channels  is  on  the  order  of  y  =  0.5.  In  Chapter  3,  we  have  presented 

2 

signal-to-noise  ratio  data  for  channels  in  which  y  ^  0.2,  primarily  for  the 
reason  that  it  is  believed  that  the  effects  of  the  synchronization  sub¬ 
system  on  the  communication  system  performance  are  no  longer  negligible 
2 

for  values  of  y  >  0.2.  For  line-of-sight  applications,  however,  other 

techniques  (e.g.  atomic  clocks,  auxiliary  synchronization  channels)  may  be 

used  to  obtain  synchronization  information  so  that  reliable  communications 

may  be  achieved.  In  this  appendix  we  present  additional  SNR  data  for 

2 

channels  in  which  y  >0.2. 

In  Figure  A1  we  have  plotted  SNR  as  a  function  of  the  correlation 
duration  of  the  channel  for  a  time-selective  fading  channel  with  a 
triangular  covariance  function.  In  plotting  Figure  A1  we  have  used  a 
maximal  connected  set  of  six  m-sequences  in  their  LSE/AO  phases  (see 
Chapter  3  for  details) .  In  Figure  A2  we  have  plotted  the  corresponding 
results  for  a  frequency-selective  fading  channel  with  a  triangular 
covariance  function.  We  see  from  these  two  plots  that  the  increased 
amount  of  fading  simply  reduces  the  average  SNR. 
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In  Figure  A3  we  have  plotted  SNR  versus  correlation  duration  for 
single  user  DS /spread-spectrum  and  PSK  systems  operating  over  a  frequency- 

selective  fading  channel  with  a  triangular  correlation  function  for  two 

2  2 
values  of  y  .  Note  that  for  these  larger  values  of  y  ,  the  SNR  at  the 

output  of  the  PSK  system  is  so  small  as  to  make  the  PSK  system  virtually 

useless,  while  for  values  of  correlation  duration  greater  than  60  the 

degredation  of  the  DS /spread-spectrum  system  from  the  ideal  AWGN  channel 

signal-to-noise  ratio  is  less  than  0.5  dB. 


lannel 


CORRELATION  DURATION  (multiples  of  t.) 

s.  correlation  duration  for  a  frequency-selective  fading  channel  with  a  triangular 
iance  function  for  single-user  DS/spread-spectrum  and  biphase  PSK  systems 
127,  K  =  6,  <f/Nn  =  10dB,  LSE/AQ  sequences). 
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APPENDIX  B 

Nth  MOMENT  SPACE  BOUNDS  ON  THE  PROBABILITY  OF  ERROR 
B.l  Introduction 

In  Chapter  4  we  discussed  the  isomorphism  theorem  regarding  relation¬ 
ships  between  arbitrary  moments  of  random  variables  and  showed  how  this 
theorem  could  be  used  to  obtain  bounds  on  the  probability  of  error  of  a 
DS/SSMA  system  with  a  Rician  fading  channel.  In  using  this  theorem,  we 
derived  the  Nth  moments  of  the  fading  terms  and  then  proceeded  to  plot  the 
probability  of  error  bounds  directly.  In  this  appendix  we  will  derive 
analytical  expressions  for  bounds  on  Pg  from  knowledge  of  the  Nth  moment 
of  the  fading  terms  and  the  normalized  distortion  D'.  The  resulting 
expressions  for  probability  of  error  bounds  were  used  in  Chapter  4  to  plot 
the  actual  error  bounds.  Many  of  the  results  to  be  presented  were  previously 
obtained  by  Yan  (1975)  but  have  been  refined  here  for  use  with  fading 
channels . 

To  begin  with,  instead  of  concerning  ourselves  with  the  details  of  a 
DS/SSMA  system  with  a  fading  channel,  we  shall  use  a  simpler  (but  more  general) 
model  of  a  digital  communication  system  with  interference  and  noise  in  deriving 
the  moment  space  error  bounds.  The  modifications  necessary  to  treat  DS/SSMA 
systems  with  fading  channels  are  straightforward  and  will  not  be  discussed 
here  (see  Chapter  4).  For  a  transmitted  data  bit  a^,  assumed  to  take  on 
value  {+l,-l}  with  equal  probability,  the  received  signal  at  the  output  of  the 
receiver  at  the  decision  instant  t  ■  T  is  given  by 

y  -  aQh  +  z  +  n  (B.l) 

where  h  may  be  thought  of  as  the  response  of  the  channel,  sampled  at  time 
t  »  T  due  to  the  transmission  of  data  bit  a^,  z  is  the  interference  component, 


The  inter- 


and  n  is  an  additive  Gaussian  noise  component  with  variance  a  • 
ference  component  is  assumed  to  have  an  unknown  distribution  function  with 
finite  distortion  D  such  that 

Izl  £  D.  (B.2) 


The  probability  of  error  at  the  output  of  the  receiver  is  given  by 
=  |Pr{y  *  0|aQ  =  +1}  +  |Pr{y  a  o|  aQ  =  -l) 


=  e{£q(^)  +  £q(^)1 


(B.3) 


where  the  expectation  is  over  the  interference  term  z. 


To  derive  Nth  moment  space  bounds  on  the  probability  of  error  of  this 


digital  communication  system,  we  will  let  k^(z)  =  z  and  let  k^(z)  = 
aQ(~^)  +  2Q(~)<  Furthermore,  we  will  assume  that  the  moment 
x  =  Efk^z)}  is  known.  Noting  that  m^  *  ECk^Cz)}  =  Pg ,  we  may  use  the 


m 


isomorphism  theorem  of  Chapter  4  to  derive  bounds  on  Pg.  In  order  to  apply 


this  theorem,  however,  we  first  need  to  evaluate  the  convex  hull  X  of  the 


curve  C*  of  k^(z)  versus  k^(z).  From  calculus,  we  have  the  following  results 


(  i)  A  function  f(x)  is  convex  U  iff  f"(x)  2:  0. 
(ii)  A  function  f(x)  is  convex  0  iff  f"(x)  £  0. 


Hence,  in  order  to  evaluate  the  convex  hull  of  C-,  we  need  to  evaluate  the 


first  and  second  derivatives  of  r(z)  =  k^Ck^Cz)).  Simple  differentiation 


shows  that 


r '  (z) 


fh-z)2  (h+z)2 
2  2 

1 _  e  2g  -e  2j  A  R  u£z). 

2no>/2TT  zn  1  zn  1 


(B.4) 


_  K  f  zu'  (z)  -  (n-l)u(z) 


where 


K  =  [Inojlri]  ~ 1 


(B.6  ) 


and 


(h-z)2  (h+z)2 

2  ~  2 
,  ,  L  la  la 

u(z)  =  e  -e 


(B.  7) 


If  we  collect  identical  terms  appearing  inside  the  brackets  of  (B.5),  r"(z) 
becomes 


r "  (z)  =  P(z)F(z) 


(B.  8) 


where 


_  (h+z)' 


P(2)  =  _ S _  e  2a4 

^  '  2  2n-l  e 

no  z 


(B.9) 


and 


2zh 

‘2 


F (z)  =*  [zh  +  z2+ (n-l)a2]  +  e  a  [zh-z2- (n-l)a2]  . 


(B. 10) 


An  alternative  form  of  writing  F(z)  is  to  expand  the  exponential  term  into 
a  Taylor  series  expansion,  multiply  by  the  appropriate  powers  of  z  appearing 
within  the  square  brackets  and  then  collect  identical  terms.  If  we  do  this, 
(B.10)  becomes 


F(z)  =  Fx(z)  +  F2(z)  (B.  11) 

where 

,2  2 

F^z)  -  2 (2-n)  [zh  +  ]  (B.12) 


,  2zh  .k 

(  o  ) 


F7(2)  *z  Z 
k-1 


k: 


C(k,n,^) 


C(k,n,|)  -  )-l>  +  \  • 

a 


(B. 13) 
(B.14) 


As  we  will  see  in  the  next  two  sections,  depending  upon  the  value  of  n,  either 
(B.8)  and  (B.10)  or  (B.8)  and  (B.ll)  is  more  convenient  to  use  ii  valuating 
the  convexity  of 
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Finally  we  will  need  to  know  the  following  results  on  the  limiting 
values  of  r' (z)  and  r"(z): 

r'(z)  2  0  for  z  2  0  (B.15) 

lim  r '  (z)  =  0  (B*  1-&) 

z->  a 


lim  r"(z) 
z—  0 


2  2  2 
(  \  -3]  .*h 


.n  =  2 


L 


GO 


n  2:  3  . 


(B. 17) 


All  of  these  results  may  be  derived  from  the  defining  equations  for  r'(z) 

and  r"(z)  and  L'Hopital's  rule.  In  the  next  two  sections  we  will  consider 

in  detail  the  convexity  of  C>  for  n  =  2  and  for  n  2  3,  respectively.  We 

will  then  use  these  results  to  obtain  analytical  expressions  for  bounds 

on  P  in  terms  of  the  Nth  moment  of  z  and  the  distortion  D. 
e 


B.2  Second  Moment  Error  Bounds 

In  the  previous  section  we  noted  that  we  could  write  an  expression  for 
r"(z)  as  either  P(z)F(z)  or  P(z)[F^(z)  +  F2(z)J.  Because  of  the  form  of 
F^(z)  it  is  more  expedient  to  use  the  latter  expression  for  r"(z)  for  n  =  2. 
Hence  for  n  *  2,  we  have  that 


r2  (2) 


P(z)F2(z) 

/ 2zh  \  k 

2*  \  c2  '  h 

P(z)z  I  — a.-r  C  (k , 2 ,  J  ). 

k-1  k-  a 


(B. 18) 


Because  k^(z)  and  k^(z)  are  both  symmetric  about  z  =  0,  in  applying  the 
isomorphism  theorem  to  bound  Pe»  we  need  only  let  I  =  [0,D],  i.e.,  only 
positive  values  of  z  need  be  considered  in  evaluating  the  convexity  of  C-. 
From  (B.9)  and  (B.18),  we  see  that  for  n  =  2  and  z  2  0,  the  sign  of  r"(z) 
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ta 


m 


is  completely  determined  by  C(k,2,  ^  ),  which  is  a  monotonically  decreasing 
function  of  k  with  a  maximum  occurring  at  k  =  1.  Hence  we  note  that 


r 


C(k>2>  ff)  =  (-1)  +  2  f(k+l)(k+2) 

O' 


£  0,  -  £  S,  Vk 

o 


(B. 19) 


s  0,  —  2  S  for  some  k  <  k 
O  1 


where  k^  is  the  smallest  integer  such  that 


2k, 


(-D  +  o  ( 


1  1  (k1+l)(k1+2) 


]  S  0. 


(B.20) 


From  (B.18)  and  (B.19)  we  have  the  result  that  if  —  <  J3  ,  C-  is 


convex  fl  for  all  values  of  z.  For  —  2  using  (B.18)  and  (B.20)  we  see 
that  for  small  values  of  z,  C-  is  convex  U.  For  larger  values  of  z,  however, 
and  for  k  ^  k^  C  is  negative  and  hence  r"(z)  may  be  less  than  zero.  From 
(B. 17)  we  see  that,  in  fact,  for  large  enough  z  and  for  —  2  ,  r"(z)  is 

negative.  Consequently  for  ^  2  ,  C,  is  initially  convex  U  and  then 

becomes  convex  fl. 

In  Figures  Bl  and  B2  we  have  plotted  Q/  for  —  <  Jz  and  —  2  Jz,  respectively 
for  n  =  2  and  for  arbitrary  values  of  D.  In  plotting  these  two  figures  we 
have  used  the  results  of  the  previous  paragraph  to  determine  the  convexity  of 
G-.  In  Figure  B2  we  have  labeled  four  points  on  O  to  aid  in  determining  the 
moment  space  probability  of  error  bounds.  The  point  z^  is  defined  to  be  the 
point  satisfying  the  equation 


r"(2l)  =  0  • 

The  point  z ^  is  the  point  on  C-  such  that  a  line  drawn  from  the  point 
(k^(0) .k^ (0) )  is  tangent  to  C  at  z^,  i.e.,  z^  is  the  solution  to  the 
equation 


(B.21) 


n  r 


s  ( 


Ei 


h  o 


« ■ 

y-  - 


h+zr 


-{Q(— — )  +  QC 


2z 
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h-z,  u 

7^)  ‘  2<^  " 

(h-z2)2  (h+z2)2 


2  a 


4 a  z2 


-e 


2a 


}  • 


(B.22) 


The  point  z^  is  defined  to  be  the  point  on  C/  such  that  a  line  drawn  from 
(k^(D) ,k2 (D))  is  tangent  to  <3  at  z^ ,  i.e.,  z 3  is  the  solution  to  the  equation 


,h-Ds 


h+z. 


h-z. 


2 (D  -z3) 


(h-a3r 

17“ 


4a 


-e 


(h+z3)‘ 

2a2 


}  • 


(B.23) 


Finally  point  z^  is  the  point  on  (2-  which  intersects  with  a  line  tangent 
at  (k^ (0) , k2 (0) ) ,  i.e.,  z^  is  the  solution  to  the  equation 


h+z 


h-z 


1  u  i  to  <  ii  1 

T-{q(— )  +  Q(-r-%  -2Q<-)5 - ; 


2a 


2z 


2av^n  a 


(B.24) 


In  (B.24)  we  have  used  the  fact  that  r'(0)  is  equal  to  the  right-hand  side 
of  (B.24),  which  follows  from  (B.4)  and  a  simple  application  of  L'Hopital's 
rule  for  derivatives. 

Note  that  all  of  the  four  points  z^  through  z^  are  solutions  to  single 
variable  nonlinear  equations.  Hence  they  may  be  solved  by  any  number  of 
methods  including  the  bisection  method,  the  f alse-position  method,  or  Newton's 
method.  Once  these  points  are  determined,  we  may  directly  write  down 
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equations  for  the  moment  space  lower  and  upper  bounds  on  the  probability  of 
error  through  the  use  of  the  isomorphism  theorem.  In  Tables  B1  and  B2  we 
summarize  these  equations  for  lower  and  upper  bounds,  respectively.  Note 
that  in  general,  these  bounds  are  functions  of  h/a,  D,  and  m^. 

B.3  Nth  Moment  Error  Bounds 

We  may  easily  extend  the  results  for  second  moment  error  bounds  to  Nth 
moment  error  bounds  for  n  even  and  n  s  2.  However  since  F^(z)  is  no  longer 
identically  zero  for  n  £  2,  we  should  not  expect  that  the  convexity  of  C-  be 
identical  to  that  of  the  n  *  2  case.  In  the  following  we  will  show  that  this 
is,  in  fact,  the  case. 

For  large  values  of  z  we  will  use  the  expression  r"(z)  *  P(z)F(z)  in 

determining  the  convexity  of  fl-.  From  the  observations  made  in  the  previous 

section  we  need  only  consider  z  2  0  in  determining  the  convex  regions  of  C-. 

From  (B. 15)  and  (B.16)  we  may  conclude  that  for  large  z,  r"(z)  is  negative 

and  hence  C-  is  convex  fl  for  large  z.  For  smaller  values  of  z  we  will  use 

the  expression  r"(z)  =  P(z)[F^(z)  +  F^(z)]  to  determine  the  convexity  of  <3. 

Since  we  are  only  considering  non-negative  values  of  z,  from  (B.9), 

P(z)  £  0  and  thus  the  sign  of  r"(z)  is  determined  by  the  sign  of  the  sum 

F^(z)  +  For  n  >  2,  F^(z)  is  always  negative  and  from  (B.14)  we  see 

that  at  least  the  first  2n-4  terms  of  F2(z)  are  negative.  Thus  for  small 

enough  z  and  — ,  the  sum  F^(z)  is  negative  and  hence  r"(z)  is  negative  for 

h 

small  values  of  z.  For  larger  values  of  z  and  —  it  may  be  possible  that 

o 

the  k>  2n-4  terms  of  F2(z)  will  dominate  the  k  S  2n-4  terms  of  F2(z)  as 

-il  as  F^(z)  so  that  F^(z)  +  F2(z)  >  0  for  large  values  of  z  and  We 

h 

may  conclude  from  this  that,  depending  upon  the  value  of  — ,  C,  is  either 
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a  a 
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2 

m-z  h+z_  h-z, 

•ij  2]  3)  +  Q<0  3» 

D  -»3 

D 
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1 

1 

Table  Bl.  Second  moment  lower  bounds  on  P 
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1 

5>^ 
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D£  z2 
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2fQ(ff  )+Q<c  >  '  2QVlml  +  Q(a} 
2z2 
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|[Q(— — h  +  Q( — -)] 

Vm^s  z2 
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Z2S  DS  Z4 
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Table  B2.  Second  moment  upper  bounds  on  P^. 
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convex  fl  for  all  values  of  z  or  <3  is  initially  convex  fl,  then  convex  U, 

and  finally  convex  0  as  z  varies  from  zero  to  infinity. 

In  order  to  determine  the  critical  value  of  —  (SNR  )  at  which  the 

a  c 

second  case  occurs,  we  will  again  consider  the  equation  r"(z)  =  P(z)F(z). 
Note  that  the  second  case  results  from  the  first  case  when  the  maximum 
of  F(z)  is  greater  than  zero.  In  equation  form  this  condition  is  equivalent 
to  requiring  that 


dF(z  ) 

—r. — =  0  when  F(z  )  =  0 
dz  '  p 


(B.25) 


where  z^  is  the  point  at  which  the  convexity  of  C-  changes.  From  (B.10)  we 


dF(z  ) 


'  l»  +  2V  +  th-2zp+p  (2ph-^  -  (n-l)<j2)]e  ° 


0  (B.26 ) 


F  (zp)  *  [zph2  +  z2  +  (n-l)cy2]  +  e 


[z  h-z2  -  (n-l)CT2] 
P  P 


=  0  .  (B.27) 


Solving  for  e  in  (B.26)  and  (B.27)  and  equating  the  two  results,  we  find 
that  Zp  is  the  solution  to  the  quartic  equation 


Zp  4-  [(2n-l)tf2  -  h2]z2  +  a^(n-l)(n-2)  =  0 


(B.28) 


Note  that  (B.28)  is  actually  a  quadratic  equation  in  the  variable  Zp  and 
hence  has  a  solution  given  by 


%  =  2^[  2  ‘  <2n_1)J  ±  K  2  ‘  -4  (n- 1)  (n-2 )  ]*  .  (B.29) 


In  order  for  a  solution  to  exist,  z  must  be  real,  i.e. 

P 


L4  -  - 

^  a  (2n-l)  +  (n-1) (n-2)  .  (B.30) 

a 

Note  that  (B.30)  is  just  the  condition  on  —  for  fl  to  be  both  convex  fl  and 

convex  U.  Alternatively,  if  we  define  SNR  as 

c 


SNR 

c 


4 


[  (2n-l)  +  V4(n-l)(n-2) 


(B. 31) 


tl 

then  fl  is  convex  fl  for  all  values  of  z  if  —  <  SNR  :  otherwise  fl  has  both 

a  c’ 

convex  fl  and  convex  U  regions. 

In  Figure  B3  we  have  plotted  fl  for  various  values  of  ^  and  D.  Using 
these  plots  and  the  value  of  m^  =  E{k^(z)},  we  may  easily  write  down  equations 
for  lower  and  upper  bounds  on  and  have  done  so  in  Tables  B4  and  B5, 
respectively.  Note  that  the  error  bounds  are  dependent  upon  the  values  of 
the  points  z^  through  z^  depicted  in  Figure  B3.  The  points  z^  and  z^  are  the 
solutions  to  the  equation 

r"(z)  -  0  (B.32) 


for  ^  ^  SNRc«  The  point  z^  is  the  point  on  fl  such  that  a  line  drawn  from 
the  point  (k^ (0) ,k^ (0) )  is  tangent  to  fl  at  Zy  i.e.,  z^  is  the  solution  to 
the  equation 


h+z. 


“  IQ(- 


)  +  Q(~~)  -  2 Q  (■“) }  =  r'(z.) 


(B.33) 
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Table  B3.  Nth  moment  lower  bounds  on  P  . 
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Table  B4.  Nth  moment  upper  bounds  on  Pg. 
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The  point  z,  is  defined  to  be  the  point  on  Q,  such  that  the  line  drawn  from 
0 

the  point  (k^ (0) , k^ (0) )  through  the  point  (k^(z^)  ,k^  (z^)  )  intersects  C>  at 

z, ,  i.e.,  z  is  the  solution  to  the  equation 
0  o 


h+z  h-z 

r'(z3)  -  —  (QC^)  +Q(-TS->  ' 

2Z6 


(B. 34) 


The  point  z ^  is  the  point  on  i  such  that  a  line  drawn  from  (k^ (D) , k^ (D) ) 
is  tangent  at  z^,  i.e.,  z^  is  the  solution  to 


h+z. 


r'(z7)  = 


Q(  *“)  +  Q(~“)  -  Q(-“0  -  Q(- 


h-z. 


2[D 


n. 

z7 


(B. 35) 


Finally  the  points  z^  and  z,.  are  defined  to  be  those  points  on  (3  such 
that  a  line  drawn  from  z^  on  C-  to  z^  on  C  is  tangent  at  both  points  and  that 
<3  does  not  intersect  the  connecting  line,  i.e.,  z^  and  z^  are  the  solutions 
to  the  equations 


and 


r'  (zL)  =  r’  (z5) 


r'  (zL) 


h+z  h-z.  h+z.  h-z. 

Q(T^)  +  -«— > 


(B.36) 


(B. 37) 


Because  (B.32)  through  (B.35)  are  one  variable  nonlinear  equations, 
they  may  be  solved  using  any  of  the  methods  listed  in  the  previous  section. 
Equations  (B.36)  and  (B.37)  however  are  two  variable  nonlinear  equations. 
One  method  of  solving  this  set  of  equations  is  to  use  a  Fibonacci  search 
([Aoki,  1971],  [Wilde,  1964])  to  search  over  z^  while  holding  z^  fixed, 
where  we  are  trying  to  minimize  the  quantity 
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